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PREFACE 



This book is based on courses MA381 and EC3080, taught at Trinity College 
Dublin since 1992. 

Comments on content and presentation in the present draft are welcome for the 
benefit of future generations of students. 

An electronic version of this book (in ETgX) is available on the World Wide 
Web at http : //pwaldron . bess . ted . ie/teaching/ma381 /notes/ 
although it may not always be the current version. 

The book is not intended as a substitute for students' own lecture notes. In particu- 
lar, many examples and diagrams are omitted and some material may be presented 
in a different sequence from year to year. 

In recent years, mathematics graduates have been increasingly expected to have 
additional skills in practical subjects such as economics and finance, while eco- 
nomics graduates have been expected to have an increasingly strong grounding in 
mathematics. The increasing need for those working in economics and finance to 
have a strong grounding in mathematics has been highlighted by such layman's 
guides as ?, ?, ? (adapted from ?) and ?. In the light of these trends, the present 
book is aimed at advanced undergraduate students of either mathematics or eco- 
nomics who wish to branch out into the other subject. 

The present version lacks supporting materials in Mathematica or Maple, such as 
are provided with competing works like ?. 

Before starting to work through this book, mathematics students should think 
about the nature, subject matter and scientific methodology of economics while 
economics students should think about the nature, subject matter and scientific 
methodology of mathematics. The following sections briefly address these ques- 
tions from the perspective of the outsider. 



What Is Economics? 

This section will consist of a brief verbal introduction to economics for mathe- 
maticians and an outline of the course. 
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x PREFACE 

What is economics? 

1 . B asic microeconomics is about the allocation of wealth or expenditure among| 
different physical goods. This gives us relative prices. 

2. Basic finance is about the allocation of expenditure across two or more time 
periods. This gives us the term structure of interest rates. 

3. The next step is the allocation of expenditure across (a finite number or a 
continuum of) states of nature. This gives us rates of return on risky assets, 
which are random variables. 

Then we can try to combine 2 and 3. 
Finally we can try to combine 1 and 2 and 3. 
Thus finance is just a subset of micoreconomics. 
What do consumers do? 

They maximise 'utility' given a budget constraint, based on prices and income. 
What do firms do? 

They maximise profits, given technological constraints (and input and output prices).| 
Microeconomics is ultimately the theory of the determination of prices by the in- 
teraction of all these decisions: all agents simultaneously maximise their objective 
functions subject to market clearing conditions. 

What is Mathematics? 

This section will have all the stuff about logic and proof and so on moved into it. 
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NOTATION 



Throughout the book, x etc. will denote points of 3r™ for n > 1 and x etc. will 

denote points of 3? or of an arbitrary vector or metric space X. X will generally 

denote a matrix. 

Readers should be familiar with the symbols V and 3 and with the expressions 

'such that' and 'subject to' and also with their meaning and use, in particular 

with the importance of presenting the parts of a definition in the correct order 

and with the process of proving a theorem by arguing from the assumptions to the 

conclusions. Proof by contradiction and proof by contrapositive are also assumed. 

There is a book on proofs by Solow which should be referred to here. 1 

3?+ = j x G 3l N : Xi > 0, i = 1, . . . ,N\ is used to denote the non-negative or- 

thant of %t N , and 3?^ + = [x G ^. N : %i > 0, i = 1, . . . , JVJ used to denote the 

positive orthant. 

T is the symbol which will be used to denote the transpose of a vector or a matrix. 



'insert appropriate discussion of all these topics here. 
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Chapter 1 
LINEAR ALGEBRA 



1.1 Introduction 

[To be written.] 

1.2 Systems of Linear Equations and Matrices 

Why are we interested in solving simultaneous equations? 
We often have to find a point which satisfies more than one equation simultane- 
ously, for example when finding equilibrium price and quantity given supply and 
demand functions. 

• To be an equilibrium, the point (Q, P) must lie on both the supply and 
demand curves. 

• Now both supply and demand curves can be plotted on the same diagram 
and the point(s) of intersection will be the equilibrium (equilibria): 



• solving for equilibrium price and quantity is just one of many examples of 
the simultaneous equations problem 

• The ISLM model is another example which we will soon consider at length. 

• We will usually have many relationships between many economic variables 
defining equilibrium. 

The first approach to simultaneous equations is the equation counting approach: 

Revised: December 2, 1998 



1.2. SYSTEMS OF LINEAR EQ UATIONS AND MATRICES 



a rough rule of thumb is that we need the same number of equations as 
unknowns 

this is neither necessary nor sufficient for existence of a unique solution, 
e.g. 

- fewer equations than unknowns, unique solution: 

x 2 + y 2 = =^ x = 0,2/ = 

- same number of equations and unknowns but no solution (dependent 
equations): 

x + y = 1 
x + y = 2 

- more equations than unknowns, unique solution: 

x = y 

x + y = 2 
x-2y + l = 



=>- x = l, y = l 

Now consider the geometric representation of the simultaneous equation problem, 
in both the generic and linear cases: 

• two curves in the coordinate plane can intersect in 0, 1 or more points 

• two surfaces in 3D coordinate space typically intersect in a curve 

• three surfaces in 3D coordinate space can intersect in 0, 1 or more points 

• a more precise theory is needed 

There are three types of elementary row operations which can be performed on a 
system of simultaneous equations without changing the solution(s): 

1 . Add or subtract a multiple of one equation to or from another equation 

2. Multiply a particular equation by a non-zero constant 

3. Interchange two equations 
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CHAPTER 1. LINEAR ALGEBRA 5 

Note that each of these operations is reversible (invertible). 

Our strategy, roughly equating to Gaussian elimination involves using elementary 

row operations to perform the following steps: 

1 . (a) Eliminate the first variable from all except the first equation 

(b) Eliminate the second variable from all except the first two equations 

(c) Eliminate the third variable from all except the first three equations 

(d) &c. 

2. We end up with only one variable in the last equation, which is easily solved. 

3. Then we can substitute this solution in the second last equation and solve 
for the second last variable, and so on. 

4. Check your solution! ! 

Now, let us concentrate on simultaneous linear equations: 
(2x2 EXAMPLE) 

x + y = 2 (1.2.1) 

2y-x = 7 (1.2.2) 

• Draw a picture 

• Use the Gaussian elimination method instead of the following 

• Solve for x in terms of y 

x = 2-2/ 

x = 2y — 7 

• Eliminate x 

2-y = 2y-7 



Find y 



Sy = 9 
V = 3 



Find x from either equation: 



x = 2-y = 2-3 = -1 
x = 2y -7 = 6-7= -1 
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6 1.2. SYSTEMS OF LINEAR EQUATIONS AND MATRICES 

SIMULTANEOUS LINEAR EQUATIONS (3x3 EXAMPLE) 
• Consider the general 3D picture . . . 



Example: 



x + 2y + Sz = 6 (1.2.3) 

4x + 5y + 6z = 15 (1.2.4) 

7x + 8y + 10z = 25 (1.2.5) 



Solve one equation (1.2.3) for x in terms of y and z: 

x = 6 — 2y — 3z 

Eliminate x from the other two equations: 

4 (6 - 2y - 32) + 5y + 6z = 15 
7 (6 - 2y - 32) + 8y + 102 = 25 

What remains is a 2 x 2 system: 

— 3y — 62 = —9 
-6y-llz = -17 



Solve each equation for y: 







y = 3-22 






17 11 






y 6 6 


liminal 


:ey: 


17 11 

3-22= 2 

6 6 


ind 2: 




1 1 

6 ~ 6 Z 

2 = 1 



• Hence y — 1 and x = 1. 
Revised: December 2, 1998 



CHAPTER 1. LINEAR ALGEBRA 7 

1.3 Matrix Operations 

We motivate the need for matrix algebra by using it as a shorthand for writing 
systems of linear equations, such as those considered above. 

• The steps taken to solve simultaneous linear equations involve only the co- 
efficients so we can use the following shorthand to represent the system of 
equations used in our example: 

This is called a matrix, i.e. — a rectangular array of numbers. 

• We use the concept of the elementary matrix to summarise the elementary 
row operations carried out in solving the original equations: 

(Go through the whole solution step by step again.) 

• Now the rules are 

- Working column by column from left to right, change all the below 
diagonal elements of the matrix to zeroes 

- Working row by row from bottom to top, change the right of diagonal 
elements to and the diagonal elements to 1 

- Read off the solution from the last column. 

• Or we can reorder the steps to give the Gaussian elimination method: 
column by column everywhere. 



1.4 Matrix Arithmetic 

• Two n x m matrices can be added and subtracted element by element. 

• There are three notations for the general 3x3 system of simultaneous linear 
equations: 



1. 'Scalar' notation: 



an^i + ai 2 x 2 + ai3^3 = h 
a 2 \Xi + a 22 x 2 + 023^3 = h 
a 31 Xi + a 32 x 2 + a 33 x 3 = b 3 
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1.4. MATRIX ARITHMETIC 



2. 'Vector' notation without factorisation: 



a n xi + a u x 2 + oti3^3 \ 


(h 


a 2 \Xi + a 2 2X 2 + ^23^3 


= h 


a 3 ixi + a 32 x 2 + 033X3 / 


\h 



3. 'Vector' notation with factorisation: 

On «12 0l3 \ { %i 

021 «22 ^23 J I ^2 
031 °32 a 33 / \ ^3 




It follows that: 



On a 12 «13 
a 21 a 22 ^23 
031 a 32 033 




anXi + a 12 x 2 + 013X3 \ 

021^1 + 0-22X 2 + ^23^3 
031^1 + O32X2 + O33X3 / 



• From this we can deduce the general multiplication rules: 

The r/th element of the matrix product AB is the product of the 
ith row of A and the jth column of B. 

A row and column can only be multiplied if they are the same 
'length.' 

In that case, their product is the sum of the products of corre- 
sponding elements. 

Two matrices can only be multiplied if the number of columns 
(i.e. the row lengths) in the first equals the number of rows (i.e. 
the column lengths) in the second. 

• The scalar product of two vectors in 5R n is the matrix product of one written 
as a row vector (lxn matrix) and the other written as a column vector (nxl 
matrix). 

• This is independent of which is written as a row and which is written as a 
column. 

So we have C = AB if and only if Cij = Y, k = l n aikbkj- 
Note that multiplication is associative but not commutative. 
Other binary matrix operations are addition and subtraction. 
Addition is associative and commutative. Subtraction is neither. 
Matrices can also be multiplied by scalars. 
Both multiplications are distributive over addition. 
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CHAPTER 1. LINEAR ALGEBRA 

We now move on to unary operations. 

The additive and multiplicative identity matrices are respectively and I n = 

—A and A -1 are the corresponding inverse. Only non-singular matrices have 

multiplicative inverses. 

Finally, we can interpret matrices in terms of linear transformations. 

• The product of an m x n matrix and an n x p matrix is an m x p matrix. 

• The product of an m x n matrix and annx 1 matrix (vector) is an m x 1 
matrix (vector). 

• So every m x n matrix, A, defines a function, known as a linear transfor- 
mation, 

T A : 9fT ->■ K m : x i-> Ax, 

which maps n— dimensional vectors to m— dimensional vectors. 

• In particular, annxn square matrix defines a linear transformation mapping 
n— dimensional vectors to n— dimensional vectors. 

• The system of n simultaneous linear equations in n unknowns 

Ax = b 

has a unique solution Vb if and only if the corresponding linear transfor- 
mation Tj^ is an invertible or bijective function: A is then said to be an 
invertible matrix. 

A matrix has an inverse if and only the corresponding linear transformation is an 
invertible function: 

• Suppose Ax = b does not have a unique solution. Say it has two distinct 
solutions, Xi and x 2 (x x ^ x 2 ): 

Axi = b 
Ax 2 = b 

This is the same thing as saying that the linear transformation T^ is not 
injective, as it maps both xi and x 2 to the same image. 

• Then whenever x is a solution of Ax = b: 

A (x + xi — x 2 ) = Ax + Axi — Ax 2 
= b + b - b 
= b, 

so x + xi — x 2 is another, different, solution to Ax = b. 
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10 1.4. MATRIX ARITHMETIC 

• So uniqueness of solution is determined by invertibility of the coefficient 
matrix A independent of the right hand side vector b. 

• If A is not invertible, then there will be multiple solutions for some values 
of b and no solutions for other values of b. 

So far, we have seen two notations for solving a system of simultaneous linear 
equations, both using elementary row operations. 

1. We applied the method to scalar equations (in x, y and z). 

2. We then applied it to the augmented matrix (A b) which was reduced to the 
augmented matrix (I x). 

Now we introduce a third notation. 

3. Each step above (about six of them depending on how things simplify) 
amounted to premultiplying the augmented matrix by an elementary ma- 
trix, say 

EgE^EgEsE! (A b) = (I x) . (1.4.1) 

Picking out the first 3 columns on each side: 

E 6 E 5 E4E 3 E2EiA = I. (1-4.2) 

We define 

A" 1 = E 6 E 5 E4E 3 E 2 Ei. (1.4.3) 

And we can use Gaussian elimination in turn to solve for each of the columns! 
of the inverse, or to solve for the whole thing at once. 

Lots of properties of inverses are listed in MJH's notes (p.A7?). 

The transpose is A T , sometimes denoted A' or A*. 

A matrix is symmetric if it is its own transpose; skewsymmetric if A T = — A. 

Note that (A T ) _1 = (A" 1 ) 7 . 

Lots of strange things can happen in matrix arithmetic. 

We can have AB = even if A ^ and B^O. 

Definition 1.4.1 orthogonal rows/columns 

Definition 1.4.2 idempotent matrix A 2 = A 

Definition 1.4.3 orthogonal 1 matrix A T = A -1 . 

Definition 1.4.4 partitioned matrices 

Definition 1.4.5 determinants 

Definition 1.4.6 diagonal, triangular and scalar matrices 



'This is what ? calls something that it seems more natural to call an orthonormal matrix. 
Revised: December 2, 1998 



CHAPTER 1. LINEAR ALGEBRA 11 

1.5 Vectors and Vector Spaces 

Definition 1.5.1 A vector is just an n x 1 matrix. 

The Cartesian product of n sets is just the set of ordered n-tuples where the ith 
component of each n-tuple is an element of the ith set. 

The ordered n-tuple (#i, #2, ■ ■ ■ , x n ) is identified with the nxl column vector 
( Xi \ 
x 2 

V Xn J 

Look at pictures of points in !ft 2 and 3ft 3 and think about extensions to 3ft n . 
Another geometric interpretation is to say that a vector is an entity which has both 
magnitude and direction, while a scalar is a quantity that has magnitude only. 

Definition 1.5.2 A real (or Euclidean) vector space is a set (of vectors) in which 
addition and scalar multiplication (i.e. by real numbers) are defined and satisfy 
the following axioms: 

1. copy axioms from simms 131 notes p. 1 

There are vector spaces over other fields, such as the complex numbers. 

Other examples are function spaces, matrix spaces. 

On some vector spaces, we also have the notion of a dot product or scalar product: 

U.V = U V 

The Euclidean norm of u is 

\/u.u =|| u II . 

A unit vector is defined in the obvious way . . . unit norm. 

The distance between two vectors is just || u — v || . 

There are lots of interesting properties of the dot product (MJH's theorem 2). 

We can calculate the angle between two vectors using a geometric proof based on 

the cosine rule. 



||v-u|| 2 = (v-u).(v-u) (1.5.1) 

= || v || 2 + || u || 2 -2v.u (1.5.2) 

= || v || 2 + || u || 2 -2 || v (HI u || cos6» (1.5.3) 

Two vectors are orthogonal if and only if the angle between them is zero. 
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12 1.6. LINEAR INDEPENDENCE 

A subspace is a subset of a vector space which is closed under addition and scalar 
multiplication. 

For example, consider row space, column space, solution space, orthogonal com- 
plement. 



1.6 Linear Independence 

Definition 1.6.1 The vectors x l7 x 2 , x 3 , . . . , x r G !R n are linearly independent if 
and only if 



y^ ajXj = =>- «j = OVi. 



«=i 



Otherwise, they are linearly dependent. 

Give examples of each, plus the standard basis. 

If r > n, then the vectors must be linearly dependent. 

If the vectors are orthonormal, then they must be linearly independent. 



1.7 Bases and Dimension 

A basis for a vector space is a set of vectors which are linearly independent and 

which span or generate the entire space. 

Consider the standard bases in 5ft 2 and 5ft n . 

Any two non-collinear vectors in 3? 2 form a basis. 

A linearly independent spanning set is a basis for the subspace which it generates. 

Proof of the next result requires stuff that has not yet been covered. 

If a basis has n elements then any set of more than n elements is linearly dependent 

and any set of less than n elements doesn't span. 

Or something like that. 

Definition 1.7.1 The dimension of a vector space is the (unique) number of vec- 
tors in a basis. The dimension of the vector space {0} is zero. 

Definition 1.7.2 Orthogonal complement 

Decomposition into subspace and its orthogonal complement. 
Revised: December 2, 1998 



CHAPTER 1. LINEAR ALGEBRA 13 

1.8 Rank 

Definition 1.8.1 

The row space of an m x n matrix A is the vector subspace ofW 1 generated by 
the m rows of A. 

The row rank of a matrix is the dimension of its row space. 

The column space ofanmxn matrix A is the vector subspace of$i m generated 
by the n columns of A. 

The column rank of a matrix is the dimension of its column space. 

Theorem 1.8.1 The row space and the column space of any matrix have the same 
dimension. 

Proof The idea of the proof is that performing elementary row operations on a 
matrix does not change either the row rank or the column rank of the matrix. 
Using a procedure similar to Gaussian elimination, every matrix can be reduced to 
a matrix in reduced row echelon form (a partitioned matrix with an identity matrix 
in the top left corner, anything in the top right corner, and zeroes in the bottom left 
and bottom right corner). 

By inspection, it is clear that the row rank and column rank of such a matrix are 
equal to each other and to the dimension of the identity matrix in the top left 
corner. 

In fact, elementary row operations do not even change the row space of the matrix. 
They clearly do change the column space of a matrix, but not the column rank as 
we shall now see. 

If A and B are row equivalent matrices, then the equations Ax = and Bx = 
have the same solution space. 

If a subset of columns of A are linearly dependent, then the solution space does 
contain a vector in which the corresponding entries are nonzero and all other en- 
tries are zero. 

Similarly, if a subset of columns of A are linearly independent, then the solution 
space does not contain a vector in which the corresponding entries are nonzero 
and all other entries are zero. 

The first result implies that the corresponding columns or B are also linearly de- 
pendent. 

The second result implies that the corresponding columns of B are also linearly 
independent. 
It follows that the dimension of the column space is the same for both matrices. 

Q.E.D. 
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14 1.9. EIGENVALUES AND EIGENVECTORS 

Definition 1.8.2 rank 

Definition 1.8.3 solution space, null space or kernel 

Theorem 1.8.2 dimension of row space + dimension of null space = number of 
columns 

The solution space of the system means the solution space of the homogenous 

equation Ax = 0. 

The non-homogenous equation Ax = b may or may not have solutions. 

System is consistent iff rhs is in column space of A and there is a solution. 

Such a solution is called a particular solution. 

A general solution is obtained by adding to some particular solution a generic 

element of the solution space. 

Previously, solving a system of linear equations was something we only did with 

non-singular square systems. 

Now, we can solve any system by describing the solution space. 

1.9 Eigenvalues and Eigenvectors 

Definition 1.9.1 eigenvalues and eigenvectors and X-eigenspaces 

Compute eigenvalues using det (A — AI) = 0. So some matrices with real entries 
can have complex eigenvalues. 

Real symmetric matrix has real eigenvalues. Prove using complex conjugate ar- 
gument. 

Given an eigenvalue, the corresponding eigenvector is the solution to a singular 
matrix equation, so one free parameter (at least). 
Often it is useful to specify unit eigenvectors. 

Eigenvectors of a real symmetric matrix corresponding to different eigenvalues 
are orthogonal (orthonormal if we normalise them). 
So we can diagonalize a symmetric matrix in the following sense: 
If the columns of P are orthonormal eigenvectors of A, and A is the matrix with 
the corresponding eigenvalues along its leading diagonal, then AP = PA so 
P _1 AP = A = P T AP as P is an orthogonal matrix. 

In fact, all we need to be able to diagonalise in this way is for A to have n linearly 
independent eigenvectors. 
P _1 AP and A are said to be similar matrices. 

Two similar matrices share lots of properties: determinants and eigenvalues in 
particular. Easy to show this. 
But eigenvectors are different. 
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1.10 Quadratic Forms 

A quadratic form is 

1.11 Symmetric Matrices 



Symmetric matrices have a number of special properties 



1.12 Definite Matrices 



Definition 1.12.1 Ann x n square matrix A is said to be 

positive definite ^=^ x T Ax >0Vx6F, x^O 
positive semi-definite 4=4> x T Ax > Vx G 3? n 

negative definite -<=>• x T Ax < Vx G !ft ra , x^O 
negative semi-definite <=^> x T Ax < Vx G W 1 



Some texts may require that the matrix also be symmetric, but this is not essential 
and sometimes looking at the definiteness of non-symmetric matrices is relevant. 
If P is an invertible n x n square matrix and A is any n x n square matrix, then 
A is positive/negative (semi-)definite if and only if P _1 AP is. 
In particular, the definiteness of a symmetric matrix can be determined by check- 
ing the signs of its eigenvalues. 

Other checks involve looking at the signs of the elements on the leading diagonal. 
Definite matrices are non- singular and singular matrices can not be definite. 
The commonest use of positive definite matrices is as the variance-covariance 
matrices of random variables. Since 



Vij = COV [fi, Tj] = Cov [fj, fi] 



(1.12.1) 



and 





N N 






T Vw 


i=l j=l 


(1.12.2) 






' N N 






= Cov 


i=l j=l 


(1.12.3) 




N 






= Var[] 


r w^i] > o 


(1.12.4) 



1=1 



a variance-covariance matrix must be real, symmetric and positive semi-definite. 
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In Theorem 3.2.4, it will be seen that the definiteness of a matrix is also an essen- 
tial idea in the theory of convex functions. 

We will also need later the fact that the inverse of a positive (negative) definite 
matrix (in particular, of a variance-covariance matrix) is positive (negative) defi- 
nite. 

Semi-definite matrices which are not definite have a zero eigenvalue and therefore 
are singular. 
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Chapter 2 

VECTOR CALCULUS 



2.1 Introduction 

[To be written.] 

2.2 Basic Topology 

The aim of this section is to provide sufficient introduction to topology to motivate 
the definitions of continuity of functions and correspondences in the next section, 
but no more. 

• A metric space is a non-empty set X equipped with a metric, i.e. a function 

d : X x X ^ [0,oo) such that 

1. d(x, y) = -<==>■ x = y. 

2. d(x, y) = d(y, x) Vx, y G X. 

3. The triangular inequality: 

d(x, z) + d(z, y) > d(x, y) Vx, y,z G X. 

• An open ball is a subset of a metric space, X, of the form 

Bt(x) = {y G X : d(y,x) < e}. 

• A subset A of a metric space is open 

Vx G A, Be > such that 5e(x) C A. 
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• A is closed -<=>■ X — A is open. (Note that many sets are neither open nor 
closed.) 

• A neighbourhood of x E X is an open set containing x. 

Definition 2.2.1 Let X = !R n . A C X is compact <=^ A is both closed and 
bounded (i.e. 3x, e such that A C Be(x)). 

We need to formally define the interior of a set before stating the separating theo- 
rem: 

Definition 2.2.2 If Z is a subset of a metric space X, then the interior of Z, 
denoted int Z, is defined by 

z E int Z -<=>• Be (z) C Z for some e > 0. 



2.3 Vector-valued Functions and Functions of Sev- 
eral Variables 

Definition 2.3.1 A function (or map) / : X — > Y from a domain X to a co- 
domain Y is a rule which assigns to each element of X a unique element ofY. 

Definition 2.3.2 A correspondence / : X — ► Y from a domain X to a co-domain 
Y is a rule which assigns to each element of X a non-empty subset ofY. 

Definition 2.3.3 The range of the function f : X —> Y is the set f(X) = {f(x) E 

Y:xeX}. 

Definition 2.3.4 The function f : X — ► Y is injective (one-to-one) 

f(x) = f(x') =>■ x = x'. 

Definition 2.3.5 The function f : X — ► Y is surjective (onto) 

f(X) = Y 

Definition 2.3.6 The function f : X —> Y is bijective (or invertiblej 

it is both injective and surjective. 
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Note that if /: X -»■ Y and A C X and 5 C F, then 

/(A)E{/(,):,eA}cy 

and 

/- 1 (S) = {xGX:/(x)G J B}CX 

Definition 2.3.7 A vector- valued function is a function whose co-domain is a sub- 
set of a vector space, say $l N . Such a function has N component functions. 

Definition 2.3.8 A function of several variables is a function whose domain is a 
subset of a vector space. 

Definition 2.3.9 The function f: X — ► Y (X C 3ft n , Y C dt) approaches the limit 

y* as x ^ x* 

Ve > 0, 35 > s.t. || x - x* ||< S => |/(x) - y*)| < e. 

This is usually denoted 

lim /(x) =y*. 

Definition 2.3.10 The function f:X^Y(XQ$t n ,YC$t)is continuous at x* 

Ve > 0, 35 > s.t. || x - x* ||< 5 =^ |/(x) - /(x*)| < e. 
This definition just says that / is continuous provided that 

lim/(x) = /(x*). 

x^x* 

? discusses various alternative but equivalent definitions of continuity. 

Definition 2.3.11 The function f : X — ► Y is continuous 

it is continuous at every point of its domain. 

We will say that a vector- valued function is continuous if and only if each of its 
component functions is continuous. 

The notion of continuity of a function described above is probably familiar from 
earlier courses. Its extension to the notion of continuity of a correspondence, 
however, while fundamental to consumer theory, general equilibrium theory and 
much of microeconomics, is probably not. In particular, we will meet it again in 
Theorem 3.5.4. The interested reader is referred to ? for further details. 
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Definition 2.3.12 I. The correspondence f:X -> Y (X C 9ft n , Y C 3?) is 
upper hemi-continuous (u.h.c.) at x* 

/or every open set N containing the set /(x*), 35 > 5.?. || x — x* ||< 
5=>f(x)CN. 

(Upper hemi-continuity basically means that the graph of the correspon- 
dence is a closed and connected set.) 

2. The correspondence f: X —> Y (X C 3ft n , Y C 3ft) « lower hemi-continuous| 
(l.h.c.) at x* 

/or every open set N intersecting the set /(x*), 35 > 5.?. || x — x* ||< 
8 =>■ /(x) intersects N. 

3. The correspondence f:X^Y(XQyt n ,YQyt) is continuous (at x*) 

it is both upper hemi-continuous and lower hemi-continuous (at x* ) 
(There are a couple of pictures from ? to illustrate these definitions.) 

2.4 Partial and Total Derivatives 

Definition 2.4.1 The (total) derivative or Jacobean of a real-valued function ofN 
variables is the N -dimensional row vector of its partial derivatives. The Jacobean 
of a vector-valued function with values in 3ft M is an M x N matrix of partial 
derivatives whose jth row is the Jacobean of the jth component function. 

Definition 2.4.2 The gradient of a real-valued function is the transpose of its Ja- 
cobean. 

Definition 2.4.3 A function is said to be differentiable at x if all its partial deriva- 
tives exist at x. 

Definition 2.4.4 The function f : X — ► Y is differentiable 
it is differentiable at every point of its domain 

Definition 2.4.5 The Hessian matrix of a real-valued function is the (usually sym- 
metric) square matrix of its second order partial derivatives. 
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Note that if /: ffl 1 — ► !R, then, strictly speaking, the second derivative (Hessian) of 
/ is the derivative of the vector- valued function 

(f'f :W l ^W l :x^ (f(x)) T . 

Students always need to be warned about the differences in notation between the 
case of n = 1 and the case of n > 1. Statements and shorthands that make sense 
in univariate calculus must be modified for multivariate calculus. 

2.5 The Chain Rule and Product Rule 

Theorem 2.5.1 (The Chain Rule) Let g: K n -► W n and /: 9ft m -► K p be contin- 
uously differentiable functions and let h: 3? n — > !R P be defined by 

/»(x) = /0/(x)). 

Then 

j^ = n^x))^(x). 

pxn pxm mxn 

Proof This is easily shown using the Chain Rule for partial derivatives. 

Q.E.D. 

One of the most common applications of the Chain Rule is the following: 

Let g: 3?" — ► 3? m and /: 3? m+n — ► W be continuously differentiable functions, let 

x G 3? n , and define h: K n -»■ 8? by: 

/ l (x) = /(^(x),x). 

The univariate Chain Rule can then be used to calculate |^- (x) in terms of partial 
derivatives of / and g for i = 1, . . . , p and j = 1, . . . , n: 

aii m flfi q k m+n ap o fc 

|-(x) = g^( 9 (x,,x)|-(x, + i i; +i |-( 9( x,,x)|_ ( x). (2 ,,) 

Note that 

— (x) = ^^{ 1 iffc = ^' 
9xj J 1 otherwise ' 

which is known as the Kronecker Delta. Thus all but one of the terms in the second 
summation in (2.5.1) vanishes, giving: 

Oil 1 m Qfi Q k Qfi 



dxj f^y dxk "' dx-i dx 
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2.5. THE CHAIN RULE AND PRODUCT RULE 



Stacking these scalar equations in matrix form and factoring yields: 



/gw 



£w\ 









tf£ 






8/ p 
9a; >r 



a^ 



(a(x),x)/ \^(x 

bfo(*)>*)\ 



£w\ 



£w/ 



ii 



dx m + 



dP 



dP 



(2.5.2) 



Now, by partitioning the total derivative of / as 



r 



px{m+n) 



DJ(-) D x f(. 



pxm 



pxn 



(2.5.3) 



we can use (2.5.2) to write out the total derivative hi (x) as a product of partitioned 
matrices: 

h' (x) = D e f (g (x) , x) g' (x) + D x f (g (x) , x) . (2.5.4) 

Theorem 2.5.2 (Product Rule for Vector Calculus) The multivariate Product Rule\ 
comes in two versions: 

1. Let f, g: K m -»■ 9£ n and define h: W n ->■ 3? £y 

MX)E(J^5(X). 
lxl lxn nxl 



Then 



^=(^/w+c/wry(x). 

lxm lxn nxra lxn nxm 



2. Le* /: K m -► 3* and g: R m -► 3? n and <fe/me ft: K m -► 3? n Zry 

ft(x) = /(x)a(x). 



nxl lxl nxl 



Then 



nxm nxl lxm lxl nxm 



Proof This is easily shown using the Product Rule from univariate calculus to 
calculate the relevant partial derivatives and then stacking the results in matrix 

form - Q.E.D. 
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2.6 The Implicit Function Theorem 

Theorem 2.6.1 (Implicit Function Theorem) Let g: 3ft n — > !ft m , where m < n. 
Consider the system ofm scalar equations in n variables, g (x*) = m . 
Partition the n- dimensional vector x as (y, z) where y = (xi, £2, • • • , x m) is Tri- 
dimensional and z = (x m +i, x m +2, ■ ■ ■ ,x n ) is (n — m) -dimensional. Similarly, 
partition the total derivative of g at x* as 

g' (x*) = [D y g D z g] g 

(m x n) (m x m) (?72 x (n — m)) 

We a/m to so/ve ?/ze5e equations for the first m variables, y, which will then be 
written as functions, h (z) of the last n — m variables, z. 

Suppose g is continuously differentiable in a neighbourhood of x*, and that the 
m x m matrix: 

( §£(**) ••• fi>*)\ 

\|M ••• £(x*)/ 

formed by the first m columns of the total derivative of g at x* is non-singular. 
Then 3 neighbourhoods Y of y * and Z of z*, and a continuously differentiable 
function h: Z —> Y such that 

1. y* = h(z*), 

2. g (h (z) , z) = Vz G Z, and 

3. h'(z*) = -(D y gy 1 D z g. 

Proof The full proof of this theorem, like that of Brouwer's Fixed Point Theorem 
later, is beyond the scope of this course. However, part 3 follows easily from 
material in Section 2.5. The aim is to derive an expression for the total derivative 
h! (z*) in terms of the partial derivatives of g, using the Chain Rule. 
We know from part 2 that 

f{z)=g{h{z),z) = O m VzeZ. 

Thus 

/' (z) = mx(n _ m) Vz G Z, 

in particular at z*. But we know from (2.5.4) that 

/' (z) = D y gti (z) + D z g. 
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Hence 

D y gti (z) + D z g = O mx ( n -m) 

and, since the statement of the theorem requires that D y g is invertible, 

h'(z*) = -(D y g)- 1 D z g J 

as required. 

Q.E.D. 

To conclude this section, consider the following two examples: 

1. the equation g (x, y) = x 2 + y 2 — 1 = 0. 

Note that g' (x, y) = (2x 2y). 

We have h(y) = \j\ — y 2 or h(y) = — \j\ — y 2 , each of which describes a 
single-valued, differentiable function on (— 1, 1). At (x,y) = (0, 1), ^| = 
and h(y) is undefined (for y > 1) or multi-valued (for y < 1) in any 
neighbourhood of y = 1. 

2. the system of linear equations g (x) = Bx = 0, where B is an m x n 
matrix. 

We have g' (x) = B Vx so the implicit function theorem applies provided 
the equations are linearly independent. 

2.7 Directional Derivatives 

Definition 2.7.1 Let X be a vector space and x ^ x' G I. Then 

1. for A G 3? and particularly for A G [0, 1], Ax + (1 — A) x' is called a convex 
combination o/x and x'. 

2. L = {Ax + (1 — A) x' : A G 3?} is the line from x', where A = 0, to x, 
where A = 1, in X. 

3. The restriction of the function /: X — ► 3? to the line L is the function 

/| L :K^3?:A^/(Ax+(l-A)x / ). 

4. If f is a differentiable function, then the directional derivative of / at x' in 
the direction from x' to x is f\' L (0). 
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• We will endeavour, wherever possible, to stick to the convention that x' 
denotes the point at which the derivative is to be evaluated and x denotes 
the point in the direction of which it is measured. 1 

• Note that, by the Chain Rule, 

f\' L (A) = /' (Ax + (1 - A) x') (x - x') (2.7.1) 

and hence the directional derivative 

/|' L (0) = /'(x')(x-x'). (2.7.2) 

• The ith partial derivative of / at x is the directional derivative of / at x in 
the direction from x to x + e i5 where e^ is the ith standard basis vector. In 
other words, partial derivatives are a special case of directional derivatives 
or directional derivatives a generalisation of partial derivatives. 

• As an exercise, consider the interpretation of the directional derivatives at a 
point in terms of the rescaling of the parameterisation of the line L. 

• Note also that, returning to first principles, 

,./ ^ / ( x ' + A (x - x')) - / (x') 
f\' L (0) = lim L± x • ( 2 - 7 - 3 ) 

• Sometimes it is neater to write x — x' = h. Using the Chain Rule, it is 
easily shown that the second derivative of /|l is 

/|I(A) = h T f(x' + Ah)h 

and 

f\l(0) = h T f (x>. 

2.8 Taylor's Theorem: Deterministic Version 

This should be fleshed out following ?. 

Readers are presumed to be familiar with single variable versions of Taylor's The- 
orem. In particular recall both the second order exact and infinite versions. 
An interesting example is to approximate the discount factor using powers of the 
interest rate: 

_ L_ = \-i + i 2 -i^ + i 4 + ... (2.8.1) 

1 + i 



'There may be some lapses in this version. 
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We will also use two multivariate versions of Taylor's theorem which can be ob- 
tained by applying the univariate versions to the restriction to a line of a function 
of n variables. 

Theorem 2.8.1 (Taylor's Theorem) Let f : X — ► 5ft be twice differentiable, 
X C !ft n . Then for any x, x' G X, 3 A G (0,1) smc/z tfwtf 

/(x) = Z(x') + /' (x')(x - x') + i (x - x') T f (x' + A(x - x'))(x - x'). (2.8.2) 

Proof Let L be the line from x' to x. 

Then the univariate version tells us that there exists A G (0, l) 2 such that 

/U(l) = /U(0) + /|' L (0) + i/|£(A). (2.8.3) 

Making the appropriate substitutions gives the multivariate version in the theorem. 

Q.E.D. 

The (infinite) Taylor series expansion does not necessarily converge at all, or to 
f(x). Functions for which it does are called analytic. ? is an example of a function 
which is not analytic. 

2.9 The Fundamental Theorem of Calculus 

This theorem sets out the precise rules for cancelling integration and differentia- 
tion operations. 

Theorem 2.9.1 (Fundamental Theorem of Calculus) The integration and dif- 
ferentiation operators are inverses in the following senses: 

1. 

7 

b f'(x)dx = f(b)-f(a) 

This can be illustrated graphically using a picture illustrating the use of integration 
to compute the area under a curve. 



Should this not be the closed interval? 
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Chapter 3 

CONVEXITY AND 
OPTIMISATION 

3.1 Introduction 

[To be written.] 

3.2 Convexity and Concavity 

3.2.1 Definitions 

Definition 3.2.1 A subset X of a vector space is a convex set 

Vx,x'gX,Ag [0,1], Ax + (1 - A)x' G X. 

Theorem 3.2.1 A sum of convex sets, such as 

X + Y = {x + y:xeX,yeY}, 
is also a convex set. 

Proof The proof of this result is left as an exercise. 



Q.E.D. 



Definition 3.2.2 Let f : X — > Y where X is a convex subset of a real vector 
space and yci Then 
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1. f is a convex function 

Vx^x'eX,Ae (0,1) 

/ (Ax + (1 - A)x') < A/(x) + (1 - A)/(x'). (3.2.1) 

(This just says that a function of several variables is convex if its restriction 
to every line segment in its domain is a convex function of one variable in 
the familiar sense.) 

2. f is a concave function 

Vx^x'eX,Ae (0,1) 

/ (Ax + (1 - A)x') > A/(x) + (1 - A)/(x'). (3.2.2) 

3. f is qffine 

f is both convex and concave. l 

Note that the conditions (3.2.1) and (3.2.2) could also have been required to hold 
(equivalently) 

Vx, x' e x, A e [o, l] 

since they are satisfied as equalities V/ when x = x', when A = and when 

A = l. 

Note that / is convex <=^> —/is concave. 

Definition 3.2.3 Again let f : X — ► Y where X is a convex subset of a real vector 
space and 7C£ Then 

I. f is a strictly convex function 

Vx^x'e X,Xe (0,1) 

/ (Ax + (1 - A)x') < A/(x) + (1 - A)/(x'). 



'A linear function is an affine function which also satisfies /(0) = 0. 
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2. / is a strictly concave function 

Vx ^ x' G X, A G (0, 1) 

/ (Ax + (1 - A)x') > A/(x) + (1 - A)/(x'). 

Note that there is no longer any flexibility as regards allowing x = x' or A = or 
A = 1 in these definitions. 

3.2.2 Properties of concave functions 

Note the connection between convexity of a function of several variables and con- 
vexity of the restrictions of that function to any line in its domain: the former is 
convex if and only if all the latter are. 

Note that a function on a multidimensional vector space, X, is convex if and only 
if the restriction of the function to the line L is convex for every line L in X, and 
similarly for concave, strictly convex, and strictly concave functions. 
Since every convex function is the mirror image of a concave function, and vice 
versa, every result derived for one has an obvious corollary for the other. In gen- 
eral, we will consider only concave functions, and leave the derivation of the 
corollaries for convex functions as exercises. 
Let / : X — ► !R and g : X —> 5ft be concave functions. Then 

1 . If a, b > 0, then af + bg is concave. 

2. If a < 0, then af is convex. 

3. min{/, g} is concave 

The proofs of the above properties are left as exercises. 

Definition 3.2.4 Consider the real-valued function f:X^Y where FCK. 

1. The upper contour sets of f are the sets {x G X : /(x) > a} (a G 3ft). 

2. The level sets or indifference curves of f are the sets {x G X : /(x) = a} 
(a G 9ty. 

3. The lower contour sets of f are the sets {x G X : /(x) < a} (a G 5R). 
In Definition 3.2.4, X does not have to be a (real) vector space. 
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Theorem 3.2.2 The upper contour sets {x G X : /(x) > a} o/a concave 
function are convex. 



Proof This proof is probably in a problem set somewhere. 



Q.E.D. 



Consider as an aside the two-good consumer problem. Note in particular the im- 
plications of Theorem 3.2.2 for the shape of the indifference curves corresponding 
to a concave utility function. Concave u is a sufficient but not a necessary condi- 
tion for convex upper contour sets. 

3.2.3 Convexity and differentiability 

In this section, we show that there are a total of three ways of characterising 
concave functions, namely the definition above, a theorem in terms of the first 
derivative (Theorem 3.2.3) and a theorem in terms of the second derivative or 
Hessian (Theorem 3.2.4). 

Theorem 3.2.3 [Convexity criterion for differentiate functions.] Let f : X — ► 5ft 
be differentiable, X C ffl 1 an open, convex set. Then: 
f is (strictly) concave 



Vx ^ x' G X, 

/(x)<(<)/(x') + / / (xO(x-x'). (3.2.3) 

Theorem 3.2.3 says that a function is concave if and only if the tangent hyperplane 
at any point lies completely above the graph of the function, or that a function is 
concave if and only if for any two distinct points in the domain, the directional 
derivative at one point in the direction of the other exceeds the jump in the value 
of the function between the two points. (See Section 2.7 for the definition of a 
directional derivative.) 

Proof (See ?.) 

1. We first prove that the weak version of inequality 3.2.3 is necessary for 
concavity, and then that the strict version is necessary for strict concavity. 

Choose x, x' G X. 
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(a) Suppose that / is concave. 
Then, for A G (0, 1), 

/ (x' + A(x - x')) > /(x') + A (/(x) - /(x')) • (3.2.4) 

Subtract J(x') from both sides and divide by A: 
/ (x' + A(x - x')) - /(x') 



A 



> /(x) - /(x'). (3-2.5) 



Now consider the limits of both sides of this inequality as A — ► 0. 
The LHS tends to /' (x') (x — x') by definition of a directional deriva- 
tive (see (2.7.2) and (2.7.3) above). The RHS is independent of A and 
does not change. The result now follows easily for concave functions. 
However, 3.2.5 remains a weak inequality even if / is a strictly con- 
cave function. 

(b) Now suppose that / is strictly concave and x ^ x'. 

Since / is also concave, we can apply the result that we have just 
proved to x' and x" = | (x + x') to show that 

/'(x')(x"-x')>/(x")-/(x'). (3.2.6) 

Using the definition of strict concavity (or the strict version of inequal- 
ity (3.2.4)) gives: 

/(x") - /(x') > \ (/(x) - /(x')) • (3.2.7) 

Combining these two inequalities and multiplying across by 2 gives 
the desired result. 

2. Conversely, suppose that the derivative satisfies inequality (3.2.3). We will 
deal with concavity. To prove the theorem for strict concavity, just replace 
all the weak inequalities (>) with strict inequalities (>), as indicated. 

Set x' = Ax + (1 — A)x". Then, applying the hypothesis of the proof in turn 
to x and x' and to x" and x' yields: 

/(x) < /(x') + f(x')(x-x') (3.2.8) 

and 

/(x") < /(x') + /'(x')(x"-x') (3.2.9) 

A convex combination of (3.2.8) and (3.2.9) gives: 
A/(x) + (1 - A)/(x") 
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< /(x') + /'(x') (A ((x - x')) + (1 - A) ((x» - x'))) 

= /(x'), (3.2.10) 

since 

A ((x - x')) + (1 - A) ((x" - x')) = Ax + (1 - A)x" - x = O n . (3.2.1 1) 
(3.2. 10) is just the definition of concavity as required. 

Q.E.D. 

Theorem 3.2.4 [Concavity criterion for twice differentiable functions.] Let f : 
X —* 3ft be twice continuously differentiable (C 2 ), X C 3t™ open and convex. 
Then: 

1. f is concave 

Vx e X, ?/ze Hessian matrix /"( x ) w negative semidefinite. 

2. /"(x) negative definite Vx e X 

/ w strictly concave. 

The fact that the condition in the second part of this theorem is sufficient but not 
necessary for concavity inspires the search for a counter-example, in other words 
for a function which is strictly concave but has a second derivative which is only 
negative semi-definite and not strictly negative definite. The standard counter- 
example is given by f(x) = x 2n for any integer n > 1. 

Proof We first use Taylor's theorem to demonstrate the sufficiency of the condi- 
tion on the Hessian matrices. Then we use the Fundamental Theorem of Calculus 
(Theorem 2.9. 1) and a proof by contrapositive to demonstrate the necessity of this 
condition in the concave case for n = 1. Then we use this result and the Chain 
Rule to demonstrate necessity for n > 1. Finally, we show how these arguments 
can be modified to give an alternative proof of sufficiency for functions of one 
variable. 

1. Suppose first that /"(x) is negative semi-definite Vx e X. Recall Taylor's 
Theorem above (Theorem 2.8.1). 

It follows that /(x) < /(x') + /'(x')(x - x'). Theorem 3.2.3 shows that / 
is then concave. A similar proof will work for negative definite Hessian and 
strictly concave function. 
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2. To demonstrate necessity, we must consider separately first functions of a 
single variable and then functions of several variables. 

(a) First consider a function of a single variable. Instead of trying to show 
that concavity of / implies a negative semi-definite {i.e. non-positive) 
second derivative Va; G I, we will prove the contrapositive. In other 
words, we will show that if there is any point x* G X where the second 
derivative is positive, then / is locally strictly convex around x* and 
so cannot be concave. 

So suppose /"(#*) > 0. Then, since / is twice continuously differ- 
entiable, f"(x) > for all x in some neighbourhood of x*, say (a, b). 
Then /' is an increasing function on (a, b). Consider two points in 
(a, b), x < x" and let x' = Ax + (1 - X)x" G X, where A G (0, 1). 
Using the fundamental theorem of calculus, 

/(*') - f(x) = f f'(t)dt < f'(x')(x' - x) 

J X 

and 



f(x")-f(x') = / f'(t)dt<f'(x')(x"-x'). 

J X 1 

Rearranging each inequality gives: 

f(x) > f(x') + f'(x')(x-x') 
and 

f( x ") > f(x') + f'(x')(x" -x'), 

which are just the single variable versions of (3.2.8) and (3.2.9). As in 
the proof of Theorem 3.2.3, a convex combination of these inequalities 
reduces to 

/(*') < \f(x) + (1 - \)f(x"), 

and hence / is locally strictly convex on (a, b). 

(b) Now consider a function of several variables. Suppose that / is con- 
cave and fix x G X and heK™. (We use an x, x + h argument instead 
of an x, x' argument to tie in with the definition of a negative definite 
matrix.) Then, at least for sufficiently small A, g{\) = /(x + Ah) also 
defines a concave function (of one variable), namely the restriction of 
/ to the line segment from x in the direction from x to x + h. Thus, 
using the result we have just proven for functions of one variable, g 
has non-positive second derivative. But we know from p. 25 above 
that g"(0) = h T /"(x)h, so /"(x) is negative semi-definite. 
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3. For functions of one variable, the above arguments can give an alternative 
proof of sufficiency which does not require Taylor's Theorem. In fact, we 
have something like the following: 

/" (x) < on (a, 6) => / locally strictly concave on (a, b) 

f" (x) < on (a, b) => f locally concave on (a, b) 

f" (x) > on (a, b) =>■ / locally strictly convex on (a, b) 

f" (x) > on (a, b) =>■ f locally convex on (a, b) 

The same results which we have demonstrated for the interval (o, b) also 
hold for the entire domain X (which of course is also just an open interval, 
as it is an open convex subset of 5ft). 

Q.E.D. 

Theorem 3.2.5 A non-decreasing twice differentiable concave transformation of 
a twice differentiable concave function (of several variables) is also concave. 

Proof The details are left as an exercise. 

Q.E.D. 

Note finally the implied hierarchy among different classes of functions: 
negative definite Hessian c strictly concave C concave = negative semidefinite 
Hessian. 

As an exercise, draw a Venn diagram to illustrate these relationships (and add 
other classes of functions to it later on as they are introduced). 
The second order condition above is reminiscent of that for optimisation and sug- 
gests that concave or convex functions will prove useful in developing theories of 
optimising behaviour. In fact, there is a wider class of useful functions, leading us 
to now introduce further definitions. 

3.2.4 Variations on the convexity theme 

Let X C $ft n be a convex set and / : X — ► 5ft a real-valued function defined on X. 
In order (for reasons which shall become clear in due course) to maintain consis- 
tency with earlier notation, we adopt the convention when labelling vectors x and 

x'that/(x')</(x). 2 



There may again be some lapses in this version. 
Revised: December 2, 1998 



CHAPTER 3. CONVEXITY AND OPTIMISATION 35 

Definition 3.2.5 Let 

C(a) = {x e X : /(x) > a}. 
Then f : X —> 5ft is quasiconcave 

Va G 5ft, C(a) « a convex set. 

Theorem 3.2.6 The following statements are equivalent to the definition ofqua- 
siconcavity: 

1. Vx ^ x' e I such that /(x') < /(x) and\/\ e (0, 1), 

/ (Ax +(1-A)x') >/(*')■ 

2. Vx, x'gI.Ag [0, 1], / (Ax + (1 - A)x') > min{/(x), /(x')}. 

3. (7/Y is differentiable,) Vx,x' e X 5wc/z ?/za? /(x) - /(x') > 0, /'(x')(x - 
x') > 0. 

Proof 1. We begin by showing the equivalence between the definition and 
condition l. 3 

(a) First suppose that the upper contour sets are convex. Let x and x' e X 
and let a = min{/(x), /(x')}. Then x and x' are in C(a). By the 
hypothesis of convexity, for any A e (0, 1), Ax + (1 — A)x' G C(a). 
The desired result now follows. 

(b) Now suppose that condition 1 holds. To show that C(a) is a convex 
set, we just take x and x' 6 C(a) and investigate whether Ax + (1 — 
A)x' e C(a). But, by our previous result, 

/ (Ax + (1 - A)x') > min{/(x), /(x')} > a 

where the final inequality holds because x and x' are in C(a). 

2. It is almost trivial to show the equivalence of conditions 1 and 2. 

(a) In the case where /(x) > /(x') or /(x') = min{/(x), /(x')}, there 
is nothing to prove. Otherwise, we can just reverse the labels x and x'. 
The statement is true for x = x' or A = or A = 1 even if / is not 
quasiconcave. 



3 This proof may need to be rearranged to reflect the choice of a different equivalent character- 
isation to act as definition. 
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(b) The proof of the converse is even more straightforward and is left as 
an exercise. 

3. Proving that condition 3 is equivalent to quasiconcavity for a differentiable 
function (by proving that it is equivalent to conditions 1 and 2) is much the 
trickiest part of the proof. 

(a) Begin by supposing that / satisfies conditions 1 and 2. Proving that 
condition 3 is necessary for quasiconcavity is the easier part of the 
proof (and appears as an exercise on one of the problem sets). Pick 
any A G (0, 1) and, without loss of generality, x and x' such that 
/( x < 7( x )- By quasiconcavity, 

/(Ax+(l-A)x')>/(x'). (3.2.12) 

Consider 

f\ L (X) = f (Ax + (1 - A)x') = / (x' + A(x - x')) . (3.2.13) 

We want to show that the directional derivative 

/|' L (0) = /'(x')(x-x')>0. (3.2.14) 

But 

/W^fe^^. (3.2,5, 

Since the right hand side is non-negative for small positive values of A 
(A < 1), the derivative must be non-negative as required. 

(b) Now the difficult part — to prove that condition 3 is a sufficient con- 
dition for quasiconcavity. 

Suppose the derivative satisfies the hypothesis of the theorem, but / is 
not quasiconcave. In other words, 3x, x', A* such that 

/(A*x+ (1 - A*)x') < min{/(x),/(x')}, (3.2.16) 

where without loss of generality /(x') < /(x). The hypothesis of the 
theorem applied first to x and A*x + (1 — A*)x' and then to x' and 
A*x+(1 - A* )x' tells us that: 

/' (A*x + (1 - A*)x') (x - (A*x + (1 - A*)x')) > (3.2.17) 
/' (A*x + (1 - A*)x') (x - (A*x + (1 - A*)x)) > 0. (3.2.18) 
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Multiplying the first inequality by (1 — A*) and the second by A* yields 
a pair of inequalities which can only be satisfied simultaneously if: 

/' (A*x + (1 - A*)x') (x - x) = 0. (3.2.19) 

Inotherwords, f\' L {\*) = 0; we already know that /| L (A*) < /|l(0) <| 
/|l(1). We can apply the same argument to any point where the value 
of f\i is less than /|l(0) to show that the corresponding part of the 
graph of f\i has zero slope, or is flat. But this is incompatible either 
with continuity of f\t, or with the existence of a point where /|l(A*) 
is strictly less than /|l(0). So we have a contradiction as required. 

Q.E.D. 

In words, part 3 of Theorem 3.2.6 says that whenever a differentiable quasicon- 
cave function has a higher value at x than at x', or the same value at both points, 
then the directional derivative of / at x' in the direction of x is non-negative. It 
might help to think about this by considering n = 1 and separating out the cases 
x > x' and x < x'. 

Theorem 3.2.7 Let f : X — ► 5ft be quasiconcave and g : 3ft — ► 3ft be increasing. 
Then g o / is a quasiconcave function. 

Proof This follows easily from the previous result. The details are left as an 
exercise. 

Q.E.D. 

Note the implications of Theorem 3.2.7 for utility theory. In particular, if pref- 
erences can be represented by a quasiconcave utility function, then they can be 
represented by a quasiconcave utility function only. This point will be considered 
again in a later section of the course. 

Definition 3.2.6 / : X — ► 3ft is strictly quasiconcave 

Vx ^ x' G X such that /(x) > /(x') and VA e (0, 1), / (Ax + (1 - A)x') > 
Z(x'). 

Definition 3.2.7 / is (strictly) quasiconvex 
— f is (strictly) quasiconcave 4 



4 EC3080 ended here for Hilary Term 1998. 
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Definition 3.2.8 / is pseudoconcave 
f is differentiable and quasiconcave and 

/(x) - /(x') > =► /'(x') (x - x') > 0. 

Note that the last definition modifies slightly the condition in Theorem 3.2.6 which 
is equivalent to quasiconcavity for a differentiable function. 
Pseudoconcavity will crop up in the second order condition for equality con- 
strained optimisation. 

We conclude this section by looking at a couple of the functions of several vari- 
ables which will crop repeatedly in applications in economics later on. 
First, consider the interesting case of the affine function 

/: 3T -> 9fc x i-> M - p T x, 

where MeK and p e 5R n . This function is both concave and convex, but neither 
strictly concave nor strictly convex. Furthermore, 

/ (Ax+(l-A)x / ) = A/(x) + (1 - A)/(x') (3.2.20) 

> min{/(x),/(x')} (3.2.21) 

and 

(-/)(Ax+(l-A)x') = A(-/)(x) + (l-A)(-/)(x') (3.2.22) 

> min{(-/)(x),(-/)(x')}, (3.2.23) 

so / is both quasiconcave and quasiconcave, but not strictly so in either case. / 
is, however, pseudoconcave (and pseudoconvex) since 

/(x)>/(x') «=* p T x<p T x / (3.2.24) 

^^ P T (x - x') < (3.2.25) 

^^ -/'(x')(x-x') < (3.2.26) 

<=> / , (x / )(x-x')>0. (3.2.27) 

Finally, here are two graphs of Cobb-Douglas functions: 5 

5 To see them you will have to have copied two . WMF files retaining their uppercase filenames 
and put them in the appropriate directory, 

C: /TCD/teaching/WWW/MA3 81 /NOTES/! 
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0.5,, 0.5 



Graph of z = x y 



0.5„,1.5 



Graph of z = x y 

3.3 Unconstrained Optimisation 

Definition 3.3.1 Let ICF,/:I^| 

Then f has a (strict) global maximum at x* 

Also f has a (strict) local maximum at x* 
£e(x*),x^x*,/(x)<(<)/(x*). 

Similarly for minima. 

Theorem 3.3.1 A continuous real-valued function on a compact subset ofW 1 at- 
tains a global maximum and a global minimum. 



Vx G X,x ^ x*,/(x) < (< 
> 3e > such that Vx G 



Proof Not given here. See ?. 



Q.E.D. 



While this is a neat result for functions on compact domains, results in calculus 
are generally for functions on open domains, since the limit of the first difference 
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of the function at x, say, makes no sense if the function and the first difference are 
not defined in some open neighbourhood of x (some Bt(x)). 
The remainder of this section and the next two sections are each centred around 
three related theorems: 

1. a theorem giving necessary or first order conditions which must be satis- 
fied by the solution to an optimisation problem (Theorems 3.3.2, 3.4.1 and 
3.5.1); 

2. a theorem giving sufficient or second order conditions under which a solu- 
tion to the first order conditions satisfies the original optimisation problem 
(Theorems 3.3.3, 3.4.2 and 3.5.2); and 

3. a theorem giving conditions under which a known solution to an optimisa- 
tion problem is the unique solution (Theorems 3.3.4, 3.4.3 and 3.5.3). 

The results are generally presented for maximisation problems. However, any 
minimisation problem is easily turned into a maximisation problem by reversing 
the sign of the function to be minimised and maximising the function thus ob- 
tained. 

Throughout the present section, we deal with the unconstrained optimisation prob- 
lem 

max/(x) (3.3.1) 

where X C 3ft n and / : X — ► 5ft is a real-valued function of several variables, 
called the objective function of Problem (3.3.1.) 

(It is conventional to use the letter A both to parameterise convex combinations 
and as a Lagrange multiplier. To avoid confusion, in this section we switch to the 
letter a for the former usage.) 

Theorem 3.3.2 Necessary (first order) condition for unconstrained maxima and 

minima. 

Let X be open and f differentiable with a local maximum or minimum at x* e X. 

Then /'(x*) = 0, or f has a stationary point at x*. 

Proof Without loss of generality, assume that the function has a local maximum 
at x*. Then 3e > such that, whenever ||h|| < e, 

/(x* + h) - /(x*) < 0. 

It follows that, for < h < e, 

/(x* + he*) - /(x* 



h 
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(where e^ denotes the ith standard basis vector) and hence that 

df , ., fix* + he,) - f(x*) 

— - x* = lim ^ — - — ^ J -^-L < 0. (3.3.2) 



Similarly, for > h > — e, 

/(x* + teO - /(x*) 



and hence 



>0, 



|/ (x . ) = lim /(^±^WM> ^ ( 33.3) 



Combining (3.3.2) and (3.3.3) yields the desired result. 



Q.E.D. 



The first order conditions are only useful for identifying optima in the interior 
of the domain of the objective function: Theorem 3.3.2 applies only to functions 
whose domain X is open. Other methods must be used to check for possible cor- 
ner solutions or boundary solutions to optimisation problems where the objective 
function is defined on a domain that is not open. 

Theorem 3.3.3 Sufficient (second order) condition for unconstrained maxima ana\ 

minima. 

Let X C ffl 1 be open and let f : X — > 5ft be a twice continuously differentiable 

function with f'{x*) = and f"(x*) negative definite. 

Then f has a strict local maximum at x*. 

Similarly for positive definite Hessians and local minima. 

Proof Consider the second order Taylor expansion used previously in the proof 
of Theorem 3.2.4: for any x e X, 3s e (0, 1) such that 

/(x) = /(x*) + /'(x*)(* - x*) + i(x - x*) T /" (x* + a (x - x*))(x - X*). 

or, since the first derivative vanishes at x*, 

/(x) = /(x*) + i(x - x*) T /" (x* + s(x - x*))(x - X*). 

Since /" is continuous, /"(x* + s(x — x*)) will also be negative definite for x 

in some open neighbourhood of x*. Hence, for x in this neighbourhood, /(x) < 

/(x*) and / has a strict local maximum at x*. 
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Q.E.D. 

The weak form of this result does not hold. In other words, semi-definiteness 
of the Hessian matrix at x* is not sufficient to guarantee that / has any sort of 
maximum at x*. For example, if f(x) = x 3 , then the Hessian is negative semi- 
definite at x = but the function does not have a local maximum there (rather, it 
has a point of inflexion). 

Theorem 3.3.4 Uniqueness conditions for unconstrained maximisation. 
If 

1. x* solves Problem (3.3.1) and 

2. f is strictly quasiconcave (presupposing that X is a convex set), 
then x* is the unique (global) maximum. 

Proof Suppose not, in other words that 3x ^ x* such that / (x) = / (x*). 
Then, for any a e (0, 1), 

/ (ax + (1 - a) x*) > min {/ (x) , / (x*)} = / (x*) , 

so / does not have a maximum at either x or x* . 

This is a contradiction, so the maximum must be unique. 

Q.E.D. 

Theorem 3.3.5 Tempting, but not quite true, corollaries of Theorem 3.3.3 are: 

• Every stationary point of a twice continuously differentiable strictly con- 
cave function is a strict global maximum (and so there can be at most one 
stationary point). 

• Every stationary point of a twice continuously differentiable strictly convex 
function is a strict global minimum. 

Proof If the Hessian matrix is positive/negative definite everywhere, then the ar- 
gument in the proof of Theorem 3.3.3 can be applied for x e X and not just for 
x G Be (x*). If there are points at which the Hessian is merely semi-definite, then 
the proof breaks down. 

Q.E.D. 

Note that many strictly concave and strictly convex functions will have no station- 
ary points, for example 

f:$^$:x^e a 
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3.4 Equality Constrained Optimisation: 
The Lagrange Multiplier Theorems 

Throughout this section, we deal with the equality constrained optimisation prob- 
lem 

max xeX / (x) 
s.t. £(x) = m (3.4.1) 

where X C 5ft n , / : X — ► 5ft is a real-valued function of several variables, called 
the objective function of Problem (3.4.1) and g : X —> $ft m is a vector-valued 
function of several variables, called the constraint function of Problem (3.4.1); or, 
equivalently, gi\ X — ► 5ft are real- valued functions for j = 1, . . . , m. In other 
words, there are m scalar constraints represented by a single vector constraint: 

/V(x) \ /o\ 



v <? m (x) ; 



V o y 



We will introduce and motivate the Lagrange multiplier method which applies to 
such constrained optimisation problems with equality constraints. We will assume 
where appropriate that the objective function / and the m constraint functions 
g 1 , . . . , g m are all once or twice continuously differentiable. 
The entire discussion here is again presented in terms of maximisation, but can 
equally be presented in terms of minimisation by reversing the sign of the objec- 
tive function. Similarly, note that the signs of the constraint function(s) can be 
reversed without altering the underlying problem. We will see, however, that this 
also reverses the signs of the corresponding Lagrange multipliers. The signifi- 
cance of this effect will be seen from the formal results, which are presented here 
in terms of the usual three theorems. 

Before moving on to those formal results, we briefly review the methodology for 
solving constrained optimisation problems which should be familiar from intro- 
ductory and intermediate economic analysis courses. 

If x* is a solution to Problem (3.4.1), then there exist Lagrange multipliers, 6 A = 
(Ai, . . . , A m ), such that 

/' (x*) + AV (x*) = n . 

Thus, to find the constrained optimum, we proceed as if optimising the Lagrangean\ 

£(x,\) = f(x*) + \ T g(x*). 
Note that 



As usual, a row of numbers separated by commas is used as shorthand for a column vector. 

Revised: December 2, 1998 



3.4. EQUALITY CONSTRAINED OPTIMISATION: 
AA THE LAGRANGE MULTIPLIER THEOREMS 



1 . C = f whenever g = and 

2. g = where C is optimised. 

Roughly speaking, this is why the constrained optimum of / corresponds to the 

optimum of C. 

The Lagrange multiplier method involves the following four steps: 

1. introduce the m Lagrange multipliers, A = (Ai, . . . , A m ). 

2. Define the Lagrangean C: X x 5ft' m — > 5ft, where 

X x 5ft m = {(x, A) : x G X, A G 5ft m } , 

by 

£(x,A)^/(x) + A T £(x). 

3. Find the stationary points of the Lagrangean, i.e. set CJ (x, A) = 0. Since 
the Lagrangean is a function of n + m variables, this gives n + m first order 
conditions. The first n are 

/' (x) + AV (x) = 

or 

d f m da j 

The last m are just the original constraints, 

5 (x) = 

or 

9 J (x) = j' = l,...,m. 

4. Finally, the second order conditions must be checked. 

As an example, consider maximisation of a utility function representing Cobb- 
Douglas preferences subject to a budget constraint. 

The first n first order or Lagrangean conditions say that the total derivative (or 
gradient) of / at x is a linear combination of the total derivatives (or gradients) of 
the constraint functions at x. 
Consider a picture with n = 2 and m = 1. 

Since the directional derivative along a tangent to a level set or indifference curve 
is zero at the point of tangency, x, (the function is at a maximum or minimum 
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along the tangent) or /' (x) (x' — x) =0, the gradient vector, /' (x) , must be 

perpendicular to the direction of the tangent, x' — x. 

At the optimum, the level sets of / and g have a common tangent, so /' (x) and 

g' (x) are collinear, or /' (x) = —Xg' (x). It can also be seen with a little thought 

that for the solution to be a local constrained maximum, A must be positive if g is 

quasiconcave and negative if g is quasiconvex (in either case, the constraint curve 

is the boundary of a convex set). 

We now consider the equality constrained optimisation problem in more depth. 

Theorem 3.4.1 First order (necessary) conditions for optimisation with equality 

constraints. 

Consider problem (3.4.1) or the corresponding minimisation problem. 

If 

1. x* solves this problem (which implies that g (x*) = 0), 

2. f and g are continuously differentiable, and 

3. the m x n matrix 

( £(**) ••• £w\ 

</(x*)= : •-. : 

vff(x*) ... £(x*)y 

is of rank m (i.e. there are no redundant constraints, both in the sense that 
there are fewer constraints than variables and in the sense that the con- 
straints which are present are 'independent'), 

then 3X* e K m such that f (x*) + X* T g' (x*) = (i.e. in R n , f (x*) is in the 
m— dimensional subspace generated by the m vectors g 1 (x*) , . . . , g m (x*)j. 

Proof The idea is to solve g (x*) = for m variables as a function of the other 
n — m and to substitute the solution into the objective function to give an uncon- 
strained problem with n — m variables. 

For this proof, we need Theorem 2.6.1 on p. 2.6.1 above, the Implicit Function 
Theorem. Using this theorem, we must find the m weights Ai, . . . , X m to prove 
that /' (x*) is a linear combination of g 1 (x*) , . . . , g m ' (x*). 
Without loss of generality, we assume that the first m columns of g' (x*) are lin- 
early independent (if not, then we merely relabel the variables accordingly). 
Now we can partition the vector x* as (y*,z*) and, using the notation of the 
Implicit Function Theorem, find a neighbourhood Z of z* and a function h defined 
on Z such that 

g (h (z) , z) = Vz G Z 
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and also 

ti(z*) = -(D y g)- 1 D z g. 

Now define a new objective function F: Z —> 3ft by 

F(z) = f(h(z),z). 

Since x* solves the constrained problem max x6 x / (x) subject to g (x) = 0, it 

follows that z* solves the unconstrained problem max ze ^F (z). (This is easily 

shown using a proof by contradiction argument.) 

Hence, z* satisfies the first order conditions for unconstrained maximisation of F, 

namely 

F'(z*) = 0. 

Applying the Chain Rule in exactly the same way as in the proof of the Implicit 
Function Theorem yields an equation which can be written in shorthand as: 

D y fh'(z) + DJ = 0. 
Substituting for h! (z) gives: 

D y f(D y g)- 1 D z g = DJ. 
We can also partition /' (x) as 

D y f{D y g)- l D y g DJ 
Substituting for the second sub-matrix yields: 

f(x) = ( D y f (Dyg)- 1 D y g D y f(D y g)- l D z g 
= D y f(D y g)- l (D y g D z g 



-A V (x) 



where we define 



■D y f (Dyg)- 1 . 



Q.E.D. 



Theorem 3.4.2 Second order (sufficient or concavity) conditions for maximisa- 
tion with equality constraints. 
If 

1. f and g are differentiable, 
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2. /' (x*) + A* g' (x*) = (i.e. the first order conditions are satisfied at x*), 

3. A* > Ofor j = 1, . . . ,m, 

4. f is pseudoconcave, and 

5. gi is quasiconcave for j = 1, . . . , m, 

then x* solves the constrained maximisation problem. 

It should be clear that non-positive Lagrange multipliers and quasiconvex con- 
straint functions can take the place of non-negative Lagrange multipliers and qua- 
siconcave Lagrange multipliers to give an alternative set of second order condi- 
tions. 

Proof Suppose that the second order conditions are satisfied, but that x* is not a 

constrained maximum. We will derive a contradiction. 

Since x* is not a maximum, 3x ^ x* such that g (x) = but / (x) > / (x*). 

By pseudoconcavity, / (x) - / (x*) > implies that /' (x*) (x - x*) > 0. 

Since the constraints are satisfied at both x and x*, we have g (x*) = g (x) = 0. 

By quasiconcavity of the constraint functions (see Theorem 3.2.6), gi (x)— gi (x*) =| 

implies that g j ' (x*) (x - x*) > 0. 

By assumption, all the Lagrange multipliers are non-negative, so 

/' (x*) (x - x*) + A*V (x*) (x - x*) > 0. 

Rearranging yields: 

[/' (x*) + A*V (x*)] (x - x*) > 0. 

But the first order condition guarantees that the LHS of this inequality is zero (not 
positive), which is the required contradiction. 

Q.E.D. 

Theorem 3.4.3 Uniqueness condition for equality constrained maximisation. 
If 

1. x* is a solution, 

2. f is strictly quasiconcave, and 

3. gi is an affine function (i.e. both convex and concave) for j = 1, . . . , m, 
then x* is the unique (global) maximum. 
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Proof The uniqueness result is also proved by contradiction. Note that it does not 
require any differentiability assumption. 

• We first show that the feasible set is convex. 
Suppose x^x* are two distinct solutions. 

Consider the convex combination of these two solutions x a = ax+(l — a)x*.| 
Since each <? J is affine and gi (x*) = gi (x) = 0, we have 

g J (xa) = ag j (x) + (1 - a) g j (x*) = 0. 

In other words, xa also satisfies the constraints. 

• To complete the proof, we find the required contradiction: 

Since / is strictly quasiconcave and / (x*) = / (x), it must be the case that 

/(x«)>/(4 

Q.E.D. 

The construction of the obvious corollaries for minimisation problems is left as 

an exercise. 

We conclude this section with 

Theorem 3.4.4 (Envelope Theorem.) Consider the modified constrained optimi- 
sation problem: 

max/(x, a.) subject to g (x, a) = 0, (3.4.2) 

where x e 5ft n , ae R q , f:W l+q -► 5ft and g:W l+q -► 5ft m (i.e. as usual f is 

the real-valued objective function and g is a vector of m real-valued constraint 

functions, but either or both can depend on exogenous or control variables a as 

well as on the endogenous or choice variables x). 

Suppose that the standard conditions for application of the Lagrange multiplier 

theorems are satisfied. 

Let x* (a.) denote the optimal choice of x for given a (x*: $l q — ► 3? n is called the 

optimal response function) and let M (a.) denote the maximum value attainable 

by f for given ol (M: 3? 9 — ► 5ft is called the envelope function). 

Then the partial derivative of M with respect to a, is just the partial derivative 

of the relevant Lagrangean, f + X g, with respect to «j, evaluated at the optimal 

value ofx. The dependence of the vector of Lagrange multipliers, X, on the vector 

a. should be ignored in calculating the last-mentioned partial derivative. 

Proof The Envelope Theorem can be proved in the following steps: 
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1. Write down the identity relating the functions M, f and x*: 

M(a) = /(x*(a),a) 

2. Use (2.5.4) to derive an expression for the partial derivatives |^ of M in 
terms of the partial derivatives of / and x*: 

M< (a) = DJ (x* (a) , a) x*' (a) + D a f (x* (a) , a) . 

3. The first order (necessary) conditions for constrained optimisation say that 

DJ (x* (a) , a) = -\ T D x g (x* (a) , a) 
and allow us to eliminate the J^ terms from this expression. 

4. Apply (2.5.4) again to the identity g (x*, a) = m to obtain 

D x # (x* (a) , a) x*' (a) + £> a # (x* (a) , a) = mxg . 

Finally, use this result to eliminate the J^- terms from your new expression 

forf^. 

Combining all these results gives: 

M' (a) = D a f (x* (a) , a) + A t Dq, 5 (x* (a) , a) , 
which is the required result. 

Q.E.D. 

In applications in economics, the most frequently encountered applications will 
make sufficient assumptions to guarantee that 

1. x* satisfies the first order conditions with each A« > 0; 

2. the Hessian /" (x*) is a negative definite (n x n) matrix; and 

3 . g is a linear function, 

so that x* is the unique optimal solution to the equality constrained optimisation 
problem. 
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3.5 Inequality Constrained Optimisation: 
The Kuhn-Tucker Theorems 

Throughout this section, we deal with the inequality constrained optimisation 
problem 

max xeX / (x) 
s.t. e/ 4 (x)>0, i = l,2,...,m (3.5.1) 

where once again X C $t n , f : X —> 5ft is a real- valued function of several 
variables, called the objective function of Problem (3.5.1) and g : X — ► 3? m 
is a vector- valued function of several variables, called the constraint function of 
Problem (3.5.1). 

Before presenting the usual theorems formally, we need some graphical motiva- 
tion concerning the interpretation of Lagrange multipliers. 
Suppose that the constraint functions are given by 

g (x, at) = at — h (x) . (3.5.2) 

The Lagrangean is 

£(x,A) = /(x) + A T (a-/i(x)). (3.5.3) 

Thus, using the Envelope Theorem, it is easily seen that the rate of change of the 
envelope function M (a) with respect to the 'level' of the ith underlying con- 
straint function h' 1 is: 

dM dC 

oai oai 
Thus 

• when Aj = 0, the envelope function is at its maximum, or the objective 
function at its unconstrained maximum, and the constraint is not binding. 

• when \ < 0, the envelope function is decreasing 

• when Xi > 0, the envelope function is increasing 

Now consider how the nature of the inequality constraint change as ot{ increases 
(as illustrated, assuming / quasiconcave as usual and K l quasiconvex or g l quasi- 
concave, so that the relationship between «; and A; is negative) 

h l (x) < a % (3.5.5) 
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Type of constraint 


Derivative of 
objective fn. 


Lagrange 
Multiplier 


Constraint 
fn. 


Binding/active 
Non-binding/inactive 


/'<0 

f = o 


II IV 

o o 


9 = 
9>0 



Table 3.1: Sign conditions for inequality constrained optimisation 

(or g % (x, a.) > 0). For values of ctj such that Aj > 0, this constraint is strictly 

binding. For values of «j such that Aj = 0, this constraint is just binding. For 

values of aj such that \ < 0, this constraint is non-binding. 

Note that in this setup, 

d 2 M d\i 

— ^ = -^ < 0, (3.5.6) 

oaf oai 

so that the envelope function is strictly concave (up to the unconstrained optimum, 

and constant beyond it). 

Thus we will find that part of the necessity conditions below is that the Lagrange 

multipliers be non-negative. (For equality constrained optimisation, the signs 

were important only when dealing with second order conditions). 

Consider also at this stage the first order conditions for maximisation of a function 

of one variable subject to a non-negativity constraint: 



They can be expressed as: 



max/(i) s.t. x > 0. 
f'{x*) = ifx>0 



(3.5.7) 
(3.5.8) 



The various sign conditions which we have looked at are summarised in Table 3.1. 

Theorem 3.5.1 Necessary (first order) conditions for optimisation with inequal- 
ity constraints. 
If 



I. x* solves Problem (3.5.1), with 
and 



0,i = 1,2,. ..,6 



5*(x*) >0,i = b+l,...,m 

(in other words, the first b constraints are binding (active) at x* and the 
last n — b are non-binding (inactive) at x*, renumbering the constraints if 
necessary to achieve this), 
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2. f and g are continuously differentiate, and 

3. the b x n submatrix of g' (x*), 

« of full rank b (i.e. there are no redundant binding constraints, both in 
the sense that there are fewer binding constraints than variables and in the 
sense that the constraints which are binding are 'independent'), 

then 3X e 3ft m such that f (x*) + X T g' (x*) = 0, with \ > Ofor i = l,2,...,m 
andg l (x*) = OifXi > 0. 

Proof The proof is similar to that of Theorem 3.4.1 for the equality constrained 
case. It can be broken into seven steps. 

1. Suppose x* solves Problem (3.5.1). 

We begin by restricting attention to a neighbourhood B e (x*) throughout 
which the non-binding constraints remain non-binding, i.e. 

g i (x) > Vx G B e (x*) , i = b + 1, . . ., m. (3.5.9) 

Such a neighbourhood exists since the constraint functions are continuous. 
Since x* solves Problem (3.5.1) by assumption, it also solves the following 
problem: 

max xg £ £ ( x .) / (x) 
s.t. ^(x)>0, £ i = 1,2, ... ,6. (3.5.10) 

In other words, since the non-binding constraints are non-binding Vx e 
B e (x*) by construction, we can ignore them if we confine our search for 
a maximum to this neighbourhood. We will return to the non-binding con- 
straints in the very last step of this proof, but until then g will be taken to 
refer to the vector of b binding constraint functions only and A to the vector 
of b Kuhn-Tucker multipliers corresponding to these binding constraints. 

2. We now introduce slack variables s = (s l7 . . ., s fe ), one corresponding to 
each binding constraint, and consider the following equality constrained 
maximisation problem: 

max xgB £ (x*),sg» 6 + /( X ) 
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s.t. G(x,s)=0 6 (3.5.11) 

where G l (x, s) = g % (x) — s i5 i = 1, 2, . . . , b. 

Since x* solves Problem (3.5.10) and all b constraints in that problem are 
binding at x*, it can be seen that (x*,0&) solves this new problem. For 
consistency of notation, we define s* = b . 

3. We proceed with Problem (3.5.1 1) as in the Lagrange case. In other words, 
we use the Implicit Function Theorem to solve the system of b equations in 
n + b unknowns, 

G(x,s) = 6 , (3.5.12) 

for the first b variables in terms of the last n. To do this, we partition the 
vector of choice and slack variables three ways: 

(x,s) = (y,z,s) (3.5.13) 

where y G ?R. b and z G $t n ~ b , and correspondingly partition the matrix of 
partial derivatives evaluated at the optimum: 

G'(y*,z*,s*) - G'(x*,s*) (3.5.14) 

D y G D Z , S G ) (3.5.15) 

D y G D Z G D S G ) (3.5.16) 

D y g D z g -l b ) . (3.5.17) 



The rank condition allows us to apply the Implicit Function Theorem and 
to find a function h : 3? n — > ?R. b such that y = h (z, s) is a solution to 

G (y, z, s) = with 

ti (z*, s*) = - (D y gy l D^ S G. (3.5.18) 

(3.5.18) can in turn be partitioned to yield 

D z h = - (Dyg)- 1 D Z G = - (Dyg)- 1 D z g (3.5.19) 

and 

DJi = - (Dyg)- 1 D S G = (Dyg)- 1 I b = (Dyg)- 1 . (3.5.20) 

4. This solution can be substituted into the original objective function / to 
create a new objective function F defined by 

F(z,s) = f(h(z,s),z) (3.5.21) 
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and another new maximisation problem where there are only (implicit) non- 
negativity constraints: 

max zeBe(z « )jSe ^ F (z, s) (3.5.22) 

It should be clear that z*, b solves Problem (3.5.22). The first order con- 
ditions for Problem (3.5.22) are just that the partial derivatives of F with 
respect to the remaining n — b choice variables equal zero (according to 
the first order conditions for unconstrained optimisation), while the partial 
derivatives of F with respect to the b slack variables must be less than or 
equal to zero. 

5. The Kuhn-Tucker multipliers can now be found exactly as in the Lagrange 
case. We know that 

D z F = D y fD z h + DJI n _ b = n _ b . 

Substituting for D z h from (3.5.19) gives: 

D y f(D y g)- 1 D z g = DJ. 

We can also partition /' (x) as 

( D y f {D y gY l D y g D z f ) 
Substituting for the second sub-matrix yields: 

/'(x) = ( D y f (Dyg)- 1 D y g D y f(D y g)- 1 D z g 
= Dyf (Dyg)- 1 ( D y g D z g) 

= -aV(x) 

where we define the Kuhn-Tucker multipliers corresponding to the binding 
constraints, A , by 

A = -Dyf (Dyg)- 1 . (3.5.23) 

6. Next, we calculate the partial derivatives of F with respect to the slack 
variables and show that they can be less than or equal to zero if and only 
if the Kuhn-Tucker multipliers corresponding to the binding constraints are 
greater than or equal to zero. This can be seen by differentiating both sides 
of (3.5.21) with respect to s to obtain: 



D S F 


= D y fD s h + D z fO (n _ b)xb 


(3.5.24) 




= Dyf (Dyg)- 1 


(3.5.25) 




= -A, 


(3.5.26) 



where we have used (3.5.20) and (3.5.23). 
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7. Finally just set the Kuhn-Tucker multipliers corresponding to the non-binding| 
constraints equal to zero. 

Q.E.D. 

Theorem 3.5.2 Second order (sufficient or concavity) conditions for optimisation 

with inequality constraints. 

If 

1. f and g are differentiable, 

2. 3X e W l such that f (x*) + \ T g' (x*) = 0, with A; > Ofor i = 1, 2, . . . , m 
and g l (x*) = ifXi > (i.e. the first order conditions are satisfied at x*J, 

3. f is pseudoconcave, and 

4. g i (x*) = 0, i = 1, 2, . . . , b; g l (x*) > 0,i = b + 1, . . . , m and g j is qua- 
siconcave for j = 1, . . . , b (i.e. the binding constraint functions are quasi- 
concave), 

then x* solves the constrained maximisation problem. 

Proof The proof just requires the first order conditions to be reduced to 

f (x*) + 5>/(x*) = (3.5.27) 

from where it is virtually identical to that for the Lagrange case, and so it is left as 
an exercise. 

Q.E.D. 

Theorem 3.5.3 Uniqueness condition for inequality constrained optimisation. 
If 

1. x* is a solution, 

2. f is strictly quasiconcave, and 

3. gi is a quasiconcave function for j = 1, . . . , m, 
then x* is the unique (global) optimal solution. 
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Proof The proof is again similar to that for the Lagrange case and is left as an ex- 
ercise. The point to note this time is that the feasible set with equality constraints 
was convex if the constraint functions were affine, whereas the feasible set with 
inequality constraints is convex if the constraint functions are quasiconcave. This 
is because the feasible set (where all the inequality constraints are satisfied simul- 
taneously) is the intersection of m upper contour sets of quasiconcave functions, 

or the intersection of m convex sets. _ ^ _^ 

Q.E.D. 

The last important result on optimisation, the Theorem of the Maximum, is closely 
related to the Envelope Theorem. Before proceeding to the statement of the theo- 
rem, the reader may want to review Definition 2.3.12. 

Theorem 3.5.4 (Theorem of the maximum) Consider the modified inequality con\ 
strained optimisation problem: 

max/(x, ex) subject to g l (x, ex) > 0,i = 1, . . . ,m (3.5.28) 

where x G 3ft n , a€ 3ft 9 , f: 3ft n+9 -► 3ft and g: 3ft n+9 -► 3ft m 

Let x* (ex) denote the optimal choice of x for given ex (x*: 3ft 9 — > 3ft"j and let 

M (ex) denote the maximum value attainable by f for given ex (M: 3ft 9 — ► 3ft). 

// 

1. f is continuous; 

2. the range of f is closed and bounded; and 

3. the constraint set is a non-empty, compact-valued, continuous correspon- 
dence of ex, 

then 

1. M is a continuous (single-valued) function; and 

2. x* is an upper hemi-continuous correspondence, and hence is continuous if 
it is a continuous (single-valued) function. 

Proof The proof of this theorem is omitted, along with some of the more technical 

material on continuity of (multi-valued) correspondences. 

Q.rL.D. 

Theorem 3.5.4 will be used in consumer theory to prove such critical results as 
the continuity of demand functions derived from the maximisation of continuous 
utility functions. 

The following are two frequently encountered examples illustrating the use of the 
Kuhn-Tucker theorems in economics. 7 



7 The calculations should be left as exercises. 
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1. The canonical quadratic programming problem. 

Find the vector x G 3?" which maximises the value of the quadratic form 
x T Ax subject to the m linear inequality constraints g' T x > «j, where 
A G 3? nxn is negative definite and g* G ffl 1 for i = 1, . . . , m. 

The objective function can always be rewritten as a quadratic form in a 
symmetric (negative definite) matrixl, since as x T Ax is a scalar, 

x T Ax = fx T Ax) T (3.5.29) 



= x T A T x (3.5.30) 

= - (x T Ax + x T A T x) (3.5.31) 

= x T ^(A + A T ))x (3.5.32) 

and | (A + A T J is always symmetric. 

Let G be the m x n matrix whose ith row is g\ G must have full rank if 
we are to apply the Kuhn-Tucker conditions. 

The Lagrangean is: 

x t Ax + A t (Gx-q). (3.5.33) 

The first order conditions are: 

2x T A + A T G = n (3.5.34) 

or, transposing and multiplying across by |A _1 : 

x = — A _1 G T A. (3.5.35) 

2 

If the constraints are binding, then we will have: 

a = -^GA" 1 G T A. (3.5.36) 

Now we need the fact that G (and hence GA _1 G T ) has full rank to solve 
for the Lagrange multipliers A: 

A = -2(gA- 1 G t )" 1 q. (3.5.37) 

Now the sign conditions tell us that each component of A must be non- 
negative. An easy fix is to let the Kuhn-Tucker multipliers be defined by: 

A* = max|o m ,-2(GA- 1 G T )~ 1 a|, (3.5.38) 
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where the max operator denotes component-by-component maximisation. 

The effect of this is to knock out the non-binding constraints (those with 
negative Lagrange multipliers) from the original problem and the subse- 
quent analysis. 

We can now find the optimal x by substituting for A in (3.5.35) the value of 
A* from (3.5.38). 

In the case in which all the constraints are binding, the solution is: 

x = A^G 7 (GA^G 7 )" 1 a (3.5.39) 



Wga^gT 1 



and the envelope function is given by: 

x T Ax = a T fGA- 1 G T V 1 GA- 1 AA- 1 G T fGA- 1 G T r 1 a 

OL 

-a T X. (3.5.40) 



The applications of this problem will include ordinary least squares and 
generalised least squares regression and the mean-variance portfolio choice 
problem in finance. 

2. Maximising a Cobb-Douglas utility function subject to a budget constraint 
and non-negativity constraints. 

The applications of this problem will include choice under certainty, choice 
under uncertainty with log utility where the parameters are reinterpreted as 
probabilities, the extension to Stone-Geary preferences, and intertemporal 
choice with log utility, where the parameters are reinterpreted as time dis- 
count factors. 

Further exercises consider the duals of each of the forgoing problems, and it is to 
the question of duality that we will turn in the next section. 

3.6 Duality 

Let ICS B and let /, g : X \— > 5ft be, respectively, pseudoconcave and pseudo- 
convex functions. Consider the envelope functions defined by the dual families of 
problems: 

M (a) = max/(x)s.t. g (x.) < a (3.6.1) 
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and 

iV ((3) = min g (x) s.t. / (x) > (3. (3.6.2) 

Suppose that these problems have solutions, say x* (a) and x" 1 " (/?) respectively, 

and that the constraints bind at these points. 

The first order conditions for the two problems are respectively: 

/'(x)-A</(x) = (3.6.3) 

and 

</(x)+/zf (x) = 0, (3.6.4) 

where A and \x are the relevant Lagrange multipliers. 

Thus if x and A* ^ solve (3.6.3), then x and jjl* = — 1/A solve (3.6.4). However, 
for the x which solves the original problem (3.6.1) to also solve (3.6.2), it must 
also satisfy the constraint, or / (x) = (5. But we know that / (x) = M (a) . This 
allows us to conclude that: 

x*(a) = x f (M(a)). (3.6.5) 

Similarly, 

xt(/3) = x*(AT(/5)). (3.6.6) 

Combining these equations leads to the conclusion that 

a = N(M(a)) (3.6.7) 

and 

f3 = M(N((3)). (3.6.8) 

In other words, the envelope functions for the two dual problems are inverse func- 
tions (over any range where the Lagrange multipliers are non-zero, i.e. where the 
constraints are binding). Thus, either a or or indeed A or ji can be used to 
parameterise either family of problems. 

We will see many examples of these principles in the applications in the next part 
of the book. In particular, duality will be covered further in the context of its 
applications to consumer theory in Section 4.6. 
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Chapter 4 

CHOICE UNDER CERTAINTY 



4.1 Introduction 

[To be written.] 

4.2 Definitions 

There are two possible types of economy which we could analyse: 

• a pure exchange economy, in which households are endowed directly with 
goods, but there are no firms, there is no production, and economic activity 
consists solely of pure exchanges of an initial aggregate endowment; and 

• a production economy, in which households are further indirectly endowed 
with, and can trade, shares in the profit or loss of firms, which can use part 
of the initial aggregate endowment as inputs to production processes whose 
outputs are also available for trade and consumption. 

Economies of these types comprise: 

• H households or consumers or agents or (later) investors or, merely, indi- 
viduals, indexed by the subscript 1 h; 

• N goods or commodities, indexed by the superscript n; and 

• (in the case of a production economy only) F firms, indexed by /. 



'Notation in what follows is probably far from consistent in regard to superscripts and sub- 
scripts and in regard to ith or nth good and needs fixing. 
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This chapter concentrates on the theory of optimal consumer choice and of equi- 
librium in a pure exchange economy. The theory of optimal production decisions 
and of equilibrium in an economy with production is mathematically similar. 
Goods can be distinguished from each other in many ways: 

• obviously, by intrinsic physical characteristics, e.g. apples or oranges. 

• by the time at which they are consumed, e.g. an Easter egg delivered before 
Easter Sunday or an Easter egg delivered after Easter Sunday. (While all 
trading takes place simultaneously in the model, consumption can be spread 
out over many periods.) 

• by the state of the world in which they are consumed, e.g. the service pro- 
vided by an umbrella on a wet day or the service which it provides on a dry 
day. (Typically, the state of the world can be any point u in a sample space 
Q.) 

The important characteristics of household h are that it is faced with the choice of 
a consumption vector or consumption plan or consumption bundle, x^ = (x\, . . . , x% 

from a (closed, convex) consumption set, Xh- Typically, Xh = 5R+ . More gener- 
ally, consumer h's consumption set might require a certain subsistence consump- 
tion of some commodities, such as water, and rule out points of 3?+ not meeting 
this requirement. The household's endowments are denoted e h G 3ft+ and can be 
traded. 
In a production economy, the shareholdings of households in firms are denoted 

A consumer's net demand is denoted z h = x h — e h E ?R. N . 

Each consumer is assumed to have a preference relation or (weak) preference 

ordering which is a binary relation on the consumption set X h (?, Chapter 7). 

Since each household will have different preferences, we should really denote 

household /i's preference relation y^, but the subscript will be omitted for the 

time being while we consider a single household. Similarly, we will assume for 

the time being that each household chooses from the same consumption set, X, 

although this is not essential. 

Recall (see ?) that a binary relation R on X is just a subset R of X x X or a 

collection of pairs (x, y) where x e X and y e X. 

If (x, y) e R, we usually just write xRy. 

Thus x y y means that either x is preferred to y or the consumer is indifferent 

between the two (i.e. that x is at least as good as y). 

The following properties of a general relation R on a general set X are often of 

interest: 
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1. A relation R is reflexive 

xRx \/x e X 

2. A relation R is symmetric 

xRy =>■ yRx 

3. A relation R is transitive 

xRy,yRz =>■ xRz 

4. A relation R is complete 

Vrc, y E X either xRy or yi?x (or both) (in other words a complete relation 
orders the whole set) 

An indifference relation, ~, and a strict preference relation, y, can be derived 
from every preference relation: 

1 . x y y means x y y but not y y x 

2. x ~ y means x y y and y y x. 

The utility function u : X — ► 5ft represents the preference relation ^ if 

u(x) > w(y) <^=^ x y y. 

If /: 5ft — ► 5ft is a monotonic increasing function and u represents the preference 
relation y, then / o u also represents y, since 

f (u(x)) > f (u(y)) <=> u(x)>u(y) <=> x^y. 

If X is a countable set, then there exists a utility function representing any pref- 
erence relation on X. To prove this, just write out the consumption plans in X in 
order of preference, and assign numbers to them, assigning the same number to 
any two or more consumption plans between which the consumer is indifferent. 
If X is an uncountable set, then there may not exist a utility function representing 
every preference relation on X. 
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4.3 Axioms 

We now consider six axioms which it are frequently assumed to be satisfied by 
preference relations when considering consumer choice under certainty. (Note 
that symmetry would not be a very sensible axiom!) Section 5.5.1 will consider 
further axioms that are often added to simplify the analysis of consumer choice 
under uncertainty. After the definition of each axiom, we will give a brief ratio- 
nale for its use. 

Axiom 1 (Completeness) A (weak) preference relation is complete. 

Completeness means that the consumer is never agnostic. 

Axiom 2 (Reflexivity) A (weak) preference relation is reflexive. 

Reflexivity means that each bundle is at least as good as itself. 

Axiom 3 (Transitivity) A (weak) preference relation is transitive. 

Transitivity means that preferences are rational and consistent. 

Axiom 4 (Continuity) The preference relation y is continuous i.e. for all con- 
sumption plans y G X the sets B y = {x e X : x y y} and W y = {x G X : y y 
x} are closed sets. 

Consider the picture when N = 2: 

We will see shortly that B y , the set of consumption plans which are better than or 
as good as y, and W y , the set of consumption plans which are worse than or as 
good as y, are just the upper contour sets and lower contour sets respectively of 
utility functions, if such exist. 
E.g. consider lexicographic preferences: 

Lexicographic preferences violate the continuity axiom. A consumer with such 
preference prefers more of commodity 1 regardless of the quantities of other com- 
modities, more of commodity 2 if faced with a choice between two consumption 
plans having the same amount of commodity 1, and so on. 
In the picture, the consumption plan y lies in the lower contour set W x * but Be (y) 
never lies completely in W x * for any e. Thus, lower contour sets are not open, and 
upper contour sets are not closed. 

Theorems on the existence of continuous utility functions have been proven by 
Gerard Debreu, Nobel laureate, whose proof used Axioms 1-4 only (see ? or ?) 
and by Hal Varian, whose proof was simpler by virtue of adding an additional 
axiom (see ?). 
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Theorem 4.3.1 (Debreu) If X is a closed and convex set and y is a complete, 
reflexive, transitive and continuous preference relation on X, then 3 a continuous 
utility function u: X — ► 3ft representing y. 

Proof For the proof of this theorem, see ?. O F D 

Axiom 5 (Greed) Greed is incorporated into consumer behaviour by assuming 
either 

1. Strong mono tonicity: 

If X = 3ft+, then y is said to be strongly monotonic iff whenever x n > y n Vn 
but x ^ y, x y y. /x = (x±, . . . , Xn) &c.]; or 

2. Local non- satiation: 
VxG^,e>0, 3x' G Be (x) s.t. x' y x 

The strong monotonicity axiom is a much stronger restriction on preferences than 

local non-satiation; however, it greatly simplifies the proof of existence of utility 

functions. 

We will prove the existence, but not the continuity, part of the following weaker 

theorem. 

Theorem 4.3.2 (Varian) If X = 3ft+ and y is a complete, reflexive, transitive, 
continuous and strongly monotonic preference relation on X, then 3 a continuous 
utility function u: X — ► 3? representing y. 

Proof (of existence only (?, p. 97)) 

Pick a benchmark consumption plan, e.g. 1 = (1,1,.. .,1). 

The idea is that the utility of x is the multiple of the benchmark consumption plan 

to which x is indifferent. 

By strong monotonicity, the sets {t e 3? : £1 y x} and {t E 3? : x y tl} are both 

non-empty. 

By continuity of preferences, both are closed (intersection of a ray through the 

origin and a closed set); and by completeness, they cover 3?. 

By connectedness of 3ft, they intersect in at least one point, w(x) say, and x ~ 

«(x)l. 

Now 

x h y <^> «(x)l y «(y)l 
^> u(x) > «(y) 

where the first equivalence follows from transitivity of preferences and the second 
from strong monotonicity. 

The assumption that preferences are reflexive is not used in establishing existence 
of the utility function, so it can be inferred that it is required to establish continuity. 
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Q.E.D. 

The rule which the consumer will follow is to choose the most preferred bundle 
from the set of affordable alternatives {budget set), in other words the bundle at 
which the utility function is maximised subject to the budget constraint, if one 
exists. 

We know that an optimal choice will exist if the utility function is continuous and 
the budget set is closed and bounded. 

If the utility function is differentiable, we can go further and use calculus to find 
the maximum. 

So we usually assume differentiability. 

If u is a concave utility function, then /on, which also represents the same prefer- 
ences, is not necessarily a concave function (unless / itself is a convex function). 
In other words, concavity of a utility function is a property of the particular repre- 
sentation and not of the underlying preferences. Notwithstanding this, convexity 
of preferences is important, as indicated by the use of one or other of the follow- 
ing axioms, each of which relates to the preference relation itself and not to the 
particular utility function chosen to represent it. 

Axiom 6 (Convexity) There are two versions of this axiom: 

1. Convexity: 

The preference relation y is convex <=^> 

x h y =► Ax + (1 - A) y y y. 

2. Strict convexity: 

The preference relation y is strictly convex -<=>• 

x h y =*> Ax + (1 - A) y y y. 

The difference between the two versions of the convexity axiom basically amounts 
to ruling out linear segments in indifference curves in the strict case. 

Theorem 4.3.3 The preference relation y is (strictly) convex if and only if every 
utility function representing y is a (strictly) quasiconcave function. 

Proof In either case, both statements are equivalent to saying that 

«(x) > «(y) => u (Ax + (1 - A) y) > (>)«(y). (4.3.1) 
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Q.E.D. 

Theorem 4.3.4 All upper contour sets of a utility function representing convex 
preferences are convex sets (i.e. indifference curves are convex to the origin for 
convex preferences). 

It can be seen that convexity of preferences is a generalisation of the two-good 
assumption of a diminishing marginal rate of substitution. 

Axiom 7 (Inada Condition) This axiom seems to be relegated to some other cat- 
egory in many existing texts and its attribution to Inada is also difficult to trace. It 
is not easy to see how to express it in terms of the underlying preference relation 
so perhaps it cannot be elevated to the status of an axiom. Should it be expressed 

as: 

du 

hm — — = oo 

or as 

hm *L = 0? 

Xi^oo Q x . 

The Inada condition will be required to rule out corner solutions in the consumer's 
problem. Intuitively, it just says that indifference curves may be asymptotic to the 
axes but never reach them. 



4.4 Optimal Response Functions: 

Marshallian and Hicksian Demand 

4.4.1 The consumer's problem 

A consumer with consumption set X h , endowment vector e h G X h , sharehold- 
ings c h G 3? F and preference ordering >z h represented by utility function u h who 
desires to trade his endowment at prices p G !R+ faces an inequality constrained 
optimisation problem: 

maxw h (x) s.t. p T x < p T e h + c^II (p) = M h (4.4.1) 

where II (p) is the vector of the F firms' maximised profits when prices are p. 
Constraining x to lie in the consumption set normally just means imposing non- 
negativity constraints on the problem. From a mathematical point of view, the 
source of income is irrelevant, and in particular the distinction between pure ex- 
change and production economy is irrelevant. Thus, income can be represented 
by M in either case. 
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Since the constraint functions are linear in the choice variables x, the Kuhn-Tucker 
theorem on second order conditions (Theorem 3.5.2) can be applied, provided that 
the utility function Uh is pseudo-concave. In this case, the first order conditions 
identify a maximum. 

The Lagrangian, using multipliers A for the budget constraint and /i. E ?R. N for the 
non-negativity constraints, is 

u h (x) + A (M - p T x) + /z T x. (4.4.2) 

The first order conditions are given by the (iV-dimensional) vector equation: 

u' h (x) + A (-p) + n = N (4.4.3) 

and the sign condition A > with A > if the budget constraint is binding. 
We also have /i. = 0^ unless one of the non-negativity constraints is binding: 
Axiom 7 would rule out this possibility. 

Now for each p e ^-+ + (ruling out bads, or goods with negative prices, and even 
(see below) free goods), e^ G Xh, and c/ t e $t F , or for each p, Mh combination, 
there is a corresponding solution to the consumer's utility maximisation problem, 
denoted x h (p, e^, c h ) or x h (p, M h ). The function (correspondence) x h is often 
called a Marshallian demand function (correspondence). 

If the utility function Uh is also strictly quasiconcave (i.e. preferences are strictly 
convex), then the conditions of the Kuhn-Tucker theorem on uniqueness (Theo- 
rem 3.5.3) are satisfied. In this case, the consumer's problem has a unique solu- 
tion for given prices and income, so that the optimal response correspondence is a 
single- valued demand function. On the other hand, the weak form of the convexity 
axiom would permit a multi- valued demand correspondence. 



4.4.2 The No Arbitrage Principle 

Definition 4.4.1 An arbitrage opportunity means the opportunity to acquire a 
consumption vector or its constituents, directly or indirectly, at one price, and 
to sell the same consumption vector or its constituents, directly or indirectly, at a 
higher price. 

Theorem 4.4.1 (The No Arbitrage Principle) Arbitrage opportunities do not ex\ 
ist in equilibrium in an economy in which at least one agent has preferences which 
exhibit local non- satiation. 2 



2 The No Arbitrage Principle is also known as the No Free Lunch Principle, or the Law of One 
Price 
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Proof If preferences exhibit local non-satiation, then Marshallian demand is not 

well-defined if the price vector permits arbitrage opportunities. 

If the no arbitrage principle doesn't hold, then any individual can increase wealth 

without bound by exploiting the available arbitrage opportunity on an infinite 

scale. Thus there is no longer a budget constraint and, since local non-satiation 

rules out bliss points, utility too can be increased without bound. 

When we come to consider equilibrium, we will see that if even one individual has 

preferences exhibiting local non-satiation, then equilibrium prices can not permit 

arbitrage opportunities. 

Q.E.D. 

Examples are usually in the financial markets, in a multi-period context, e.g. cov- 
ered interest parity, term structure of interest rates, &c. The most powerful appli- 
cation is in the derivation of option-pricing formulae, since options can be shown 
to be identical to various synthetic portfolios made up of the underlying security 
and the riskfree security. 

The simple rule for figuring out how to exploit arbitrage opportunities is 'buy 
low, sell high'. With interest rates and currencies, for example, this may be a 
non-trivial calculation. 

Exercise: If the interest rate for one-year deposits or loans is r\ per annum com- 
pounded annually, the interest rate for two-year deposits or loans is r 2 per annum 
compounded annually and the forward interest rate for one year deposits or loans 
beginning in one year's time is f 12 per annum compounded annually, calculate the 
relationship that must hold between these three rates if there are to be no arbitrage 
opportunities. 
Solution: (1 + n) (1 + / 12 ) = (1 + r 2 ) 2 . 

4.4.3 Other Properties of Marshallian demand 

Other noteworthy properties of Marshallian demand include the following: 

1 . If preferences exhibit local non-satiation, then the budget constraint is bind- 
ing. 

This is because no consumption vector in the interior of the budget set can 
maximise utility as some nearby consumption vector will always be both 
preferred and affordable. At the optimum, on the budget hyperplane, the 
nearby consumption vector which is preferred will not be affordable. 

2. If p includes a zero price (p n = for some n), then x/j (p, Mh) may not be 
well defined. 
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This is because, at least in the case of strongly monotonic preferences, the 
consumer will seek to acquire and consume an infinite amount of the free 
good, thereby increasing utility without bound. For this reason, it is neater 
to define Marshallian demand only on the open non-negative orthant in $l N , 
namely 3ft+ + . 

3. The demand x h (p, M h ) is independent of the representation u h of the un- 
derlying preference relation y h which is used in the statement of the con- 
sumer's problem. 

4. x/j (p, M h ) is homogenous of degree in p, M h . 

In other words, if all prices and income are multiplied by a > 0, then 
demand does not change: 

x h (ap, aM h ) = x h (p, M h ) . {AAA) 

5. Demand functions are continuous. 

This follows from the theorem of the maximum (Theorem 3.5.4). It fol- 
lows that small changes in prices or income will lead to small changes in 
quantities demanded. 

4.4.4 The dual problem 

Consider also the (dual) expenditure minimisation problem: 

minp T x s.t. Uh (x) > u. (4.4.5) 

In other words, what happens if expenditure is minimised subject to a certain level 
of utility, u, being attained? 

The solution (optimal response function) is called the Hicksian or compensated 
demand function (or correspondence) and is usually denoted h/j (p, u). 
There should really be a more general discussion of duality, based on the mean- 
variance problem as well as the utility maximisation/expenditure minimisation 
problems. 

If the local non-satiation axiom holds, then the constraints are binding in both 
the utility-maximisation and expenditure-minimisation problems, and we have a 
number of duality relations. In particular, there will be a one-to-one correspon- 
dence between income M and utility u for a given price vector p. The expenditure 
function and the indirect utility function will then act as a pair of inverse envelope 
functions mapping utility levels to income levels and vice versa respectively. 
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A full page table setting out exactly the parallels between the two problems is 
called for here. 

The following duality relations (or fundamental identities as ? calls them) will 
prove extremely useful later on: 



e(p,v{p,M)) 


= M 


(4.4.6) 


v(p,e(p,u)) 


= u 


(4.4.7) 


x(p,M) 


= h(p,v(p,M)) 


(4.4.8) 


h(p,w) 


= x(p,e(p,«)) 


(4.4.9) 



These are just equations (3.6.5)-(3.6.8) adapted to the notation of the consumer's 
problem. 

4.4.5 Properties of Hicksian demands 

As in the Marshallian approach, if preferences are strictly convex, then any solu- 
tion to the expenditure minimisation problem is unique and the Hicksian demands 
are well-defined single valued functions. 

It's worth going back to the uniqueness proof with this added interpretation. If two 
different consumption vectors minimise expenditure, then they both cost the same 
amount, and any convex combination of the two also costs the same amount. But 
by strict convexity, a convex combination yields higher utility, and nearby there 
must, by continuity, be a cheaper consumption vector still yielding utility u. 
If preferences are not strictly convex, then Hicksian demands may be correspon- 
dences rather than functions. 
Hicksian demands are homogenous of degree in prices: 

h h (ap, u) = h h (p, u) . {A A. 10) 

4.5 Envelope Functions: 

Indirect Utility and Expenditure 

Now consider the envelope functions corresponding to the two approaches: 

1. The indirect utility function: 

v h (p,M)=u h (yi h (p,M)) (4.5.1) 

2. The expenditure function: 

e h (p,u) = p T h h (p,u). (4.5.2) 

Sometimes we meet two other related functions: 
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3. The money metric utility function 

m h {p,x) = e h (p,u h (x)) (4.5.3) 

is the (least) cost at prices p of being as well off as with the consumption 
vector x. 

4. The money metric indirect utility function 

Mp;q,M) = e h (p,Mq,M)) (4-5.4) 

is the (least) cost at prices p of being as well off as if prices were q and 
income was M. 

The following are interesting properties of the indirect utility function: 

1. By the Theorem of the Maximum, the indirect utility function is continuous 
for positive prices and income. 

2. The indirect utility function is non-increasing in p and non-decreasing in 

M. 

3. The indirect utility function is quasi-convex in prices. 

To see this, let B (p) denote the budget set when prices are p and let pA = 

Ap + (1-A)p'. 

Then B (pA) C (B (p) U B (p')). 

Suppose this was not the case, i.e. for some x, pA T x < M but p T x > 
M and p /T x > M. Then taking a convex combination of the last two 
inequalities yields 

Ap T x + (1 - A) p /T x > M, 

which contradicts the first inequality. 

It follows that the maximum value of u h (x) on the subset B (pA) is less 
than or equal to its maximum value on the superset B (p) U B (p'). 

In terms of the indirect utility function, this says that 

v h (pA, M) < max {v h (p, M) , v h (p\ M)} , 
or that Vh is quasiconvex. 

4. Vh (p, M) is homogenous of degree zero in p, M, or 

v h (Xp,XM) = v h (p,M). 
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The following are interesting properties of the expenditure function: 

1. The expenditure function is continuous. 

2. The expenditure function itself is non-decreasing in prices, since raising the 
price of one good while holding the prices of all other goods constant can 
not reduce the minimum cost of attaining a fixed utility level. 

3. The expenditure function is concave in prices. 

To see this, we just fix two price vectors p and p' and consider the value of 
the expenditure function at the convex combination pA = Ap + (1 — A) p'. 

e(pA,w) = (pA) T h(pA,IZ) (4.5.5) 

= Ap T h(pA,fZ) + (l- A)(p) T h(pA,M) (4.5.6) 

> Ap T h(p,fZ) + (l-A)(p , ) T h(p , ,fZ) (4.5.7) 

= Ae(p,«) + (1-A)e(p',«), (4.5.8) 

where the inequality follows because the cost of a suboptimal bundle for 
the given prices must be greater than the cost of the optimal (expenditure- 
minimising) consumption vector for those prices. 

4. The expenditure function is homogenous of degree 1 in prices: 

e h (ap, u) = ae h (p, u) . (4.5.9) 

4.6 Further Results in Demand Theory 

In this section, we present four important theorems on demand functions and the 
corresponding envelope functions. Shephard's Lemma will allow us to recover 
Hicksian demands from the expenditure function. Similarly, Roy's Identity will 
allow us to recover Marshallian demands from the indirect utility function. The 
Slutsky symmetry condition and the Slutsky equation provide further insights into 
the properties of consumer demand. 

Theorem 4.6.1 (Shephard's Lemma.) 

<)<h (p,u) = ^-(p T x + AK(x)-^)) (4.6.1) 



dp n dp n 

= x n (4.6.2) 

which, when evaluated at the optimum, is just h.% (p, u). 

In other words, the partial derivatives of the expenditure function with respect to 

prices are the corresponding Hicksian demand functions. 
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Proof By differentiating the expenditure function with respect to the price of 
good n and applying the envelope theorem (Theorem 3.4.4), we obtain Shephard's 
Lemma: 

(To apply the envelope theorem, we should be dealing with an equality constrained 
optimisation problem; however, if we assume local non- satiation, we know that 
the budget constraint or utility constraint will always be binding, and so the in- 
equality constrained expenditure minimisation problem is essentially and equality 
constrained problem.) 

Q.E.D. 

Theorem 4.6.2 (Roy's Identity.) Marshallian demands may be recovered from 
the indirect utility function using: 



"° (p,M) 

m (p» M ) ' 



* n (P,M) = -^^. (4.6.3) 



Proof For Roy's Identity, see ?. 

It is obtained by differentiating equation (4.4.7) with respect to p n , using the Chain 

Rule: 

v(p,e(p,u)) =u 

implies that 

dv dv dc 

— (p, e (p, «)) + — (p, e (p, «)) — (p, u) = (4.6.4) 

and using Shephard's Lemma gives: 

^- n (p, e (p, «)) + ^ (p, e (p, «)) h" (p, u) = (4.6.5) 

Hence 

|^(p,e(p,w)) 
h"(p^) = - t - (4-6-6) 

and expressing this last equation in terms of the relevant level of income M rather 
than the corresponding value of utility u: 



% (P, M) 

OM (P» M ) ' 



x"( P ,M) = - t n ;r;;; . (4.6.7) 



Q.E.D. 
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Theorem 4.6.3 (Slutsky symmetry condition.) All cross-price substitution effects^ 
are symmetric: 

dhl _ dh% 

Qpin Qpn 

Proof From Shephard's Lemma, we can easily derive the Slutsky symmetry condi- 
tions, assuming that the expenditure function is twice continuously differentiable, 
and hence that 

^ ^ (4.6.9) 



Opm Qpn dp n dp Tl 

_ _ h 

h dp m h dp n 



Since hY 1 = |% and hJl = |%, and the result follows. 



Q.E.D. 

The next result doesn't really have a special name of its own. 

Theorem 4.6.4 Since the expenditure function is concave in prices {seep. 75), the 
corresponding Hessian matrix is negative semi-definite. In particular, its diagonal 
entries are non-positive, or 

^<0, n = l,...,N. (4.6.10) 

Using Shephard's Lemma, it follows that 

dh n h 



Qpn 



< 0, n=l,..., AT. (4.6.11) 



In other words, Hicksian demand functions, unlike Marshallian demand functions, 
are uniformly decreasing in own price. Another way of saying this is that own 
price substitution effects are always negative. 

Theorem 4.6.5 (Slutsky equation.) The total effect of a price change on (Mar- 
shallian) demand can be decomposed as follows into a substitution effect and an 
income effect: 

Q r nn f)h m Q T ra 

W^ M) = W iP ^~^^ M)hn{P ^ h (4 - 6 - '^ 

where u = V (p, M). 

Before proving this, let's consider the signs of the various terms in the Slutsky 
equation and look at what it means in a two-good example. 
By Theorem 4.6.4, we know that own price substitution effects are always non- 
positive. 
[This is still on a handwritten sheet.] 
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Proof Differentiating both sides of the /th component of (4.4.9) with respect to 
p n , using the Chain Rule, will yield the so-called Slutsky equation which decom- 
poses the total effect on demand of a price change into an income effect and a 
substitution effect. 
Differentiating the RHS of (4.4.9) with respect to p n yields: 

(p,e(p,«)) H (p,e(p,«)) (p,w). (4.6.13) 

dp n u*> vp, )) t dM \v, \v, )) dpn kv, ) v > 

To complete the proof: 

1. set this equal to |^ (p, u) 

2. substitute from Shephard's Lemma 

3. define M = e (p, u) (which implies that u = V (p, M)) 

Q.E.D. 

4.7 General Equilibrium Theory 

4.7.1 Walras' law 

Walras . . . 3 

4.7.2 Brouwer's fixed point theorem 

4.7.3 Existence of equilibrium 

4.8 The Welfare Theorems 

4.8.1 The Edgeworth box 

4.8.2 Pareto efficiency 

Definition 4.8.1 A feasible allocation X = (xi, . . . , x#) is Pareto efficient if 
there does not exist any feasible way of reallocating the same initial aggregate 
endowment, J2h=i x /i> which makes one individual better off without making any 
other worse off. 

Definition 4.8.2 X is Pareto dominated by X' = (x' 1? . . . , x'^) if Y,h=i x h = 
J2h=i x 'h> x 'h —h x h V/i and x.' h y h x h for at least one h. 



3 This material still exists only in handwritten form in Alan White's EC3080 notes from 1991- 
2. One thing missing from the handwritten notes is Kakutani's Fixed Point Theorem which should 
be quoted from ?. 
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4.8.3 The First Welfare Theorem 

(See ?.) 

Theorem 4.8.1 (First Welfare Theorem) If the pair (p, X) is an equilibrium (for\ 
given preferences, >Zh> which exibit local non-satiation and given endowments, e/,, 
h = 1 , . . . , H ), then X is a Pareto efficient allocation. 

Proof The proof is by contradiction. Suppose that X is an equilibrium allocation 
which is Pareto dominated by a feasible allocation X'. 

If individual h is strictly better off under X' or x' h >-/,, x^, then it follows that 
individual h cannot afford x' h at the equilibrium prices p or 

p T x' fe > p T X/l = p T e h . (4.8.1) 

The latter equality is just the budget constraint, which is binding since we have 
assumed local non-satiation. 

Similarly, if individual h is indifferent between X and X' or x' h ~^ x^, then it 
follows that 

p T x' /t > p T x h = p T e h , (4.8.2) 

since if x! h cost strictly less than x/j, then by local non- satiation some consumption 
vector near enough to x^ to also cost less than x^ would be strictly preferred to 
Xh and Xfe would not maximise utility given the budget constraint. 
Summing (4.8.1) and (4.8.2) over households yields 

H H H 

p T £ x; > p T J2 x, = p T £ e h , (4.8.3) 

h=l h=l h=l 

(where the equality is essentially Walras' Law). 

But since X' is feasible we must have for each good n 

H H 

h=l h=l 

and, hence, multiplying by prices and summing over all goods, 

TT TT 

p T £x;<p T £e ft . (4.8.4) 

h=l h=l 



But (4.8.4) contradicts the inequality in (4.8.3), so no such Pareto dominant allo- 

Q.E.D. 



cation X^ can exist. 



Before proceeding to the second welfare theorem, we need to say a little bit about 
separating hyperplanes. 
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4.8.4 The Separating Hyperplane Theorem 

Definition 4.8.3 The set jz e $l N : p T z = p T z* | is the hyperplane through z* 
with normal p. 

Note that any hyperplane divides $t N into two closed half-spaces, 

{zef:p T z<p T zj 

and 

{zef :p T z>pV}. 

The intersection of these two closed half-spaces is the hyperplane itself. 

In two dimensions, a hyperplane is just a line; in three dimensions, it is just a 

plane. 

The idea behind the separating hyperplane theorem is quite intuitive: if we take 

any point on the boundary of a convex set, we can find a hyperplane through 

that point so that the entire convex set lies on one side of that hyperplane. We 

will essentially be applying this notion to the upper contour sets of quasiconcave 

utility functions, which are of course convex sets. We will interpret the separating 

hyperplane as a budget hyperplane, and the normal vector as a price vector, so that 

at those prices nothing giving higher utility than the cutoff value is affordable. 

Theorem 4.8.2 (Separating Hyperplane Theorem) If Z is a convex subset of 
R N and z* G Z, z* g int Z, then 3p* f in $l N such that p* T z* < p* T z Vz e Z, 
or Z is contained in one of the closed half-spaces associated with the hyperplane 
through z* with normal p*. 

Proof Not given. See ? 

Q.E.D. 

4.8.5 The Second Welfare Theorem 

(See ?.) 

We make slightly stronger assumptions than are essential for the proof of this 

theorem. This allows us to give an easier proof. 

Theorem 4.8.3 (Second Welfare Theorem) If all individual preferences are strictly^ 
convex, continuous and strictly monotonic, and if X* is a Pareto efficient al- 
location such that all households are allocated positive amounts of all goods 
(x* h 9 > \/g = 1, . . . ,N;h = 1, . . . , H), then a reallocation of the initial ag- 
gregate endowment can yield an equilibrium where the allocation is X*. 
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Proof There are four main steps in the proof. 

1. First we construct a set of utility-enhancing endowment perturbations, and 
use the separating hyperplane theorem to find prices at which no such en- 
dowment perturbation is affordable. 

We need to use the fact (Theorem 3.2. 1) that a sum of convex sets, such as 

X + Y = {x + y:xeX,yeY}, 

is also a convex set. 

Given an aggregate initial endowment x* = S^Li X L we interpret any vec- 
tor of the form z = S^Li *-h — x* as an endowment perturbation. Now 
consider the set of all ways of changing the aggregate endowment without 
making anyone worse off: 



Z = 

1 z e ft N : 3x1 > V 5 , h s.t. u h (x h ) > u h (x* h ) & z 

I h=l 

(4.8.5) 



h=i > 



Z is a sum of convex sets provided that preferences are assumed to be 
convex: 

z=j:x h -{x*} 

h=l 

where 

X h = {x h : u h (x h ) > u h (x* h )} . 

2. Next, we need to show that the zero vector is in the set Z, but not in the 
interior of Z. 

To show that e Z, we just set x^ = x£ and observe that = Y,h=i x t ~ 

x* 4 

The zero vector is not, however, in the interior of Z, since then Z would 
contain some vector, say z*, in which all components were strictly negative. 
In other words, we could take away some of the aggregate endowment of 
every good without making anyone worse off than under the allocation X* . 
But by then giving — z* back to one individual, he or she could be made 
better off without making anyone else worse off, contradicting Pareto opti- 
mality, again using the assumption that preferences are strictly monotonic. 



4 Note also (although I'm no longer sure why this is important) that budget constraints are 
binding and that there are no free goods (by the monotonicity assumption). 
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So, applying the Separating Hyperplane Theorem with z* = 0, we have a 
price vector p* such that = p* T < p* T z Vz e Z. 

Since preferences are monotonic, the set Z must contain all the standard 
unit basis vectors ((1, 0, . . . , 0), &c). This fact can be used to show that 
all components of p* are non-negative, which is essential if it is to be inter- 
preted as an equilibrium price vector. 

3. Next, we specify one way of redistributing the initial endowment in order 
that the desired prices and allocation emerge as a competitive equilibrium. 
All we need to do is value endowments at the equilibrium prices, and redis- 
tribute the aggregate endowment of each good to consumers in proportion 
to their share in aggregate wealth computed in this way. 

4. Finally, we confirm that utility is maximised by the given Pareto efficient 
allocation, X*, at these prices. As usual, the proof is by contradiction: the 
details are left as an exercise. 

Q.E.D. 

4.8.6 Complete markets 

The First Welfare Theorem tells us that competitive equilibrium allocations are 
Pareto optimal if markets are complete. If there are missing markets, then com- 
petitive trading may not lead to a Pareto optimal allocation. We can use the Edge- 
worth Box diagram to illustrate the simplest possible version of this principle. 

4.8.7 Other characterizations of Pareto efficient allocations 

There are a total of five equivalent characterisations of Pareto efficient allocations. 

Theorem 4.8.4 Each of the following is an equivalent description of the set of 
allocations which are Pareto efficient: 

1. by definition, feasible allocations such that no other allocation strictly in- 
creases at least one individual's utility without decreasing the utility of any 
other individual; 

2. by the Welfare Theorems, equilibrium allocations for all possible distribu- 
tions of the fixed initial aggregate endowment; 

3. in two dimensions, allocations lying on the contract curve in the Edgeworth 
box; 
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4. allocations which solve: 5 



H 



max V X h [u h (x h )] (4.8.6) 

{x h :h=l,...,H} h=1 

subject to the feasibility constraints 

H H 

5> h = E e ^ ( 4 - 8 - 7 ) 

h=\ h=l 

for some non-negative weights {\h\h=v 
5. allocations which maximise the utility of a representative agent given by 

H 

£ X h KK)] (4-8.8) 

h=\ 

where {\h}h=i are again any non-negative weights. 

Proof If an allocation is not Pareto efficient, then the Pareto-dominating alloca- 
tion gives a higher value of the objective function in the above problem for all 
possible weights. 

If an allocation is Pareto efficient, then the relative weights for which the above 
objective function is maximized are the ratios of the Lagrange multipliers from 
the problem of maximizing any individual's utility subject to the constraint that 
all other individuals' utilities are unchanged: 

max«i(xi) (4.8.9) 

s.t. 

u h (x h ) = u h (x* h )h = 2,...,H (4.8.10) 

since these two problems will have the same necessary and sufficient first order 
conditions. 

The absolute weights corresponding to a particular allocation are not unique, as 
they can be multiplied by any positive constant without affecting the maximum. 
Different absolute weights (or Lagrange multipliers) arise from fixing different 
individuals' utilities in the last pr oblem, but the relative weights will be the same. 
5 The solution here would be unique if the underlying utility function were concave, since linear 
combinations of concave functions with non-negative weights are concave, and the constraints 
specify a convex set on which the objective function has a unique optimum. This argument can 
not be used with merely quasiconcave utility functions. 
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Q.E.D. 

Note that corresponding to each Pareto efficient allocation there is at least one: 

1. set of non-negative weights defining 

(a) the objective function in 4. and 

(b) the representative agent in 5. 

and 

2. initial allocation leading to the competitive equilibrium in 2. 

4.9 Multi-period General Equilibrium 

In Section 4.2, it was pointed out that the objects of choice can be differentiated 
not only by their physical characteristics, but also both by the time at which they 
are consumed and by the state of nature in which they are consumed. These 
distinctions were suppressed in the intervening sections but are considered again 
in this section and in Section 5.4 respectively. 

The multi-period model should probably be introduced at the end of Chapter 4 
but could also be left until Chapter 7. For the moment this brief introduction is 
duplicated in both chapters. 

Discrete time multi-period investment problems serve as a stepping stone from 
the single period case to the continuous time case. 

The main point to be gotten across is the derivation of interest rates from equilib- 
rium prices: spot rates, forward rates, term structure, etc. 

This is covered in one of the problems, which illustrates the link between prices 
and interest rates in a multiperiod model. 
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Chapter 5 

CHOICE UNDER UNCERTAINTY 



5.1 Introduction 

[To be written.] 

5.2 Review of Basic Probability 

Economic theory has, over the years, used many different, sometimes overlapping, 
sometimes mutually exclusive, approaches to the analysis of choice under uncer- 
tainty. This chapter deals with choice under uncertainty exclusively in a single 
period context. Trade takes place at the beginning of the period and uncertainty 
is resolved at the end of the period. This framework is sufficient to illustrate the 
similarities and differences between the most popular approaches. 
When we consider consumer choice under uncertainty, consumption plans will 
have to specify a fixed consumption vector for each possible state of nature or 
state of the world. This just means that each consumption plan is a random vector. 
Let us review the associated concepts from basic probability theory: probability 
space; random variables and vectors; and stochastic processes. 
Let Cl denote the set of all possible states of the world, called the sample space. 
A collection of states of the world, A C Q, is called an event. 
Let A be a collection of events in Q. The function P : A — ► [0, 1] is a probability 
function if 

1. (a) fie A 

(b) AeA => ft- AeA 

(c) AiGAfori = l,...,oo=> U^i A e A 
(i.e. A is a sigma-algebra of events) 
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and 

2. (a) P (Q) = 1 

(b) P (Q, - A) = 1 - P (A) \/A G A (redundant assumption) 

(c) P(U~iA) = T,TLiP(Ai) when A 1 ,A 2 ,... are pairwise disjoint 
events in A. 

(fi, v4, P) is then called a probability space 

Note that the certainty case we considered already is just the special case of un- 
certainty in which the set fl has only one element. 

We will consider these concepts in more detail when we come to intertemporal 
models. 

Suppose we are given such a probability space. 

The function x : fl — > 5ft is a random variable (r.v.) if \/x G 5ft {uj G Vl : £ (a;) < 
x} G .A, z.e. a function is a random variable if we know the probability that the 
value of the function is less than or equals any given real number. 
The function F~ x : 5ft -»■ [0, 1] : x i-> Pr (5 < x) = P ({w G fi : x (uj) < x}) is 
known as the cumulative distribution function (c.d.f.) of the random variable x. 
The convention of using a tilde over a letter to denote a random variable is com- 
mon in financial economics; in other fields capital letters may be reserved for ran- 
dom variables. In either case, small letters usually denote particular real numbers 
{i.e. particular values of the random variable). 

A random vector is just a vector of random variables. It can also be thought of as 
a vector-valued function on the sample space fl. 

A stochastic process is a collection of random variables or random vectors indexed 
by time, e.g. {x t : t G T} or just {x t } if the time interval is clear from the context. 
For the purposes of this part of the course, we will assume that the index set 
consists of just a finite number of times i.e. that we are dealing with discrete time 
stochastic processes. 

Then a stochastic process whose elements are TV-dimensional random vectors is 
equivalent to an N \T\ -dimensional random vector. 

The (joint) c.d.f. of a random vector or stochastic process is the natural extension 
of the one-dimensional concept. 

Random variables can be discrete, continuous or mixed. The expectation (mean, 
average) of a discrete r.v., x, with possible values x\, x 2 , x%, ... is given by 

oo 

E \x\ = y^ XiPr (x = Xi) . 

i=l 

For a continuous random variable, the summation is replaced by an integral. 
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The covariance of two random variables x and y is given by 

Cov[x,y} = E[(x - e[x])(y - E[y})}. 

The covariance of a random variable with itself is called its variance. 

The expectation of a random vector is just the vector of the expectations of the 

component random variables. 

The variance (variance-covariance matrix) of a random vector is the (symmetric, 

positive semi-definite) matrix of the covariances between the component random 

variables. 

Given any two random variables x and y, we can define a third random variable~e 

by 

~e = y — a — fix. (5.2.1) 

To specify ~e completely, we can either specify a and /3 explicitly or fix them 
implicitly by imposing (two) conditions on~e. We do the latter by insisting 

1. ~e and x are uncorrelated (this is not the same as assuming statistical inde- 
pendence, except in special cases such as bivariate normality) 

2. E[e] = 

It follows that: 

Cov \x, y] 

P = v r-i (5 - 2 - 2) 

Var \x\ 



and 



a = E[y\- pE [x] . (5.2.3) 



But what about the conditional expectation, E [y\x = x]l This is not equal to 
a + fix, as one might expect, unless E [e\x = x] =0. This requires statistical 
independence rather than the assumed lack of correlation. Again, a sufficient 
condition is multivariate normality. 

The notion of the (3ofy with respect to x as given in (5.2.2) will recur frequently. 
The final concept required from basic probability theory is the notion of a mixture 
of random variables. For lotteries which are discrete random variables, with pay- 
offs #i, X2, £3, . . . occuring with probabilities 7Ti, 7T2, tt3, • • • respectively, we will 
use the notation: 

7TiXi © ir 2 x 2 © 7T 3 X 3 © . . . 

Similar notation will be used for compound lotteries (mixtures of random vari- 
ables) where the payoffs themselves are further lotteries. 
This might be a good place to talk about the MVN distribution and Stein's lemma. 
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5.3 Taylor's Theorem: Stochastic Version 

We will frequently use the univariate Taylor expansion as applied to a function 
of a random variable expanded about the mean of the random variable. 1 Taking 
expectations on both sides of the Taylor expansion: 

OO 1 

/(*) = f(E[x]) + J2-J {n) (E[xm-E[x]) n (5.3.1) 



«=i n! 



yields: 



OO 1 

E[f(x)} = f(E[x]) + J2^f {n) (E[x])m n (x), (5.3.2) 



„=2 n! 



where 

In particular, 



and 



m n (x) = E [(x - E[x]) n ] . (5.3.3) 

m\x) = E[(x-E[x]Y] =0 (5.3.4) 

m 2 (x) = E[(x- E[x]) 2 ] = Var [x] (5.3.5) 

m 3 (x) = E[(x- E[x]) 3 ] = Skew [x] (5.3.6) 

m 4 (x) = E[(x- E[x]) 4 ] = Kurt [x] , (5.3.7) 



which allows us to start the summation in (5.3.2) at n = 2 rather than n = 1. 
Indeed, we can rewrite (5.3.2) as 

E[f(x)} = f(E[x}) + h"(E[x])Vax [x] + V(£[x])Skew [x] 

2 b 

I OO I 

+-f""(E[x])Kmt [x] + Y, ^f {n) (E[x])m n (x). (5.3.8) 

n=5 

5.4 Pricing State-Contingent Claims 

This part of the course draws on ?, ? and ?. 

The analysis of choice under uncertaintly will begin by reinterpreting the general 
equilibrium model of Chapter 4 so that goods can be differentiated by the state 
of nature in which they are consumed. Specifically, it will be assumed that the 



'This section will eventually have to talk separately about fcth order and infinite order Taylor 
expansions. 
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underlying sample space comprises a finite number of states of nature. A more 
thorough analysis of choice under uncertainty, allowing for infinite and continuous 
sample spaces and based on additional axioms of choice, follows in Section 5.5. 
Consider a world with M possible states of nature (distinguished by a first sub- 
script), markets for N securities (distinguished by a second subscript) and H con- 
sumers (distinguished by a superscript). 2 

Definition 5.4.1 A state contingent claim or Arrow-Debreu security is a random 
variable or lottery which takes the value I in one particular state of nature and 
the value in all other states. 

Definition 5.4.2 A complex security is a random variable or lottery which can 
take on arbitrary values. The payoffs of a typical complex security will be repre- 
sented by a column vector, y_,- G !K M , where i/ij is the payoff in state i of security 
3- 

The set of all complex securities on a given finite sample space is an M-dimensional| 
vector space and the M possible Arrow-Debreu securities constitute the standard 
basis for this vector space. 

State contingent claims prices are determined by the market clearing equations in 
a general equilibrium model: 

Aggregate consumption in state i = Aggregate endowment in state i. 

Each individual will have an optimal consumption choice depending on endow- 
ments and preferences and conditional on the state of the world. Optimal future 
consumption is denoted 

I X * \ 



x 2 



(5.4.1) 



\ X N I 

If there are N complex securities, then the investor must find a portfolio w 

(wi, . . . , Wn) whose payoffs satisfy 



N 

3=1 

,3 



Let Y be the M x N matrix whose jth column contains the payoffs of the jth 
complex security in each of the M states of nature, i.e. 

Y = ( yi ,y 2 ,...y N ). (5.4.2) 



2 Check for consistency in subscripting etc in what follows. 
3 Or maybe I mean its transpose. 
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Theorem 5.4.1 If there are M complex securities (M = N) and the pay off matrix 
Y is non-singular, then markets are complete. 

Proof Suppose the optimal trade for consumer i state j is Xy — e^. Then can 
invert Y to work out optimal trades in terms of complex securities. 



Q.E.D. 



An (JV + l)st security would be redundant. 

Either a singular square matrix or < N complex securities would lead to incom- 
plete markets. 

So far, we have made no assumptions about the form of the utility function, written 
purely as 

u(x ,x 1 ,x 2 ,.. ■ ,x N ) , 

where x represents the quantity consumed at date and Xi (i > 0) represents the 
quantity consumed at date 1 if state i materialises. 



5.4.1 Completion of markets using options 

Assume that there exists a state index portfolio, Y, yielding different non-zero 
payoffs in each state (i.e. a portfolio with a different payout in each state of nature, 
possibly one mimicking aggregate consumption). WLOG we can rank the states 
so that Yi < Yj if i < j. 

We now present some results, following ?, showing conditions under which trad- 
ing in a state index portfolio and in options on the state index portfolio can lead to 
the Pareto optimal complete markets equilibrium allocation. Now consider com- 
pletion of markets using options on aggregate consumption. 
In real-world markets, the number of linearly independent corporate securities is 
probably less than M. 

However, options on corporate securities may be sufficient to form complete mar- 
kets, and thereby ensure allocational (Pareto) efficiency for arbitrary preferences. 
Further assume that 3 M — 1 European call options on Y with exercise prices 
Yi, Y 2 , . . . ,Ym~i- 

A European call option with exercise price K is an option to buy a security for 
K on a fixed date. 

An American call option is an option to buy on or before the fixed date. 
A put option is an option to sell. 

Revised: December 2, 1998 



CHAPTER 5. CHOICE UNDER UNCERTAINTY 



91 



Here, the original state index portfolio and the M — 1 European call options yield 
the payoff matrix: 



fyi 





2/2 

2/2 -2/1 



2/3 
2/3-2/1 • 
2/3-2/2 • 


2/M 
2/M ~ 2/1 
2/M - 2/2 


Vo 








• 2/M - 2/M- 



\ 



/ 



security Y 
call option 1 
call option 2 



\ 



(5.4.3) 



i / V call option M - 1 / 



and as this matrix is non-singular, we have constructed a complete market. 
Instead of assuming a state index portfolio exists, we can assume identical proba- 
bility beliefs and state-independent utility and complete markets in a similar man- 
ner (see below). 

5.4.2 Restrictions on security values implied by allocational ef- 
ficiency and covariance with aggregate consumption 

J Cuj = aggregate consumption in state uj 
Let \ ft k = {ujeft:CLO = k} 
Let 4>{k) be the value of the claim with payoffs 



Vku 



1 if Cuj = k 
otherwise 



(5.4.4) 



and let the agreed probability of the event Q k (i.e. of aggregate consumption 
taking the value k) be: 

7r(jfc) = ]T itlu (5.4.5) 

u>Gfl k 

By time-additivity and state-independence of the utility function: 

7rum-(c it 



<j)U) 

The no arbitrage condition implies 

m = 



<o( c jo) 



Vw g n 



<{h{k)) ^ 
— — — — 2^ ^^ 



<(/i(fe)) 
Wio(Cio) 



TT(Jfc) 



(5.4.6) 

(5.4.7) 
(5.4.8) 
(5.4.9) 



where fi(k) denotes the i-th individual's equilibrium consumption in those states 
where aggregate consumption equals k. 
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x(0) 


x(l) 


x(2) 


(7 = 1 
(7 = 2 
(7 = 3 


1 

2 
3 



1 

2 




1 


C = L 


L 


L- 1 


L-2 



Table 5.1: Payoffs for Call Options on the Aggregate Consumption 

State-independence of the utility function is required for fi(k) to be well-defined. 
Therefore, an arbitrary security x has value: 



x = e (j)ojxu 


(5.4.10) 


wen 




= E E ( t )ujxuj 


(5.4.11) 


k uj^Qk 




^ u 'i(fi( k )) v- 

= 2^ — ti — ^ Z^ kujxoj 


(5.4.12) 


- £<K*0 £ _ m ^ 


(5.4.13) 




(5.4.14) 



5.4.3 Completing markets with options on aggregate consump- 
tion 

Let x(k) be the vector of payoffs in the various possible states on a European call 

option on aggregate consumption with one period to maturity and exercise price 

k. 

Let {1, 2, . . . , L} be the set of possible values of aggregate consumption C(u>). 

Then payoffs are as given in Table 5.1. 

This all assumes 

1. identical probability beliefs 

2. time-additivity of u 

3. state-independent u 
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5.4.4 Replicating elementary claims with a butterfly spread 

Elementary claims against aggregate consumption can be constructed as follows, 
for example, for state 1, using a butterfly spread: 



[x(0)-x(l)]- [x(l)-x(2)] 



(5.4.15) 



yields the payoff: 



7i\ 

2 
3 


— 


( ° \ 

1 

2 


— 


1 ° ^ 

1 

2 

U-iJ 


— 


( ° \ 


1 

\L-2J 


= 


1 
1 

U/ 


— 


1 
1 



/1\ 




w 



(5.4.16) 



i.e. this replicating portfolio pays 1 iff aggregate consumption is 1, and other- 
wise. 

The prices of this, and the other elementary claims, must, by no arbitrage, equal 
the prices of the corresponding replicating portfolios. 

5.5 The Expected Utility Paradigm 

5.5.1 Further axioms 

The objects of choice with which we are concerned in a world with uncertainty 
could still be called consumption plans, but we will acknowledge the additional 
structure now described by terming them lotteries. 

If there are k physical commodities, a consumption plan must specify a £>dimensional| 
vector, x E yt k , for each time and state of the world. 
We assume a finite number of times, denoted by the set T. 
The possible states of the world are denoted by the set ft. 

So a consumption plan or lottery is just a collection of \T\ fc-dimensional random 
vectors, i.e. a stochastic process. 

Again to distinguish the certainty and uncertainty cases, we let C denote the col- 
lection of lotteries under consideration; X will now denote the set of possible 
values of the lotteries in C. 
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Preferences are now described by a relation on C. We will continue to assume that 
preference relations are complete, reflexive, transitive, and continuous. 
X can be identified with a subset of C, in that each sure thing in X can be identified 
with the trivial lottery that pays off that sure thing with probability (w.p.) 1. 
Although we have moved from a finite-dimensional to an infinite-dimensional 
problem by explicitly allowing a continuum of states of nature, it can be shown 
that the earlier theory of choice under certainty carries through to choice under 
uncertainty, in particular a preference relation can always be represented by a 
continuous utility function on C 

However, we would like utility functions to have a stronger property than conti- 
nuity, namely the expected utility property. 

Axiom 8 (Substitution or Independence Axiom) If a E (0, 1] andp y q, then 

ap © (1 — a) f y aq © (1 — a) f Vf G C. 

Axiom 9 (Archimedian Axiom) Ifp y q y r then 

3a, b E (0, 1) s.t. ap © (1 - a) f y q y bp © (1 - 6) r. 

(The Archimedian axiom is just a generalisation of the continuity axiom.) 

Axiom 10 (Sure Thing Principle) If probability is concentrated on a set of sure 
things which are preferred to q, then the associated consumption plan is also 
preferred to q. 

(The Sure Thing Principle is just a generalisation of the Substitution Axiom.) 

Now let us consider the Allais paradox. 

Suppose 

lflra. >- 0.1£5ra. © 0.89£lm. © O.OlfO. 

Then, unless the substitution axiom is contradicted: 

10 1 

l£lm. > £5rn. © — £0. 

Finally, by the substitution axiom again, 

0.11£lm. © 0.89£0 y 0.1£5m. © 0.9£0. 

If these appears counterintuitive, then so does the independence axiom above. 
One justification for persisting with the independence axiom is provided by ?. 
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5.5.2 Existence of expected utility functions 

A function u: !ft fcx l T l — ► 5ft can be thought of as a utility function on sure things. 

Definition 5.5.1 Let V : C — > 5ft be a utility function representing the preference 
relation y. 

Then y is said to have an expected utility representation if there exists a utility 
function on sure things, u, such that 

V({x t }) = E[u({x t })} 

u ({x t }) dF {it} ({x t }) 

Such a representation will often be called a Von Neumann-Morgenstern (or VNM) 
utility function, after its originators (?), or just an expected utility function. 

Any strictly increasing transformation of a VNM utility function represents the 

same preferences. 

However, only increasing affine transformations: 

/ (x) = a + bx (b > 0) 

retain the expected utility property. 

Proof of this is left as an exercise. 

We will now consider necessary and sufficient conditions on preference relations 

for an expected utility representation to exist. 

Theorem 5.5.1 If X contains only a finite number of possible values, then the 
substitution and Archimidean axioms are necessary and sufficient for a preference 
relation to have an expected utility representation. 

Proof We will just sketch the proof that the axioms imply the existence of an 

expected utility representation; the proof of the converse is left as an exercise. 

For full details, see ?. 

Since X is finite, and unless the consumer is indifferent among all possible choices,! 

there must exist maximal and minimal sure things, say p + and p_ respectively. 

By the substitution axiom, and a simple inductive argument, these are maximal 

and minimal in C as well as in X. (If X is not finite, then an inductive argument 

can no longer be used and the Sure Thing Principle is required.) 

From the Archimedean axiom, it can be deduced that for every other lottery, p, 

there exists a unique V (p) such that 

p~V(p)p + ® (1- V(p))p_. 
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It is easily seen that V represents y. 

Linearity can be shown as follows: 

We leave it as an exercise to deduce from the axioms that if x ~ y and z ~ i then 

itx © (1 — 7r) 2 ~ iry © (1 — 7r) t. 

Define 5 = 7r5 © (1 — 7r) y. 

Then, using the definitions of V (x) and V (y), 

Z ~ ITX © (1 — 7r) £/ 

- vr (V(5)j9 + © (1 - V (x))p_) © (1 - tt) (V (jf)p+ © (1 - V (y))p-) 

= (ttV (x) + (1 - vr) V (y))p + © (tt (1 - V (x)) + (1 - tt) (1 - V (y)))P- 

It follows that 

V (nx © (1 - it) y) = ttV (x) + (1 - vr) V (y) . 

This shows linearity for compound lotteries with only two possible outcomes: by 
an inductive argument, every lottery can be reduced recursively to a two-outcome 
lottery when there are only a finite number of possible outcomes altogether^ p, ^ 

Theorem 5.5.2 For more general C, to these conditions must be added some tech- 
nical conditions and the Sure Thing Principle. 

Proof We will not consider the proof of this more general theorem. It can be 
found in ?. Q>E>D> 

Note that expected utility depends only on the distribution function of the con- 
sumption plan. 

Two consumption plans having very different consumption patterns across states 
of nature but the same probability distribution give the same utility. E.g. if wet 
days and dry days are equally likely, then an expected utility maximiser is indiffer- 
ent between any consumption plan and the plan formed by switching consumption 
between wet and dry days. 

The basic objects of choice under expected utility are not consumption plans but 
classes of consumption plans with the same cumulative distribution function. 
Chapter 6 will consider the problem of portfolio choice in considerable depth. 
This chapter, however, must continue with some basic analysis of the choice be- 
tween one riskfree and one risky asset, following ?. 
Such an example is sufficient to show several things: 

1. There is no guarantee that the portfolio choice problem has any finite or 
unique solution unless the expected utility function is concave. 

2. probably local risk neutrality and stuff like that too. 
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5.6 Jensen's Inequality and SiegePs Paradox 

Theorem 5.6.1 (Jensen's Inequality) The expected value of a (strictly) concave 

function of a random variable is (strictly) less than the same concave function of 

the expected value of the random variable. 

E u (w) < u ( E \W J when u is concave 

Similarly, the expected value of a (strictly) convex function of a random variable 

is (strictly) greater than the same convex function of the expected value of the 

random variable. 

Proof There are three ways of motivating this result, but only one provides a fully 
general and rigorous proof. Without loss of generality, consider the concave case. 

1. One can reinterpret the defining inequality (3.2.1) in terms of a discrete 
random vector x taking on the value x with probability ir and x' with prob- 
ability 1 — 7r: 

Vx^x' EX,ir e (0,1) 

/ (^x + (1 - tt)x') > tt/(x) + (1 - 7r)/(x'), (5.6.1) 

which just says that 

/(£[x]) >£[/(£)]. (5.6.2) 

An inductive argument can be used to extend the result to all discrete r.v.s 
with a finite number of possible values, but runs into problems if the number 
of possible values is either countably or uncountably infinite. 

2. Using a similar approach to that used with the Taylor series expansion in 
(5.3.2), take expectations on both sides of the first order condition for con- 
cavity (3.2.3) given by Theorem 3.2.3, where the two vectors considered are 
the mean E [x] and a generic value x: 

/(x) <f(E [x]) + /' (E [x]) (k-E [x]) . (5.6.3) 

Taking expectations on both sides, the first order term will again disappear, 
once more yielding: 

/(£[*]) >£[/(*)]■ (5-6.4) 

3. One can also appeal to the second order condition for concavity, and the 
second order Taylor series expansion of / around E [x] : 

E[f(x)} = f(E[x]) + ±f"(x*)Vax[x], (5.6.5) 

Revised: December 2, 1998 



98 5.6. JENSEN'S INEQUALITY AND SIEGEL'S PARADOX 

for some x* in the support of x. However, this supposes that x* is fixed, 
whereas in fact it varies with the value taken on by x, and is itself a random 
variable, correlated with x. 

Accepting this (wrong) approximation, if / is concave, then by Theorem 3.2. 4| 
the second derivative is non-positive and the variance is non-negative, so 

E[f{x)\ < f(E[x}). (5.6.6) 

The arguments for convex functions, strictly concave functions and strictly con- 
cave functions are almost identical. 

Q.E.D. 

This result is often useful with functions such as x i— ► In x and £!—►-. 

X 

To get a feel for the extent to which E[f(x)] differs from f(E[x]), we can again 
use the following (wrong) second order Taylor approximation based on (5.3.8): 

E[f(x)} * f(E[x}) + i/"(£[x])Var [x] . (5.6.7) 

This shows that the difference is larger the larger is the curvature of / (as measured 
by the second derivative at the mean of x) and the larger is the variance of x. 
One area in which this idea can be applied is the computation of present values 
based on replacing uncertain future discount factors with point estimates derived 
from expected future interest rates. 
Another nice application of Jensen's Inequality in finance is: 

Theorem 5.6.2 (Siegel's Paradox) Current forward (relative) prices can not all 
equal expected future spot prices. 

Proof Let F t be the current forward price and S t+1 the unknown future spot price. 
If 

E t [S t+1 ] = F t , 

then Jensen's Inequality tells us that 



1 



Ft E t [S t+1 



<E, 



5*, 



m 



except in the degenerate case where S t +i is known with certainty at time t. 

But since the reciprocals of relatives prices are also relative prices, we have shown 

that our initial hypothesis is untenable in terms of a different numeraire. 
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Q.E.D. 

The original and most obvious application of Siegel's paradox is in the case of 
currency exchange rates. In that case, I F t \ and I S t \ are stochastic processes 
representing forward and spot exhange rates respectively. It seems reasonable to 
assume that forward exchange rates are good predictors of spot exchange rates in 
the future, say: 

F t = E t [§ t+1 ] . 

But this is an internally inconsistent hypothesis. 

However, Siegel's paradox applies equally well to any theory which uses current 
prices as a predictor of future values. Another such theory which is enormously 
popular is the Efficient Markets Hypothesis of ?. In its general form, this says 
that current prices should fully reflect all available information about (expected) 
future values. Attempts to make the words fully reflect in any way mathematically 
rigorous quickly run into problems. 

5.7 Risk Aversion 

An individual is risk averse if he or she is unwilling to accept, or indifferent to, 
any actuarially fair gamble. 

(An individual is strictly risk averse if he or she is unwilling to accept any actu- 
arially fair gamble.) 

Figure 1.17.1 goes here. 

i.e. VW and Vp, h\, h 2 such that phi + (1 — p) h 2 = 

W h(y)W +phi + (l-p)h 2 

i.e. 

u (W ) > (>)pu (W + hi) + (l-p)u (W + h 2 ) 

The following interpretation of the above definition of risk aversion is based on 
Jensen's Inequality (see Section 5.6). 

In other words, a (strictly) risk averse individual is one whose VNM utility func- 
tion is (strictly) concave. 

Similarly, a (strictly) risk loving individual is one whose VNM utility function is 
(strictly) convex. 
Finally, a risk neutral individual is one whose VNM utility function is affine. 
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Most functions do not fall into any of these categories, and represent behaviour 
which is locally risk averse at some wealth levels and locally risk loving at other 
wealth levels. 

However, in most of what follows we will find it convenient to assume that indi- 
viduals are globally risk averse. 

Individuals who are globally risk averse will never gamble, in the sense that they 
will never have a bet unless they believe that the expected return on the bet is 
positive. Thus assuming global risk aversion (and rational expectations) rules out 
the existence of lotteries (except with large rollovers) and most other forms of 
betting and gaming. 

We can distinguish between local and global risk aversion. An individual is locally 
risk averse at w if u"(w) < and globally risk averse if u"(w) < \/w. 
Individuals who gamble are not globally risk averse but may still be locally risk 
averse around their current wealth level. 

Cut and paste relevant quotes from Purfield-Waldron papers in here. 
Some people are more risk averse than others; some functions are more concave 
than others; how do we measure this? 

The importance and usefulness of the Arrow-Pratt measures of risk aversion which 
we now define will become clearer as we proceed, in particular from the analysis 
of the portfolio choice problem: 

Definition 5.7.1 (The Arrow-Pratt coefficient of) absolute risk aversion is: 

R A {w) = -u"(w)/u'(w) 

which is the same for u and au + b. 

Note that this varies with the level of wealth. 

u' (w) alone is meaningless, as u and v! can be multiplied by any positive constant 
and still represent the same preferences. However, the above ratio is independent 
of the expected utility function chosen to represent the preferences. 

Definition 5.7.2 (The Arrow-Pratt coefficient of) relative risk aversion is: 

R R (w) = wR A (w) 

The utility function u exhibits increasing (constant, decreasing) absolute risk 
aversion (IARA, CARA, DARA) <^> 

R' A {w) > (=,<)0Vw. 

The utility function u exhibits increasing (constant, decreasing) relative risk aver- 
sion (IRRA, CRRA, DRRA) ^^ 
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R' R (w) > (=,<)0Vw. 
Note: 

• CARA or IARA ^ IRRA 

• CRRA or DRRA =► DARA 

Here are some examples of utility functions and their risk measures: 

• Quadratic utility (IARA, IRRA): 



b> 0; 



u(w) = 


b 2 

w w , 

2 


u'{w) = 


1 — bw; 


u"{w) = 


-b<0 


Ra{w) = 


b 


1 — bw 


<IRa{w) 


b 2 



dw (1 — bw) 2 

In this case, marginal utility is positive and utility increasing if and only if 
w < 1/6. 1/6 is called the bliss point of the quadratic utility function. For 
realism, 1/6 should be rather large and thus 6 rather small. 

Negative exponential utility (CARA, IRRA): 



u{w) 


= 


-e~ bw , 6 > 


u'{w) 


= 


be~ bw > 0; 


u"(w) 


= 


-b 2 e~ hw < 


Ra(w) 


= 


6 


dR A {w) 




n 



dw 



Narrow power utility (CRRA, DARA): 



r> 

u(w) = -w 1 '^, w > 0, B > 0, B ^ 1 

B — 1 
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The proofs in this case are left as an exercise. The solution is roughly as 
follows: 4 

u(z) = - — -z 1 '^, z > 0,B > 
B — 1 

u'(z) = Z~z 

1 

B 



u"(z) = -^z-^- 1 



Ra{z) = \,z~ x 



Rr{z) 



B 
1 

B 



dR A (z) = _]_ z -2 q 

dz B 

cIRr(z) 



dz 
Extended power utility (): 







.5 



Note that a risk neutral investor will seek to invest his or her entire wealth in the 
asset with the highest expected return. If all investors are risk neutral, then prices 
will adjust in equilibrium so that all securities have the same expected return. If 
there are risk averse or risk loving investors, then there is no reason for this result 
to hold, and in fact it almost certainly will not. 



5.8 The Mean- Variance Paradigm 

Three arguments are commonly used to motivate the mean-variance framework 
for analysis of the portfolio choice problem: 

1. Taylor- approximated utility functions (see Section 2.8): 

U (W) = u(E[W])+u'(E[W])(W -E[W]) 

+ -u"(E[W])(W-E[W]) 2 + R 3 (5.8.1) 

OO 1 

R 3 = y -u^(E\W])(W - E\W]) n (5.8.2) 



4 Add comments: u'(w) > 0, u"(w) < and u'(w) — > 1/w as B — * 1, so that u(w) —> Inw. 
5 Check ? for details of this one. 
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which implies: 

E[u(W)} = u(E[W]) + )-u"{E[W\)a 2 {W) + E[R z ] (5.8.3) 
where 

OO 1 

E[R 3 ] = Y -u {n \E\W])m n (W) (5.8.4) 

n=3 nl ' 

It follows that the sign of the nth derivative of the utility function determines 
the direction of preference for the nth central moment of the probability 
distribution of terminal wealth. For example, a positive third derivative im- 
plies a preference for greater skewness. It can be shown fairly easily that 
an increasing utility function which exhibits non-increasing absolute risk 
aversion has a non-negative third derivative. 



2. Quadratic utility: 

E[u(W)] = E[W] - ^E[W 2 ] (5.8.5) 

= E[W]- b -((E[W]) 2 + a 2 (W)) (5.8.6) 

= u(E[W}) - -Var[W] (5.8.7) 

3. Normally distributed asset returns. 

Note that the expected utility axioms are neither necessary nor sufficient to guar- 
antee that the Taylor approximation to n moments is a valid representation of the 
utility function. Some counterexamples of both types are probably called for here, 
or maybe can be left as exercises. 

Extracts from my PhD thesis can be used to talk about signing the first three 
coefficients in the Taylor expansion, and to speculate about further extensions to 
higher moments. 



5.9 The Kelly Strategy 

In a multi-period, discrete time, investment framework, investors will be con- 
cerned with both growth (return) and security (risk). There will be a trade-off 
between the two, and investors will be concerned with finding the optimal trade- 
off. This, of course, depends on preferences, but some useful benchmarks exist. 
There are three ways of measuring growth: 

1 . the expected wealth at time t 
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2. the expected rate of growth of wealth up to time t 

3. the expected first passage time to reach a critical level of wealth, say £ lm 
Likewise, there are three ways of measuring security: 

1. the probability of reaching a critical level of wealth, say £lm, by a specific 
timet 

2. the probability that the wealth process lies above some critical path, say 
above b t at time t, t = 1, 2, . . . 

3. the probability of reaching some goal U before falling to some low level of 
wealth L 

It can be shown that the ? strategy both maximises the long run exponential 
growth rate and minimises the expected time to reach large goals. 
All this is also covered in ? which is reproduced in ?. 

5.10 Alternative Non-Expected Utility Approaches 

Those not happy with the explanations of choice under uncertainty provided within| 
the expected utility paradigm have proposed various alternatives in recent years. 
These include both vague qualitative waffle about fun and addiction and more for- 
mal approaches looking at maximum or minimum possible payoffs, irrational ex- 
pectations, etc. Before proceeding, the reader might want to review Section 3.2. 
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Chapter 6 



PORTFOLIO THEORY 



6.1 Introduction 

Portfolio theory is an important topic in the theory of choice under uncertainty. It 
deals with the problem facing an investor who must decide how to distribute an 
initial wealth of, say, Wq among a number of single -period investment opportuni- 
ties, called securities or assets. 

The choice of portfolio will depend on both the investor's preferences and his 
beliefs about the uncertain payoffs of the various securities. A mutual fund is just 
a special type of (managed) portfolio. 

The chapter begins by considering some issues of definition and measurement. 
Section 6.3 then looks at the portfolio choice problem in a general expected utility 
context. Section 6.4 considers the same problem from a mean-variance perspec- 
tive. This leads on to a discussion of the properties of equilibrium security returns 
in Section 6.5. 



6.2 Notation and preliminaries 

6.2.1 Measuring rates of return 

Good background reading for this section is ?. 

A rate of interest (growth, inflation, &c.) is not properly defined unless we state 

the time period to which it applies and the method of compounding to be used. 

2% per annum is very different from 2% per month. 

Table 6.1 illustrates what happens to £100 invested at 10% per annum as we 

change the interval of compounding. The final calculation in the table uses the 
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Compounded 

Annually £100 -> £110 

Semi-annually £100 -► £100 x (1.05) 2 = £110.25 

Quarterly £100 -► £100 x (1.025) 4 = £110.381. . . 



Monthly £100 -► £100 x (l + ^f) 12 = £110.471. 



Weekly £100 -► £100 x (l + § ) =£110.506... 

/ ln \ 365 

Daily £100 -► £100 x (l + |gj =£110.515... 

Continuously £100 -► £100 x e 010 = £110.517. . . 

Table 6.1: The effect of an interest rate of 10% per annum at different frequencies 
of compounding. 



fact that 



lim 1 + - 

n->oo \ n 



where e » 2.7182. . . 

This is sometimes used as the definition of e but others prefer to start with 

r 2 r 3 ™ri 

e = 1+r+ 2T + 37 + - = gjf 

where n\ = n (n — 1) (n — 2) . . . 3.2.1 and 0! = 1 by convention. 
There are five concepts which we need to be familiar with: 

1. Discrete compounding: 

1 + n) P °- 
We can solve this equation for any of five quantities given the other four: 

(a) present value 

(b) final value 

(c) implicit rate of return 

(d) time 

(e) interval of compounding 
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2. Continuous compounding: 

Pt = e rt P . 

We can solve this equation for any of four quantities given the other three 
((a)-(d) above). 

Note also that the exponential function is its own derivative; that y = 1 +r is 
the tangent to y = e r at r = and hence that e r > 1 + r Vr. In other words, 
given an initial value, continuous compounding yields a higher terminal 
value than discrete compounding for all interest rates, positive and negative, 
with equality for a zero interest rate only. Similarly, given a final value, 
continuous discounting yields a higher present value than does discrete. 

3. Aggregating and averaging returns across portfolios (f w = w T f) and over 
time: 

Pt = e rt P 

Pt + At = e rAt P t 

lnP t+At -lnP t 



r(s) 



At 
dlnP, 



ds 



r(s)ds = I dlnP s = lnP t -lnP 
Jo 

P t = P Q eti r(s)ds 



Discretely compounded rates aggregate nicely across portfolios; continu- 
ously compounded rates aggregate nicely across time. 

4. The (net) present value (NPV) of a stream of cash flows, 

Po,Pi, ■ ■ ■ ,Pt 
is 

Pq , Pi , Pi . Pt 



NPV(r) 



\0 ' /-i , \1 ' /-i , \2 i - - - i T - 



(l + r) u (1+r) 1 (l + r) z (l + r) J 

5. The internal rate of return (IRR) of the stream of cash flows, 

Pq, Pi, ■ ■ ■ , Pt, 

is the solution of the polynomial of degree T obtained by setting the NPV 
equal to zero: 

Pq P\ , P2 Pt 



\0 ~ /-> , \1 ~ /1 , \2 ~ ■ ■ ■ ~ /-, , nT 



0. 



(l+r) u (1 + r) 1 (l+r) z (l + r) J 
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W = the investor's initial wealth 

W\ = the investor's desired expected final wealth 

iV = number of risky assets 

77 = return on the riskfree asset 

f j G 3? = return on jth risky asset 

e = E[t] G Jt^ = vector of expected returns 

V G 3? 7Vx7V = variance-covariance matrix of returns 

1 = (1,1,..., 1) 

= iV-dimensional vector of Is 

Wj G 3? = amount invested in jth risky asset 

weK 1 = (wi,...,w n ) 

f w = w T f = return on the portfolio w 

fi = E[f w ] = H 3 - G 5ft = the investor's desired expected return 



Table 6.2: Notation for portfolio choice problem 

In general, the polynomial defining the IRR has T (complex) roots. Condi- 
tions have been derived under which there is only one meaningful real root 
to this polynomial equation, in other words one corresponding to a positive 
IRR. 1 

Consider a quadratic example. 

Simple rates of return are additive across portfolios, so we use them in one period 

cross sectional studies, in particular in this chapter. 

Continuously compounded rates of return are additive across time, so we use them 

in multi-period single variable studies, such as in Chapter 7. 

Consider as an example the problem of calculating mortgage repayments. 

6.2.2 Notation 

The investment opportunity set for the portfolio choice problem will generally 
consist of N risky assets. From time to time, we will add a riskfree asset. The 
notation used throughout this chapter is set out in Table 6.2. The presentation 
is in terms of a single period problem, and the unconditional distribution of re- 
turns. The analysis of the multiperiod, infinite horizon, discrete time problem, 
concentrating on the conditional distribution of the next period's returns given 
this period's, is quite similar. 2 



'These conditions are discussed in ?. 
2 Make this into an exercise. 
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Definition 6.2.1 w is said to be a unit cost or normal portfolio if its weights sum 
to 1 (w T l = 1). 

The portfolio held by an investor with initial wealth Wo can be thought of either 
as a w with w T l = Wo or as the corresponding normal portfolio, ^w. It will 
hopefully be clear from the context which meaning of 'portfolio' is intended. 

Definition 6.2.2 w is said to be a zero cost or hedge portfolio if its weights sum 
toO(w T l = 0). 

The vector of net trades carried out by an investor moving from the portfolio w 
to the portfolio wi can be thought of as the hedge portfolio w x — w . 

Definition 6.2.3 Short-selling a security means owning a negative quantity of it. 

In practice short-selling means promising (credibly) to pay someone the same cash 
flows as would be paid by a security that one does not own, always being prepared, 
if required, to pay the current market price of the security to end the arrangement. 
Thus when short-selling is allowed w can have negative components; if short- 
selling is not allowed, then the portfolio choice problem will have non-negativity 
constraints, u>i > for i = 1, . . . , N. 
A number of further comments are in order at this stage. 

1. In the literature, initial wealth is often normalised to unity (Wo = 1), but 
the development of the theory will be more elegant if we avoid this. 

2. Since we will be dealing on occasion with hedge portfolios, we will in future 
avoid the concepts of rate of return and net return, which are usually thought 
of as the ratio of profit to initial investment. These terms are meaningless 
for a hedge portfolio as the denominator is zero. 

Instead, we will speak of the gross return on a portfolio or the portfolio pay- 
off. This can be defined unambiguously as follows. There is no ambiguity 
about the payoff on one of the original securities, which is just the gross 
return per pound invested. The payoff on a unit cost or normal portfolio 
is equivalent to the gross return. It is just w T r. The payoff on a zero cost 
portfolio can also be defined as w T f . 

Note that where we initially worked with net rates of return (-p- — 1), we 
will deal henceforth with gross rates of return (|r). 
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6.3 The Single-period Portfolio Choice Problem 

6.3.1 The canonical portfolio problem 

Unless otherwise stated, we assume that individuals: 

1. have von Neumann-Morgenstern (VNM) utilities: 

i.e. preferences have the expected utility representation: 

v(z) = E[u(z)] 

u{z)dF- z {z) 

where v is the utility function for random variables (gambles, lotteries) 
and u is the utility function for sure things. 

2. prefer more to less 



i.e. u is increasing: 

3. are (strictly) risk-averse 
i.e. u is strictly concave: 



u'(z) > V2 (6.3.1) 



u"(z) <0 Vz (6.3.2) 



Date investment: 

• Wj (pounds) in jth risky asset, j = 1, . . . , N 

• {Wo _ J2j Wj) in risk free asset 
Date 1 payoff: 

• Wjfj from jth risky asset 

• {Wo — J2j w j) r f from risk free asset 

It is assumed here that there are no constraints on short-selling or borrowing 

(which is the same as short-selling the riskfree security). 

The solution is found as follows: 

Choose WjS to maximize expected utility of date 1 wealth, 



W = {W -Y J Wj) r f + Y. u '- r 
j j 

= W r f + Y, w j{rj- r f 



3' 3 
3 



Revised: December 2, 1998 



CHAPTER 6. PORTFOLIO THEORY 111 

i.e. 

max / (wi, . . .,w N ) = E[u(W r f + J2 W J^J ~ r /))] 

The first order conditions are: 

E[u'(W)(fj - r f )\ = Vj. (6.3.3) 

The Hessian matrix of the objective function is: 

A = E [a" (W) (f - r/1) (f - r f l) T ] . (6.3.4) 

h T Ah < Vh ^ Oat if and only if u" < 0. Since we have assumed that investor 
behaviour is risk averse, A is a negative definite matrix and, by Theorem 3.2.4, 
/ is a strictly concave (and hence strictly quasiconcave) function. Thus, under 
the present assumptions, Theorems 3.3.3 and 3.3.4 guarantee that the first order 
conditions have a unique solution. The trivial case in which the random returns 
are not really random at all can be ignored. 

The rest of this section should be omitted until I figure out what is going on. 
Another way of writing (6.3.3) is: 

E[u'{W)fj\ = E[u'(W)]r f \fj, (6.3.5) 

or 

Cov [w'(W0,fJ + E[u'(W)]E[fj] = E[u'(W)]r f Vj, (6.3.6) 

or 



Cov \u'(W),f u 

E[rj ~ r f ] = \ - J Vj, (6.3.7) 

Suppose pj is the price of the random payoff Xj. Then fj = — and 



Pj = E 



u'{W) 

E[u'(W)]r f 



Vj. (6.3.8) 



In other words, payoffs are valued by taking their expected present value, using 
the stochastic discount factor " ' ;- 1 — , which ends up being the same for all 

J E[u'(W)]rf r ° 

investors. Practical corporate finance and theoretical asset pricing models to a 
large extent are (or should be) concerned with analysing this discount factor. 
(We could consider here the explicit example with quadratic utility from the prob- 
lem sets.) 
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6.3.2 Risk aversion and portfolio composition 

For the moment, assume only one risky asset (N = 1). 

We first consider the concept of local risk neutrality. The optimal investment in 

the risky asset is positive 

The objective function is increasing at a = 

/'(o)>0 (6.3.9) 

E[u (W r f ) (f - 77)] > (6.3.10) 

u'(W o r f )E[(r-r f )]>0 (6.3.11) 



E[r] > E[r f ] = r f 

This is the property of local risk neutrality — a risk averse investor will always 
prefer a little of a risky asset paying a higher expected return than ry to none of 
the risky asset. 

Definition 6.3.1 Let f:X—> 3? ++ be a positive-valued function defined on X C 
3ft+ + - Then the elasticity of / with respect to Xi at x* is 

x* df 



f (x) dx,, 



(x*). 



Roughly speaking, the elasticity is just Q ^£ or the slope of the graph of the 
function on log-log graph paper. 

A function is said to be inelastic when the absolute value of the elasticity is less 
than unity; and elastic when the absolute value of the elasticity is greater than 
unity. The borderline case is called a unit elastic function. 

One useful application of elasticity is in analysing the behaviour of the total rev- 
enue function associated with a particular inverse demand function, P(Q). We 
have: 

dP (Q) Q dP 

-^r - q ^ +p ( "- 12) 

= FI1 + -). (6.3.13) 
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Hence, total revenue is constant or maximised or minimised where elasticity equals| 
— 1; increasing when elasticity is less than —1 (demand is elastic); and decreasing 
when elasticity is between and —1 (demand is inelastic). 
Now consider other properties of asset demands (assuming that E [f] > r/): 

• DARA => risky asset normal ( Jp- > 0) 

• CARA =► (^ = 0) 

• IARA =>• risky asset inferior ( -^- < 0] 

Define the wealth elasticity of demand for the risky asset to be 

V = WTT- (6.3.14) 



a dW< 



Then we have 



d (4) W ° 







da 
dW 



Note that 



dW Wl 

a 



sign f^-j =sign (77 - 1) (6.3.15) 



since by assuming a positive expected risk premium on the risky asset we guaran- 
tee (by local risk-neutrality) that a is positive. 

• DRRA ^> increasing proportion of wealth invested in the risky asset (77 > 1) 

• CRRA =>■ constant proportion of wealth invested in the risky asset (77 = 1) 

• IRRA ^> decreasing proportion of wealth invested in the risky asset (77 < 1) 

Theorem 6.3.1 DARA => RISKY ASSET NORMAL 

Proof By implicit differentiation of the now familiar first order condition (6.3.3), 
which can be written: 

E[u'(W r f + a(r - r f ))(r - 77)] = 0, (6.3.16) 

we have 

da _ E{u"(W)(r-r f )]r f ^^ 



dW -E[u"(W)(f-r f ) 2 ] 
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By concavity, the denominator is positive. Therefore: 

sign (da/dWo) = sign {E[u"(W)(f - 77)]} (6.3.18) 

We will show that both are positive. 
For decreasing absolute risk aversion: 3 

r>r f =* R A (W) < R A (W r f ) 
r<r f =► i^W > ^(Wor/) 

Multiplying both sides of each inequality by —u'(W)(r — 77) gives respectively: 

u"(W)(r - 77) > -R A (W r f )u'(W)(f - 77) (6.3.19) 

in the event that f > 77, and 

u"(W)(r - 77) > -i? A (W r / )w / (iy)(f - r 7 ) (6.3.20) 

(the same result) in the event that f < 77 
Integrating over both events implies: 

E[u"(W0(f - r/)] > -i?^(W r / )E[w / (iy)(f - r/)], (6.3.21) 

provided that f > 77 with positive probability. 

The RHS of inequality (6.3.21) is at the optimum, hence the LHS is positive as 

claimed. 

Q.E.D. 

The other results are proved similarly (exercise!). 

6.3.3 Mutual fund separation 

Commonly, investors delegate portfolio choice to mutual fund operators or man- 
agers. We are interested in conditions under which large groups of investors will 
agree on portfolio composition. For example, all investors with similar utility 
functions might choose the same portfolio, or all investors with similar probabil- 
ity beliefs might choose the same portfolio. More realistically, we may be able to 
define a group of investors whose portfolio choices all lie in a subspace of small 
dimension (say 2) of the iV-dimensional portfolio space. The first such result is 
due to ?. 



Think about whether separating out the case of f = 77 is necessary. 
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Theorem 6.3.2 3 Two fund monetary separation 

i.e. Agents with different wealths (but the same increasing, strictly concave, VNM 

utility) hold the same risky unit cost portfolio, p* say, (but may differ in the mix of 

the riskfree asset and risky portfolio) 

i.e. V portfolios p, wealths Wo, 3X s.t. 

E [u (W r f + XW p* T (f - 77I))] > E [u (w r f + p T (f - 77I))] (6.3.22) 



Risk-tolerance (1/ Ra{z)) is linear (including constant) 

i.e. 3 Hyperbolic Absolute Risk Aversion (HARA, incl. CARA) 

i.e. the utility function is of one of these types: 

• Extended power: u(z) = ( C +i )B (A + Bz) c+1 

• Logarithmic: u{z) = ln(A + Bz) 

• Negative exponential: u(z) = — ^exp{Bz} 

where A, B and C are chosen to guarantee v! > 0, u" < 0. 
i.e. marginal utility satisfies 

u'{z) = {A + Bz) c or u'{z) = Aexp{Bz} (6.3.23) 

where A, B and C are again chosen to guarantee v! > 0, u" < 0. 

Proof The proof that these conditions are necessary for two fund separation is 
difficult and tedious. The interested reader is referred to ?. 
We will show that u'(z) = (A + Bz) c is sufficient for two-fund separation. 
The optimal dollar investments u>j are the unique solution to the first order condi- 
tions: 

= E[u'(W)—] (6.3.24) 

OWi 

= E[{A + BWf{fi-r f )} (6.3.25) 

= E[(A + BW r f + Y, Bwjifj - r f )f(f t - 77)], (6.3.26) 

3 

or equivalently to the system of equations 

E[{1 + ? A + BWorf ^ ~ Tf))C{h ~ r/)] = ° (63 - 27) 
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or 

E[{1 + ^x j (f j - r f )) c {f t - r f )\ = (6.3.28) 

3 

where 

Bwj 

Xj — 



A + BW r f 

The unique solutions for Xj are clearly independent of Wo which does not appear 
in (6.3.28). 

Since A and B do not appear either, the unique solutions for Xj are also indepen- 
dent of those parameters. 
However, they do depend on C. 
But the risky portfolio weights are 



Wi _ Bw t /(A + BW r f ) 



(6.3.29) 
(6.3.30) 



T. jWj EjBwj/iA + BWorf) 

x% 

Z^j Xj 

and so are also independent of initial wealth. 

Since the dollar investment in the j\h risky asset satisfies: 

w 3 = x j(-g+ W orf) (6.3.31) 

we also have in this case that the dollar investment in the common risky portfolio 

is a linear function of the initial wealth. The other sufficiency proofs are similar 

and are left as exercises. 

Q.E.D. 

Some humorous anecdotes about Cass may now follow. 

6.4 Mathematics of the Portfolio Frontier 

6.4.1 The portfolio frontier in $l N : 
risky assets only 

The portfolio frontier 

Definition 6.4.1 The (mean- variance) portfolio frontier is the set of solutions to 
the mean-variance portfolio choice problem faced by a (risk-averse) investor with 
an initial wealth of Wo who desires an expected final wealth of at least W\ = fiWo 
(or, equivalently, an expected rate of return of p,), but with the smallest possible 
variance of final wealth. 
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The mean- variance frontier can also be called the two -moment portfolio frontier, 
in recognition of the fact that the same approach can be extended (with difficulty) 
to higher moments. 

The mean- variance portfolio frontier is a subset of the portfolio space, which is 
just R N . However, introductory treatments generally present it (without proof) as 
the envelope function, in mean- variance space or mean-standard deviation space 
(3? + x 5ft), of the variance minimisation problem. 

Definition 6.4.2 w is a frontier portfolio 

its return has the minimum variance among all portfolios that have the same cost, 
w T l, and the same expected payoff , w T e. 

We will begin by supposing that all assets are risky. Formally, the frontier port- 
folio corresponding to initial wealth W and expected return \x (expected terminal 
wealth jjWq) is the solution to the quadratic programming problem: 

minw Vw (6.4.1) 

w 

subject to the linear constraints: 

w T l = W (6.4.2) 



and 



w T e > W 1 = fiW . (6.4.3) 



The first constraint is just the budget constraint, while the second constraint states 

that the expected rate of return on the portfolio is at least the desired mean return 

fi. 

The frontier in this case is the set of solutions for all values of W and W-y (or //) to 

this variance minimisation problem, or to the equivalent maximisation problem: 

max — w Vw (6.4.4) 

w 

subject to the same linear constraints (6.4.2) and (6.4.3). 

The properties of this two-moment frontier are well known, and can be found, 
for example, in ? or ?. The notation here follows ?. The derivation of the mean- 
variance frontier is generally presented in the literature in terms of portfolio weight 
vectors or, equivalently, assuming that initial wealth, Wo, equals 1. This assump- 
tion is not essential and will be avoided. 
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The solution 

The inequality constrained maximisation problem (6.4.4) is just a special case 
of the canonical quadratic programming problem considered at the end of Sec- 
tion 3.5, except that it has explicitly one equality constraint and one inequality 
constraint. 
To avoid degeneracies, we require: 

1. that not every portfolio has the same expected return, i.e. 

e^E[h]l, (6.4.5) 

and in particular that N > 1. 

2. that the variance-covariance matrix, V, is (strictly) positive definite. We 
already know from (1.12.4) that V must be positive semi-definite, but we 
require this slightly stronger condition. To see why, suppose 

3w ^ N s.t. w T Vw = (6.4.6) 

Then 3 a portfolio whose return w T f = f w has zero variance. 

This implies that f w = r (say) w.p.l or, essentially, that this portfolio is 
riskless. 

Arbitrage will force the returns on all riskless assets to be equal in equilib- 
rium, so this situation is equivalent economically to the introduction of a 
riskless asset later. 

In the portfolio problem, the place of the matrix A in the canonical quadratic 

programming problem is taken by the (symmetric) negative definite matrix, —V, 

which is just the negative of the variance-covariance matrix of asset returns; g 1 = 

1 T and a.\ = W ; and g 2 = e T and a 2 = W\. (6.4.5) guarantees that the 2 x A 

matrix G is of full rank 2. 

The parallels are a little fuzzy in the case of the budget constraint since it is really 

an equality constraint. 

(3.5.39) says that the optimal w is a linear combination of the two columns of the 

A x 2 matrix 

V- 1 G T fGV" 1 G TX " 



with columns weighted by initial wealth W and expected final wealth, W\. 
We will call these columns g and h and write the solution as 

w = I¥ g + Wih = W (g + (ib.) . (6.4.7) 
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The components of g and h are functions of the means and variances of security 
returns. Thus the vector of optimal portfolio proportions, 

— w = g + yUh, (6.4.8) 

Wo 

is independent of the initial wealth W . 

It is easy to see the economic interpretation of g and h: 

• g is the frontier portfolio corresponding to W = 1 and W\ = 0. In other 
words, it is the normal portfolio which would be held by an investor whose 
objective was to (just) go bankrupt with minimum variance. 

• Similarly, h is the frontier portfolio corresponding to W = and W\ = 1. 
In other words, it is the hedge portfolio which would be purchased by a 
variance-minimising investor in order to increase his expected final wealth 
by one unit. 

Alternatively, (3.5.35) says that the optimal w is a linear combination of the two 
columns of the TV x 2 matrix 

^V" 1 G T = ( IV-l IV^e 

with columns weighted by the Lagrange multipliers corresponding to the two con- 
straints. We will call the Lagrange multipliers 27/C and 2X/A respectively, where 
we define: 

A = l T V _1 e = e T V _1 l (6.4.9) 

B = e T \~ 1 e > (6.4.10) 

C = 1 T V _1 1 > (6.4.11) 

and 

D = BC-A 2 (6.4.12) 

and the inequalities follow from the fact that V -1 (like V) is positive definite. 
This allows the solution to be written as: 

w = ^(V- 1 l) + ^(V- 1 e). (6.4.13) 

^(V _1 1) and i(V _1 e) are both unit portfolios, so 7 + A = W . We know that 
for the portfolio which minimises variance for a given initial wealth, regardless 
of expected final wealth, the corresponding Lagrange multiplier, A = 0. Thus 
^(V -1 1) is the global minimum variance portfolio with cost Wo (which in fact 
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equals 7 in this case) and ^(V _1 l) is the global minimum variance unit cost 

portfolio, which we will denote w M vp- 

In fact, we can combine (6.4.7) and (6.4.13) and write the solution as: 

w = W (w MV p +(u-—jh\. (6.4.14) 

The details are left as an exercise. 4 

The set of solutions to this quadratic programming problem for all possible (Wo, Wi)| 
combinations (including negative Wo) is the vector subspace of the portfolio space,| 
which is generated either by the vectors g and h or by the vectors w MV p and h (or 
by any pair of linearly independent frontier portfolios). 

In R N , the set of unit cost frontier portfolios is the line passing thru g, parallel to 
h. 

It follows immediately that the frontier (like any straight line in $l N ) is a convex 
set, and can be generated by linear combinations of any pair of frontier portfolios 
with weights of the form a and (1 — a). 

An important exercise at this stage is to work out the means, variances and co- 
variances of the returns on w M vp, g and h. They will drop out of the portfolio 
decomposition below. 
The portfolio weight vectors g and h are 

g = I^V-^-ACV-'e)] (6.4.15) 

h = hc(y- 1 e)-A(V- 1 l)} (6.4.16) 



We have 



zr 



Var[f g ] = g T Vg = ^ (6.4.17) 

Var[f h ] = h T Vh=^ (6.4.18) 



from which it follows that D > 0. 



Orthogonal decomposition of portfolios 



At this stage, we must introduce a scalar product on the portfolio space, namely 
that based on the variance-covariance matrix V. Since V is a non-singular, pos- 
itive definite matrix, it defines a well behaved scalar product and all the standard 
results on orthogonal projection (&c.) from linear algebra are valid. 



4 At least for now. 
Revised: December 2, 1998 



CHAPTER 6. PORTFOLIO THEORY 121 

Two portfolios wi and w 2 are orthogonal with respect to this scalar product 

w[Vw 2 = 

Cov wj~f, wjf = 

the random variables representing the returns on the portfolios are uncorrected. 
Thus, the terms 'orthogonal' and 'uncorrected' may legitimately, and shall, be 
applied interchangeably to pairs of portfolios. Furthermore, the squared length of 
a weight vector corresponds to the variance of its returns. 

Note that w M vp and h are orthogonal vectors in this sense. In fact, we have the 
following theorem: 

Theorem 6.4.1 Ifw is a frontier portfolio and u is a zero mean hedge portfolio, 
then w and u are uncorrelated. 

Proof There is probably a full version of this proof lost somewhere but the fol- 
lowing can be sorted out. 

Since w M vp is collinear with V _1 l, it is orthogonal to all portfolios w for which 
w T VV _1 l = or in other words to all portfolios for which w T l = 0. But these 
are precisely all hedge portfolios, including h. 

Similarly, any portfolio collinear with V _1 e is orthogonal to all portfolios with 
zero expected return, since w T VV"'e = or in other words w T e = 0. 

Q.E.D. 

Some pictures are in order at this stage. 

For A^ = 3, in the set of portfolios costing Wq (the Wo simplex), the iso-variance 
curves are concentric ellipses, the iso-mean curves are parallel lines, and the solu- 
tions for different /is (or Wis) are the tangency points between these ellipses and 
lines, which themselves lie on a line orthogonal (in the sense defined above) to 
the iso-mean lines. The centre of the concentric ellipses is at the global minimum 
variance portfolio corresponding to Wo, W^oWmvp- A similar geometric interpre- 
tation can be applied in higher dimensions. 

? has some nice pictures of the frontier in portfolio space, as opposed to mean- 
variance space. 

At this stage, recall the definition of f3 in (5.2.2). 

We will now derive an orthogonal decomposition of a portfolio q into two frontier 
portfolios and a zero-mean zero-cost portfolio and prove that the coefficients on 
the two frontier portfolios are the /3s of q with respect to those portfolios and sum 
to unity. 
We can always choose an orthogonal basis for the portfolio frontier. 
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For any frontier portfolio p ^ w M vp, there is a unique unit cost frontier portfolio 

z p which is orthogonal to p. 

Another important exercise is to figure out the relationship between E [f p ] and 

Any two frontier portfolios span the frontier, in particular any unit cost p ^ w MV p 
and z p (or the original basis, w MV p and h). 

Any (frontier or non-frontier) portfolio q with non-zero cost W can be written in 
the form f q + u q where 

f q = VMg + £[r q ]h) (6.4.19) 

= W (/3 qp p + (1 - /? qp )z p ) (say) (6.4.20) 

is the frontier portfolio with expected return E[f^\ and cost Wq and u q is a hedge 
portfolio with zero expected return. Geometrically, this decomposition is equiva- 
lent to the orthogonal projection of q onto the frontier. 

Theorem 6.4.1 shown that any portfolio sharing these properties of u q is uncorre- 
cted with all frontier portfolios. 5 

If p is a unit cost frontier portfolio {i.e. the vector of portfolio proportions) and 
q is an arbitrary unit cost portfolio, then the following decomposition therefore 
holds: 

q = f q + U q = /? qp p + (1 - /3 qp ) Z p + U q (6.4.21) 

where the three components (i.e. the vectors p, z p and u q ) are mutually orthogo- 
nal. 
We can extend this decomposition to cover 

1. portfolio proportions (orthogonal vectors) 

2. portfolio proportions (scalars/components) 

3. returns (uncorrected random variables) 

4. expected returns (numbers) 

Note again the parallel between orthogonal portfolio vectors and uncorrelated 

portfolio returns/payoffs. 

We will now derive the relation: 

E[r n ] - E[r Zp ] = /3 qp (£[f p ] - E[r Zp }) (6.4.22) 



5 Aside: For the frontier portfolio f q to second degree stochastically dominate the arbitrary 
portfolio q, we will need zero conditional expected return on u q , and will have to show that 

Cov [f Uq) f fq ] = =► E[f u Jf t J = 

The normal distribution is the only case where this is true. 
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which may be familiar from earlier courses in financial economics and which is 

quite general and neither requires asset returns to be normally distributed nor any 

assumptions about preferences. 

Since Cov ^uqj^p = Gov f Zp ,f p = 0, taking covariances with f p in (6.4.21) 

gives: 

Cov [f q ,f p ] = Cov [ff q ,f p ] = /3 qp Var[f p ] (6.4.23) 

or 

Thus (5 in (6.4.21) has its usual definition from probability theory, given by (5.2.2). 6 | 
Reversing the roles of p and z p , it can be seen that 

/3 qZp = 1 - /? qp (6.4.25) 

Taking expected returns in (6.4.21) yields again: 

#[r q ] = P^E[r p ] + (1 - /3 qp )£[f Zp ], (6.4.26) 

which can be rearranged to obtain (6.4.22). 

The Global Minimum Variance Portfolio 

Var[f g+Mh ] = g T Vg + 2/x(g T Vh) + // 2 (h T Vh) (6.4.27) 

which has its minimum at 

g T Vh 

The latter expression reduces to A/C and the minimum value of the variance is 
1/C. The global minimum variance portfolio is denoted mvp. 

t- / K T Vh \ 

Cov[f h ,f MV p] = h V (g-|rVh h ) (6A29) 

T g T Vh.h T Vh 

= h v g- , TU , = (6.4.30) 

h Vh 

i.e. the returns on the portfolio with weights h and the minimum variance portfolio 
are uncorrelated. 



6 Assign some problems involving the construction of portfolio proportions for various desired 
/3s. Also problems working from prices for state contingent claims to returns on assets and port- 
folios in both single period and multi -period worlds. 
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Further, if p is any portfolio, the mvp is the minimum variance combination of 
itself and p, i.e. a = solves: 

1 
mm -Var[r ap+ (i_ a )MVpJ (6.4.31) 

which has necessary and sufficient first order condition: 

aVar[f p ] + (1 — 2a)Cov [f p ,f mvp] — (1 — a)Var[f M vp] = (6.4.32) 

Hence, setting o = 0: 

Cov [f p , f MV p] - Var[f MV p] = (6.4.33) 

and the covariance of any portfolio with mvp is 1/C. 

6.4.2 The portfolio frontier in mean-variance space: 
risky assets only 

The portfolio frontier in mean- variance and in mean-standard deviation space| 

We now move on to consider the mean-variance relationship along the portfolio 

frontier. 

The mean, \x, and variance, a 2 , of the rate of return associated with each point on 

the frontier are related by the quadratic equation: 

{a 2 - Var[w MVP f]) = <£(// - E[w MVP r]) 2 , (6.4.34) 

where the shape parameter = C/D represents the variance of the (gross) return 
on the hedge portfolio, h. The two-moment frontier is generally presented as the 
graph in mean- variance space of this parabola, showing the most desirable distri- 
butions attainable, but the frontier can also be thought of as a plane in portfolio 
space or as a line in portfolio weight space. The latter interpretations are far more 
useful when it comes to extending the analysis to higher moments. 
The equations of the frontier in mean- variance and mean-standard deviation space 
can be derived heuristically using the following stylized diagram illustrating the 
portfolio decomposition. 

Figure 3A goes here. 

Applying Pythagoras' theorem to the triangle with vertices at 0, p and mvp yields: 

( A\ 2 
a 2 = Var[f p ] = Var[f MV p] + [» - ^J Var[f h ] (6.4.35) 
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Recall from the coordinate geometry of conic sections that 

Var[f p ] = Var[f MVP ] + (fi - E[r MWP ]) 2 Var[f h ] (6.4.36) 

or 

is a quadratic equation in \x. 

i.e. the equation of the parabola with vertex at 

Var[f p ] = Var[f MV p] = ^ (6.4.38) 

fi = E[r M wp] = £ (6.4.39) 

Thus in mean-variance space, the frontier is a parabola. 

Figure 3.1 1.2 goes here: indicate position of g on figure. 

Similarly, in mean-standard deviation space, the frontier is a hyperbola. To see 
this, recall that: 

a 2 = Var[f MV p] + [ft - -J Var[f h ] (6.4.40) 

is the equation of the hyperbola with vertex at 



a = sJVax[f MWP \ = d- (6.4.41) 

li = ^ (6-4.42) 

centre at a = 0, \x = A/C and asymptotes as indicated. 

Figure 3.1 1.1 goes here: indicate position of g on figure. 

The other half of the hyperbola (a < 0) has no economic meaning. 

Recall two other types of conic sections: 

Var[f h ] < (impossible) gives a circle with center (1/C, A/C). 

Var [f MV p] = (the presence of a riskless asset) allows the square root to be taken 

on both sides: 

A' 



a = ±^--j VVarffh] (6.4.43) 

i.e. the conic section becomes the pair of lines which are its asymptotes otherwise. 
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Portfolios on which the expected return, fi, exceeds Wf^ vp e are termed efficient, 
since they maximise expected return given variance; other frontier portfolios min- 
imise expected return given variance and are inefficient. 
A frontier portfolio is said to be an efficient portfolio iff 

its expected return exceeds the minimum variance expected return A/C = -E[r M vp]-| 
The set of efficient portfolios in $t N (or efficient frontier) is the half-line emanat- 
ing from mvp in the direction of h, and hence is also a convex set. 
Convex combinations (but not all linear combinations with weights summing to 
1) of efficient portfolios are efficient. 

We now consider zero-covariance (zero-beta) portfolios. In portfolio weight space,| 
can easily construct a frontier portfolio having zero covariance with any given 
frontier portfolio: 

Figure 3B goes here. 

Algebraically, the expected return /x on the zero-covariance frontier portfolio of 
a frontier portfolio with expected return fi solves: 

Cov [f M vp + {fi - E[r MyP })r h , f M vp + {fio ~ E[f M yp\)fh\ = (6.4.44) 

or, since r^ and f MV p are uncorrected: 

Var[f MV p] + (A* - E[f MYP ])(fi - £[f M vp])Var[f h ] = (6.4.45) 

To make this true, we must have 

{fi - E[r MWP ])(ii - E[r MWP }) < (6.4.46) 

or n and fi on opposite sides of E[r M vp] as shown. 

There is a neat trick which allows zero-covariance portfolios to be plotted in mean- 
standard deviation space. 
Implicit differentiation of the fi — a relationship (6.4.35) along the frontier yields: 

da a 

— = — (6.4.47) 

da (fi- £[r MV p])Var[f h ] 

so the tangent at (a, fi) intercepts the fi axis at 

dn a 2 
fi-a-f = ix-- — — (6.4.48) 

do (fi - £[r M vpJ)Var[r h J 

= »--, Jr^Tv FT - (" - E ^V (6A49) 

(li - £[r MVP J)Var[r h J 

= ^[^vp]- ? M- r[r Tv N ( 6A5 °) 

{fi - £[r MV pJ)Var[r h J 
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where we substituted for a 2 from the definition of the frontier. 

A little rearrangement shows that expression (6.4.50) satisfies the equation (6.4.45)| 

defining the return on the zero-covariance portfolio. 

In mean-standard deviation space the picture is like this: 

Figure 3.15.1 goes here. 

To find z p in mean-variance space, note that the line joining (a 2 ,fi) to the mvp 
intercepts the fi axis at: 

V- a —2 — aTT^ T = ^~ a 7 i^F \2m r- i (6.4-51) 

a 2 - Var[r M vpj (A 4 ~ #[r M vp) Var[r h ] 

After cancellation, this is exactly the first expression (6.4.48) for the zero-covariance| 
return we had on the previous page. 

Figure 3.15.2 goes here. 

Alternative derivations 

The treatment of the portfolio frontier with risky assets only concludes with some 
alternative derivations following closely ?. They should probably be omitted alto- 
gether at this stage. 

1. The variance minimisation solution from first principles. 
It can be seen that w is the solution to: 

min L = -w T Vw + A(/x-w T e) +7(1^0 -w T l) (6.4.52) 

{w, , } 2 

which has necessary and sufficient first order conditions: 

dL 

— = Vw - Ae - 7 1 = (6.4.53) 

aw 

dL T 

= / u- w T e = (6.4.54) 

oX 

dL 

= W - W T 1 = (6.4.55) 

The solution can be found by premultiplying the FOC (6.4.13) in turn by 
e T and 1 T and using the constraints yields: 

H = A(e T V- 1 e)+ 7 (e T V- 1 l) (6.4.56) 

1 = A(l T V- 1 e)+ 7 (l T V- 1 l) (6.4.57) 
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The solutions for A and 7 are: 

A = °^ (6.4.58) 

B - An 
7 = D (6.4.59) 

2. Derivation of (6.4.22). 

If we only have frontier portfolio p and interior portfolio q, we get a frontier 
(in fi-a space) entirely within the previous frontier and tangent to it at p. 

The frontiers must have the same slope at p: 
Figure 3C goes here. 

We already saw that the outer frontier has slope , p Zp . 

At the point on the inner frontier with w^ invested in q and (1 — Wq) in p, 

ti = £[f p ] + w q (£[f q - f p ]) (6.4.60) 

a 2 = tOqVar[f q ] 

+2w q (l - «; q )Cov [r p ,f a ] + (1 - «; q ) 2 Var[f p ] (6.4.61) 

Differentiating these w.r.t. w^. 

d/.i 



dlUq 

da 



E[f a - f p ] (6.4.62) 



2<T- — = 2w; q Var[r q ] 

dWq 

+2(1 - 2w q )Cov [f p ,f q ] - 2(1 - w; q )Var[f p j;6.4.63) 

Taking the ratio and setting w q = gives the slope of the inner frontier at 

P- 

d/d _ E[r„ - f p ] 



d a 2Cov[fp,f q ]-2Vax[r p ] 

2^/VarIfp] 



(6.4.64) 



Equating this to the slope of the outer frontier, setting 

Covffrj, f a 1 

A* = ,/r i (6-4.65) 



and rearranging yields: 

E[f a ] ~ E[f Zp ] = /3 qp (£[f p ] - £[f Zp ]) (6.4.66) 
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6.4.3 The portfolio frontier in $l N : 
riskfree and risky assets 

We now consider the mathematics of the portfolio frontier when there is a riskfree 

asset. 

In this case, the frontier portfolio solves: 

min - w Vw (6.4.67) 

w 2 

s.t. w T e+(l - w T l)r/ = y u (6.4.68) 

There is no longer a restriction on portfolio weights, and whatever is not invested 
in the iV risky assets is assumed to be invested in the riskless asset. 
The solution (which can be left as an exercise) is by a similar method to the case 
where all assets were risky: 

w p = V- 1 (e-r f l) ti ~ rf (6.4.69) 

where 

H = (e - lr / ) T V" 1 (e - lr f ) = B - 2Ar f + Cr f 2 > V77 (6.4.70) 

Along the frontier, we have: 



A 4- r f •£ 

~W lf ^ <r 7< 



(6.4.71) 



6.4.4 The portfolio frontier in mean- variance space: 
riskfree and risky assets 

We can now establish the shape of the mean-standard deviation frontier with a 

riskless asset. 

Graphically, in mean-standard deviation space, combining any portfolio p with 

the riskless asset in proportions a and (1 — o) gives a portfolio with expected 

return 

aE[r p ] + (1 - a)r f = 77 + a{E[f p ] - 77) 



and standard deviation of returns aJVar[f p ]. 

i.e. these portfolios trace out the ray in a-fi space emanating from (0, 77) and 

passing through p. 

For each a the highest return attainable is along the ray from 77 which is tangent 

to the frontier generated by the risky assets. 
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On this ray, the riskless asset is held in combination with the tangency portfolio t. 

This only makes sense for 77 < A/C = E[r mvp \. 

Above t, there is a negative weight on the riskless asset — i.e. borrowing. 

Figure 3D goes here. 

Limited borrowing 

Unlimited borrowing as allowed in the preceding analysis is unrealistic. 
Consider what happens 

1. with margin constraints on borrowing: 

Figure 3E goes here. 

The frontier is the envelope of all the finite rays through risky portfolios, 
extending as far as the borrowing constraint allows. 

2. with differential borrowing and lending rates: 

Figure 3F goes here. 

There is a range of expected returns over which a pure risky strategy pro- 
vides minimum variance; 

lower expected returns are achieved by riskless lending; 

and higher expected returns are achieved by riskless borrowing. 

6.5 Market Equilibrium and the Capital Asset Pric- 
ing Model 

6.5.1 Pricing assets and predicting security returns 

Need more waffle here about prediction and the difficulties thereof and the prop- 
erties of equilibrium prices and returns. 

We are looking for assumptions concerning probability distributions that lead to 
useful and parsimonious asset pricing models. The CAPM restrictions are the 
best known. At a very basic level, they can be expressed by saying that every 
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investor has mean-variance preferences. This can be achieved either by restricting 
preferences to be quadratic or the probability distribution of asset returns to be 
normal. CAPM is basically a single-period model, but can be extended by assum- 
ing that return distributions are stable over time. ? and ? have generalised the 
distributional conditions. 

Recall also the limiting behaviour of the variance of the return on an equally 
weighted portfolio as the number of securities included goes to infinity. If se- 
curities are added in such a way that the average of the variance terms and the 
average of the covariance terms are stable, then the portfolio variance approaches 
the average covariance as a lower bound. 

6.5.2 Properties of the market portfolio 

Let 

rrij = weight of security j in the market portfolio m 

Wq(> 0) = individual* 's initial wealth 

Wij = proportion of individual z's initial wealth 
invested in j-th security 
Then total wealth is defined by 

i 

W m0 = J2 W o (6-5.1) 

i=i 

and in equilibrium the relation 

i 

y £,w ij Wi = m j W m0 Vj (6.5.2) 

i=i 

must hold. Dividing by W m0 yields: 

1 W i 

and thus in equilibrium the market portfolio is a convex combination of individual 
portfolios. 

6.5.3 The zero-beta CAPM 

Theorem 6.5.1 (Zero-beta CAPM theorem) If every investor holds a mean-variance^ 
frontier portfolio, then the market portfolio, m, is a mean-variance frontier port- 
folio, and hence, Vq, the CAPM equation 

E [f q ] = (1 - qm ) E [f Zm ] + (3 nm E [f m ] (6.5.4) 

holds. 
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Theorem 6.5.2 All strictly risk-averse investors hold frontier portfolios if and 
only if 

E\fuJ'\] =0 Vq (6.5.5) 

Note the subtle distinction between uncorrelated returns (in the definition of the 
decomposition) and independent returns (in this theorem). They are the same only 
for the normal distribution and related distributions. 

We can view the market portfolio as a frontier portfolio under two fund separation. 
If p is a frontier portfolio, then we showed earlier that for purely mathematical 
reasons in the definition of a frontier portfolio: 

#[r q ] = (1 " Pn P )E[f Zp ] + P np E[r p ] (6.5.6) 

If two fund separation holds, then individuals hold frontier portfolios. 
Since the market portfolio is then on the frontier, it follows that: 

£[f q ] = (1 " Apo W-J + /? qm £[f m ] (6.5.7) 

where 

N 

r m = X^ m /i (6.5.8) 

3=1 



Cov [f q , r m ] 
q> " Var[f 



/3qm = w L r q ' r (6.5.9) 



This implies for any particular security, from the economic assumptions of equi- 
librium and two fund separation: 

E[fj\ = (1 - P jm )E[r z J + (3 jm E[r m \ (6.5.10) 

This relation is the ? Zero-Beta version of the Capital Asset Pricing Model 
(CAPM). 

6.5.4 The traditional CAPM 

Now we add the risk free asset, which will allow us to determine the tangency 

portfolio, t, and to talk about Capital Market Line (return v. standard deviation) 

and the Security Market Line (return v. (5). 

Normally in equilibrium there is zero aggregate supply of the riskfree asset. 

Recommended reading for this part of the course is ?, ?, ? and ?. 

Now we can derive the traditional CAPM. Note that by construction 

r f = E[f Zt ]. (6.5.11) 
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Theorem 6.5.3 (Separation Theorem) The risky asset holdings of all investors 
who hold mean-variance frontier portfolios are in the proportions given by the 
tangency portfolio, t. 

Theorem 6.5.4 (Traditional CAPM Theorem) If every investor holds a mean- 
variance frontier portfolio, then the market portfolio of risky assets, m, is the 
tangency portfolio, t, and hence, Vq, the traditional CAPM equation 

E [r q ] = (1 - /3 qm ) r s + P nm E [f m ] (6.5. 12) 

holds. 

Theorem 6.5.4 is sometimes known as the Sharpe-Lintner Theorem. 
The riskless rate is unique by the No Arbitrage Principle, since otherwise a greedy 
investor would borrow an infinite amount at the lower rate and invest it at the 
higher rate, which is impossible in equilibrium. 

We can also think about what happens the CAPM if there are different riskless 
borrowing and lending rates (see ?). If all individuals face this situation in equi- 
librium, realism demands that both riskless assets are in zero aggregate supply and 
hence that all investors hold risky assets only. 

Note that the No Arbitrage Principle also allows us to rule out correlation matri- 
ces for risky assets which permit the construction of portfolios with zero return 
variance, i.e. synthetic riskless assets. 
Assume that a riskless asset exists, with return rj < i?[f mvp ]. 
If the distributional conditions for two fund separation are satisfied, then the tan- 
gency portfolio, t, must be the market portfolio of risky assets in equilibrium. We 
know then that for any portfolio q (with or without a riskless component): 

£[f q ] -r f = (3 nm (E[f m ] - r f ) (6.5.13) 

This is the traditional Sharpe-Lintner version of the CAPM. 

Figure 4A goes here. 

The next theorem relates to the mean- variance efficiency of the market portfolio. 

Theorem 6.5.5 If 

1. the distributional conditions for two fund separation are satisfied; 

2. risky assets are in strictly positive supply; and 

3. investors have strictly increasing (concave) utility functions 
then the market/tangency portfolio is efficient. 
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Proof By Jensen's inequality and monotonicity, the riskless asset dominates any 
portfolio with 

E[f] < r f (6.5.14) 

for 



E[u(W (l + f))] < u(E[W (l + r)]) 
< u(W (l + r f )) 



(6.5.15) 
(6.5.16) 



Hence the expected returns on all individuals' portfolios exceed 77. 

It follows that the expected return on the market portfolio must exceed 77. 



Q.E.D. 



Now we can calculate the risk premium of the market portfolio. 

CAPM gives a relation between the risk premia on individual assets and the risk 

premium on the market portfolio. 

The risk premium on the market portfolio must adjust in equilibrium to give 

market-clearing. 

In some situations, the risk premium on the market portfolio can be written in 

terms of investors' utility functions. 

Assume there is a riskless asset and returns are multivariate normal (MVN). Recall 

the first order conditions for the canonical portfolio choice problem: 



= Eiu^W^rj-rf)} \/i,j 

= e k (m)] e [fj - 77] + cov [«;( Wi) , fj 

= £[«;(Wi)]£fo - r f ] + E[u'!(Wi)]Cov \W U 



(6.5.17) 
(6.5.18) 
(6.5.19) 

using the definition of covariance and Stein's lemma for MVN distributions. Re- 
arranging: 

E[fj - 77] 



where 



0i 
0* 



= Cov [Wi,?, 

-E[<m)} 



E[u'(W t )] 
is the i-th investor's global absolute risk aversion. Since 



jV 



Wl = W*(l + r f + Y,w lk (h-r f )) 



(6.5.20) 
(6.5.21) 

(6.5.22) 



fc=i 



we have, dropping non-stochastic terms, 



Cov^fj-J =Cov 
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W 0^2 w ikh,fj 



fc=l 



(6.5.23) 



CHAPTER 6. PORTFOLIO THEORY 135 

Hence, 



^- r /! = cov r * 



6,, 



WQ^Wikfk,' 



(6.5.24) 



fc=i 

Summing over i, this gives (since we have J2i W^w^ = W™w m k by market- 
clearing and J2k w mk?k = f m by definition): 

i 
E[fj-r f ](£0i l ) = W m0 Cov [f m ,f,] (6.5.25) 

or E[fj-r f ] = (J2 e i l )~ lW moCov[r m ,f 3 ] (6.5.26) 

i=i 

i.e., in equilibrium, the risk premium on the j-th asset is the product of the aggre- 
gate relative risk aversion of the economy and the covariance between the return 
on the j-th asset and the return on the market. 
Now take the average over j weighted by market portfolio weights: 

i 
E[f m - r f ] = (E -1 WmoVax[(]f m ) (6.5.27) 

i.e., in equilibrium, the risk premium on the market is the product of the aggre- 
gate relative risk aversion of the economy and the variance of the return on the 
market. Equivalently, the return to variability of the market equals the aggregate 
relative risk aversion. 
We conclude with some examples. 

1. Negative exponential utility: 

uAz) = exp{— diz} Oj > (6.5.28) 

CLi 

implies: 

(XX 1 )" 1 = (XX 1 )" 1 > (6.5.29) 

i=i i=i 

and hence the market portfolio is efficient. 

2. Quadratic utility: 

Ui(z) = a,iZ — -^-z 2 di,bi > (6.5.30) 

(E^ 1 )-^^^-^^])) (6.5.31) 

This result can also be derived without assuming MVN and using Stein's 
lemma. 
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Chapter 7 
INVESTMENT ANALYSIS 



7.1 Introduction 

[To be written.] 

7.2 Arbitrage and the Pricing of Derivative Securi- 
ties 

7.2.1 The binomial option pricing model 

This still has to be typed up. It follows very naturally from the stuff in Section 5.4. 

7.2.2 The Black-Scholes option pricing model 

Fischer Black died in 1995. In 1997, Myron Scholes and Robert Merton were 
awarded the Nobel Prize in Economics 'for a new method to determine the value 
of derivatives.' See 

http : //www . nobel . se /announcement- 97 /economy 97 . html 
Black and Scholes considered a world in which there are three assets: a stock, 
whose price, S t , follows the stochastic differential equation: 

dSt = jJiStdt + aStdzt, 

where {z t } t =o * s a Brownian motion process; a bond, whose price, B t , follows the 
differential equation: 

dB t = rB t dt; 

and a call option on the stock with strike price X and maturity date T. 
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They showed how to construct an instantaneously riskless portfolio of stocks and 

options, and hence, assuming that the principle of no arbitrage holds, derived the 

Black-Scholes partial differential equation which must be satisfied by the option 

price. 

The option pays (St — X) + = max {S T — X, 0} at maturity. 

Let the price of the call at time t be C t . Guess that C t = C(S t , t). By Ito's lemma: 

~ (dC dC ~ ; ld 2 C 2 ^\,, M dC ~ 

The no arbitrage principle yields the partial differential equation: 

dC ld 2 C 2q ,dC 

-ob + — — rb — rC = 



dt 2dS 2 dS 

subject to the boundary condition 

C(S,T) = (S-X) + . 

Let r = T — t be the time to maturity. 

Then we claim that the solution to the Black-Scholes equation is: 



C(S, t) = SN (d (S, r)) - Xe~ rT N (d (S, r) - a^r) , 

where N (•) is the cumulative distribution function of the standard normal distri- 
bution and 

In f + (V - \a 2 ) t 
d(S,r) = — K ^ J +o^ (7.2.1) 

0\/T 

In | + (r + \a 2 ) r 

-!—. (7.2.2) 



We can check that this is indeed the solution by calculating the various partial 
derivatives and substituting them in the original equation. 
Note first that 

N(z)= f Z -Le-^ 2 d£ 

J — oo V 27T 

and hence by the fundamental theorem of calculus 

N' (z) = -^e-^ z \ 
v/27 
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which of course is the corresponding probability density function. For the last step 
in this proof, we will need the partials of d (S, r) with respect to S and t, which 
are: 



dd(S,r) dd(S,r) _ { r + \° 2 ) 



dt dr 2a^ 

and 

dd(S,r) 1 



(7.2.3) 
(7.2.4) 



dS Sa-y/r 

Note also that 

N'(d(S,r)-aV^) = e-^ 2T e d{s ' T) ^N'(d(S,T)) (7.2.5) 



e-l^ 



e( r+ ^ 2 ) T \N'{d{S,T)) (7. 2.6) 



S 
^e"N'(d(S,T)). (7.2.7) 



Using these facts and the chain rule, we have: 

dC 
dt 

= SN'(d(S,r)) dd ^ T) -Xe~ rr x 

U> (d (S, r) - aVT) (^fl -^) +rN ( d ^ T ) ~ ^ 



dc 

dS 



(7.2.8) 
-SN' (d (S, t)) -?-= - Xe~ rT rN (d (S, t) - ay/r) (7.2.9) 

2\JT V J 



= SN'(d(S,r))-^p- + N(d(S,r)) 

-Xe^N' (d (S, r) - a^) 9d ^f± (7.2.10) 

= N(d(S,r)) (7.2.11) 

d 2 C 



OS 2 



= N'(d(S,r))^^. (7.2.12) 

ob 

Substituting these expressions in the original partial differential equation yields: 

dC ld 2 C 9 „ dC 



-a 



2 S 2 + ^-rS - rC 



dt 2dS 2 dS 
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N (d (S, r)) (rS - rS) + N' (d (S, r)) (^ + I a ^ 2 ^^ 




+iV (d (S, t) - aVr) [-Xe~ rT r + rXe~ rT ) 
= 0. 

The boundary condition should also be checked. As r — > 0, d (5, r) — > ±oo 
according as 5 > X or 5 < X. In the former case, C (S,T) = S — X; and in the 
latter case, C (S, T) = 0, so the boundary condition is indeed satisfied. 

7.3 Multi-period Investment Problems 

In Section 4.2, it was pointed out that the objects of choice can be differentiated 
not only by their physical characteristics, but also both by the time at which they 
are consumed and by the state of nature in which they are consumed. These 
distinctions were suppressed in the intervening sections but are considered again 
in this section and in Section 5.4 respectively. 

The multi-period model should probably be introduced at the end of Chapter 4 
but could also be left until Chapter 7. For the moment this brief introduction is 
duplicated in both chapters. 

Discrete time multi-period investment problems serve as a stepping stone from 
the single period case to the continuous time case. 

The main point to be gotten across is the derivation of interest rates from equilib- 
rium prices: spot rates, forward rates, term structure, etc. 

This is covered in one of the problems, which illustrates the link between prices 
and interest rates in a multiperiod model. 

7.4 Continuous Time Investment Problems 
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